R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(1418
+ ,26
+ ,869
+ ,20
+ ,1530
+ ,19
+ ,2172
+ ,19
+ ,901
+ ,20
+ ,463
+ ,25
+ ,3201
+ ,25
+ ,371
+ ,22
+ ,1192
+ ,26
+ ,1583
+ ,22
+ ,1439
+ ,17
+ ,1764
+ ,22
+ ,1495
+ ,19
+ ,1373
+ ,24
+ ,2187
+ ,26
+ ,1491
+ ,21
+ ,4041
+ ,13
+ ,1706
+ ,26
+ ,2152
+ ,20
+ ,1036
+ ,22
+ ,1882
+ ,14
+ ,1929
+ ,21
+ ,2242
+ ,7
+ ,1220
+ ,23
+ ,1289
+ ,17
+ ,2515
+ ,25
+ ,2147
+ ,25
+ ,2352
+ ,19
+ ,1638
+ ,20
+ ,1222
+ ,23
+ ,1812
+ ,22
+ ,1677
+ ,22
+ ,1579
+ ,21
+ ,1731
+ ,15
+ ,807
+ ,20
+ ,2452
+ ,22
+ ,829
+ ,18
+ ,1940
+ ,20
+ ,2662
+ ,28
+ ,186
+ ,22
+ ,1499
+ ,18
+ ,865
+ ,23
+ ,1793
+ ,20
+ ,2527
+ ,25
+ ,2747
+ ,26
+ ,1324
+ ,15
+ ,2702
+ ,17
+ ,1383
+ ,23
+ ,1179
+ ,21
+ ,2099
+ ,13
+ ,4308
+ ,18
+ ,918
+ ,19
+ ,1831
+ ,22
+ ,3373
+ ,16
+ ,1713
+ ,24
+ ,1438
+ ,18
+ ,496
+ ,20
+ ,2253
+ ,24
+ ,744
+ ,14
+ ,1161
+ ,22
+ ,2352
+ ,24
+ ,2144
+ ,18
+ ,4691
+ ,21
+ ,1112
+ ,23
+ ,2694
+ ,17
+ ,1973
+ ,22
+ ,1769
+ ,24
+ ,3148
+ ,21
+ ,2474
+ ,22
+ ,2084
+ ,16
+ ,1954
+ ,21
+ ,1226
+ ,23
+ ,1389
+ ,22
+ ,1496
+ ,24
+ ,2269
+ ,24
+ ,1833
+ ,16
+ ,1268
+ ,16
+ ,1943
+ ,21
+ ,893
+ ,26
+ ,1762
+ ,15
+ ,1403
+ ,25
+ ,1425
+ ,18
+ ,1857
+ ,23
+ ,1840
+ ,20
+ ,1502
+ ,17
+ ,1441
+ ,25
+ ,1420
+ ,24
+ ,1416
+ ,17
+ ,2970
+ ,19
+ ,1317
+ ,20
+ ,1644
+ ,15
+ ,870
+ ,27
+ ,1654
+ ,22
+ ,1054
+ ,23
+ ,937
+ ,16
+ ,3004
+ ,19
+ ,2008
+ ,25
+ ,2547
+ ,19
+ ,1885
+ ,19
+ ,1626
+ ,26
+ ,1468
+ ,21
+ ,2445
+ ,20
+ ,1964
+ ,24
+ ,1381
+ ,22
+ ,1369
+ ,20
+ ,1659
+ ,18
+ ,2888
+ ,18
+ ,1290
+ ,24
+ ,2845
+ ,24
+ ,1982
+ ,22
+ ,1904
+ ,23
+ ,1391
+ ,22
+ ,602
+ ,20
+ ,1743
+ ,18
+ ,1559
+ ,25
+ ,2014
+ ,18
+ ,2143
+ ,16
+ ,2146
+ ,20
+ ,874
+ ,19
+ ,1590
+ ,15
+ ,1590
+ ,19
+ ,1210
+ ,19
+ ,2072
+ ,16
+ ,1281
+ ,17
+ ,1401
+ ,28
+ ,834
+ ,23
+ ,1105
+ ,25
+ ,1272
+ ,20
+ ,1944
+ ,17
+ ,391
+ ,23
+ ,761
+ ,16
+ ,1605
+ ,23
+ ,530
+ ,11
+ ,1988
+ ,18
+ ,1386
+ ,24
+ ,2395
+ ,23
+ ,387
+ ,21
+ ,1742
+ ,16
+ ,620
+ ,24
+ ,449
+ ,23
+ ,800
+ ,18
+ ,1684
+ ,20
+ ,1050
+ ,9
+ ,2699
+ ,24
+ ,1606
+ ,25
+ ,1502
+ ,20
+ ,1204
+ ,21
+ ,1138
+ ,25
+ ,568
+ ,22
+ ,1459
+ ,21
+ ,2158
+ ,21
+ ,1111
+ ,22
+ ,1421
+ ,27
+ ,2833
+ ,24
+ ,1955
+ ,24
+ ,2922
+ ,21
+ ,1002
+ ,18
+ ,1060
+ ,16
+ ,956
+ ,22
+ ,2186
+ ,20
+ ,3604
+ ,18
+ ,1035
+ ,20
+ ,1417
+ ,3261
+ ,1587
+ ,1424
+ ,1701
+ ,1249
+ ,946
+ ,1926
+ ,3352
+ ,1641
+ ,2035
+ ,2312
+ ,1369
+ ,1577
+ ,2201
+ ,961
+ ,1900
+ ,1254
+ ,1335
+ ,1597
+ ,207
+ ,1645
+ ,2429
+ ,151
+ ,474
+ ,141
+ ,1639
+ ,872
+ ,1318
+ ,1018
+ ,1383
+ ,1314
+ ,1335
+ ,1403
+ ,910
+ ,616
+ ,1407
+ ,771
+ ,766
+ ,473
+ ,1376
+ ,1232
+ ,1521
+ ,572
+ ,1059
+ ,1544
+ ,1230
+ ,1206
+ ,1205
+ ,1255
+ ,613
+ ,721
+ ,1109
+ ,740
+ ,1126
+ ,728
+ ,689
+ ,592
+ ,995
+ ,1613
+ ,2048
+ ,705
+ ,301
+ ,1803
+ ,799
+ ,861
+ ,1186
+ ,1451
+ ,628
+ ,1161
+ ,1463
+ ,742
+ ,979
+ ,675
+ ,1241
+ ,676
+ ,1049
+ ,620
+ ,1081
+ ,1688
+ ,736
+ ,617
+ ,812
+ ,1051
+ ,1656
+ ,705
+ ,945
+ ,554
+ ,1597
+ ,982
+ ,222
+ ,1212
+ ,1143
+ ,435
+ ,532
+ ,882
+ ,608
+ ,459
+ ,578
+ ,826
+ ,509
+ ,717
+ ,637
+ ,857
+ ,830
+ ,652
+ ,707
+ ,954
+ ,1461
+ ,672
+ ,778
+ ,1141
+ ,680
+ ,1090
+ ,616
+ ,285
+ ,1145
+ ,733
+ ,888
+ ,849
+ ,1182
+ ,528
+ ,642
+ ,947
+ ,819
+ ,757
+ ,894)
+ ,dim=c(2
+ ,289)
+ ,dimnames=list(c('pageviews'
+ ,'I1')
+ ,1:289))
> y <- array(NA,dim=c(2,289),dimnames=list(c('pageviews','I1'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par6 = '0.0'
> par5 = 'unpaired'
> par4 = 'two.sided'
> par3 = '0.95'
> par2 = '2'
> par1 = '1'
> main = 'Two Samples'
> par6 <- '0.0'
> par5 <- 'unpaired'
> par4 <- 'two.sided'
> par3 <- '0.95'
> par2 <- '2'
> par1 <- '1'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication:
> #Technical description:
> par1 <- as.numeric(par1) #column number of first sample
> par2 <- as.numeric(par2) #column number of second sample
> par3 <- as.numeric(par3) #confidence (= 1 - alpha)
> if (par5 == 'unpaired') paired <- FALSE else paired <- TRUE
> par6 <- as.numeric(par6) #H0
> z <- t(y)
> if (par1 == par2) stop('Please, select two different column numbers')
> if (par1 < 1) stop('Please, select a column number greater than zero for the first sample')
> if (par2 < 1) stop('Please, select a column number greater than zero for the second sample')
> if (par1 > length(z[1,])) stop('The column number for the first sample should be smaller')
> if (par2 > length(z[1,])) stop('The column number for the second sample should be smaller')
> if (par3 <= 0) stop('The confidence level should be larger than zero')
> if (par3 >= 1) stop('The confidence level should be smaller than zero')
> (r.t <- t.test(z[,par1],z[,par2],var.equal=TRUE,alternative=par4,paired=paired,mu=par6,conf.level=par3))
Two Sample t-test
data: z[, par1] and z[, par2]
t = 7.7577, df = 576, p-value = 3.953e-14
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
429.1065 720.0492
sample estimates:
mean of x mean of y
1199.4464 624.8685
> (v.t <- var.test(z[,par1],z[,par2],conf.level=par3))
F test to compare two variances
data: z[, par1] and z[, par2]
F = 1.0871, num df = 288, denom df = 288, p-value = 0.4788
alternative hypothesis: true ratio of variances is not equal to 1
95 percent confidence interval:
0.8625208 1.3702425
sample estimates:
ratio of variances
1.087135
> (r.w <- t.test(z[,par1],z[,par2],var.equal=FALSE,alternative=par4,paired=paired,mu=par6,conf.level=par3))
Welch Two Sample t-test
data: z[, par1] and z[, par2]
t = 7.7577, df = 574.998, p-value = 3.963e-14
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
429.1060 720.0497
sample estimates:
mean of x mean of y
1199.4464 624.8685
> (w.t <- wilcox.test(z[,par1],z[,par2],alternative=par4,paired=paired,mu=par6,conf.level=par3))
Wilcoxon rank sum test with continuity correction
data: z[, par1] and z[, par2]
W = 58054, p-value = 4.735e-16
alternative hypothesis: true location shift is not equal to 0
> (ks.t <- ks.test(z[,par1],z[,par2],alternative=par4))
Two-sample Kolmogorov-Smirnov test
data: z[, par1] and z[, par2]
D = 0.3495, p-value = 8.882e-16
alternative hypothesis: two-sided
Warning message:
In ks.test(z[, par1], z[, par2], alternative = par4) :
p-values will be approximate in the presence of ties
> m1 <- mean(z[,par1],na.rm=T)
> m2 <- mean(z[,par2],na.rm=T)
> mdiff <- m1 - m2
> newsam1 <- z[!is.na(z[,par1]),par1]
> newsam2 <- z[,par2]+mdiff
> newsam2 <- newsam2[!is.na(newsam2)]
> (ks1.t <- ks.test(newsam1,newsam2,alternative=par4))
Two-sample Kolmogorov-Smirnov test
data: newsam1 and newsam2
D = 0.2837, p-value = 1.572e-10
alternative hypothesis: two-sided
Warning message:
In ks.test(newsam1, newsam2, alternative = par4) :
p-values will be approximate in the presence of ties
> mydf <- data.frame(cbind(z[,par1],z[,par2]))
> colnames(mydf) <- c('Variable 1','Variable 2')
> postscript(file="/var/fisher/rcomp/tmp/1qqtn1354642186.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> boxplot(mydf, notch=TRUE, ylab='value',main=main)
Warning message:
In bxp(list(stats = c(7, 474, 1220, 1743, 3604, 7, 20, 25, 1161, :
some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device
1
> postscript(file="/var/fisher/rcomp/tmp/259n11354642186.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(z[,par1],main='Normal QQplot - Variable 1')
> qqline(z[,par1])
> dev.off()
null device
1
> postscript(file="/var/fisher/rcomp/tmp/3xbam1354642186.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(z[,par2],main='Normal QQplot - Variable 2')
> qqline(z[,par2])
> dev.off()
null device
1
>
> #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,paste('Two Sample t-test (',par5,')',sep=''),2,TRUE)
> a<-table.row.end(a)
> if(!paired){
+ a<-table.row.start(a)
+ a<-table.element(a,'Mean of Sample 1',header=TRUE)
+ a<-table.element(a,r.t$estimate[[1]])
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Mean of Sample 2',header=TRUE)
+ a<-table.element(a,r.t$estimate[[2]])
+ a<-table.row.end(a)
+ } else {
+ a<-table.row.start(a)
+ a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE)
+ a<-table.element(a,r.t$estimate)
+ a<-table.row.end(a)
+ }
> a<-table.row.start(a)
> a<-table.element(a,'t-stat',header=TRUE)
> a<-table.element(a,r.t$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'df',header=TRUE)
> a<-table.element(a,r.t$parameter[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,r.t$p.value)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'H0 value',header=TRUE)
> a<-table.element(a,r.t$null.value[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Alternative',header=TRUE)
> a<-table.element(a,r.t$alternative)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI Level',header=TRUE)
> a<-table.element(a,attr(r.t$conf.int,'conf.level'))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI',header=TRUE)
> a<-table.element(a,paste('[',r.t$conf.int[1],',',r.t$conf.int[2],']',sep=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'F-test to compare two variances',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'F-stat',header=TRUE)
> a<-table.element(a,v.t$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'df',header=TRUE)
> a<-table.element(a,v.t$parameter[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,v.t$p.value)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'H0 value',header=TRUE)
> a<-table.element(a,v.t$null.value[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Alternative',header=TRUE)
> a<-table.element(a,v.t$alternative)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI Level',header=TRUE)
> a<-table.element(a,attr(v.t$conf.int,'conf.level'))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI',header=TRUE)
> a<-table.element(a,paste('[',v.t$conf.int[1],',',v.t$conf.int[2],']',sep=''))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/4s2nn1354642186.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,paste('Welch Two Sample t-test (',par5,')',sep=''),2,TRUE)
> a<-table.row.end(a)
> if(!paired){
+ a<-table.row.start(a)
+ a<-table.element(a,'Mean of Sample 1',header=TRUE)
+ a<-table.element(a,r.w$estimate[[1]])
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Mean of Sample 2',header=TRUE)
+ a<-table.element(a,r.w$estimate[[2]])
+ a<-table.row.end(a)
+ } else {
+ a<-table.row.start(a)
+ a<-table.element(a,'Difference: Mean1 - Mean2',header=TRUE)
+ a<-table.element(a,r.w$estimate)
+ a<-table.row.end(a)
+ }
> a<-table.row.start(a)
> a<-table.element(a,'t-stat',header=TRUE)
> a<-table.element(a,r.w$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'df',header=TRUE)
> a<-table.element(a,r.w$parameter[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,r.w$p.value)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'H0 value',header=TRUE)
> a<-table.element(a,r.w$null.value[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Alternative',header=TRUE)
> a<-table.element(a,r.w$alternative)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI Level',header=TRUE)
> a<-table.element(a,attr(r.w$conf.int,'conf.level'))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'CI',header=TRUE)
> a<-table.element(a,paste('[',r.w$conf.int[1],',',r.w$conf.int[2],']',sep=''))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/5wa5i1354642186.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,paste('Wicoxon rank sum test with continuity correction (',par5,')',sep=''),2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'W',header=TRUE)
> a<-table.element(a,w.t$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,w.t$p.value)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'H0 value',header=TRUE)
> a<-table.element(a,w.t$null.value[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Alternative',header=TRUE)
> a<-table.element(a,w.t$alternative)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Kolmogorov-Smirnov Test to compare Distributions of two Samples',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'KS Statistic',header=TRUE)
> a<-table.element(a,ks.t$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,ks.t$p.value)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Kolmogorov-Smirnov Test to compare Distributional Shape of two Samples',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'KS Statistic',header=TRUE)
> a<-table.element(a,ks1.t$statistic[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,ks1.t$p.value)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/fisher/rcomp/tmp/6q6p11354642186.tab")
>
> try(system("convert tmp/1qqtn1354642186.ps tmp/1qqtn1354642186.png",intern=TRUE))
character(0)
> try(system("convert tmp/259n11354642186.ps tmp/259n11354642186.png",intern=TRUE))
character(0)
> try(system("convert tmp/3xbam1354642186.ps tmp/3xbam1354642186.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
1.878 0.562 2.420