R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing 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(36361044.798 + ,44942851.780 + ,52792117.596 + ,58987856.120 + ,9817684.280 + ,2686941.048 + ,27004830.074 + ,22200509.400 + ,13875207.231 + ,50632268.368 + ,40847986.206 + ,17509460.320 + ,57657199.379 + ,51114646.772 + ,33135196.435 + ,18875537.457 + ,7039635.987 + ,2777207.566 + ,8579609.420 + ,16238483.496 + ,29057116.106 + ,3265808.236 + ,32920447.802 + ,37083270.836 + ,27438431.426 + ,58957430.634 + ,54556562.144 + ,17615360.627 + ,16958873.113 + ,37995371.498 + ,57160689.248 + ,13156891.138 + ,30942147.512 + ,19252361.685 + ,28153012.700 + ,15338592.208 + ,36119785.686 + ,57924863.630 + ,9845637.610 + ,48981537.260 + ,19012326.326 + ,12163566.577 + ,41872234.131 + ,45935987.523 + ,33038681.382 + ,21533386.919 + ,1429002.988 + ,17637213.595 + ,36442251.000 + ,29200765.169 + ,19142553.224 + ,14681808.420 + ,10306108.416 + ,77208823.856 + ,32552845.124 + ,68340235.240 + ,27342318.448 + ,6778449.060 + ,9927903.168 + ,22285123.234 + ,19962132.886 + ,64824345.086 + ,10076578.476 + ,21168774.205 + ,34556371.040 + ,22878477.938 + ,14816666.092 + ,21333969.025 + ,19990601.747 + ,2474425.280 + ,23729529.338 + ,52386211.904 + ,12057368.840 + ,43979764.583 + ,11805589.848 + ,42301947.814 + ,3890720.530 + ,13209990.380 + ,53180992.296 + ,38066037.224 + ,18186921.663 + ,27546948.302 + ,32914142.490 + ,44498674.354 + ,32411086.961 + ,26493078.291 + ,27707941.848 + ,71726744.785 + ,23247527.828 + ,10497804.428 + ,3033894.172 + ,37013819.338 + ,7283147.700 + ,44237870.722 + ,28571004.280 + ,38622407.926) + ,dim=c(2 + ,49) + ,dimnames=list(c('bcc3' + ,'nas') + ,1:49)) > y <- array(NA,dim=c(2,49),dimnames=list(c('bcc3','nas'),1:49)) > 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 = -1.8748, df = 96, p-value = 0.06387 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -13629818.0 389244.4 sample estimates: mean of x mean of y 25614545 32234832 > (v.t <- var.test(z[,par1],z[,par2],conf.level=par3)) F test to compare two variances data: z[, par1] and z[, par2] F = 0.5534, num df = 48, denom df = 48, p-value = 0.04294 alternative hypothesis: true ratio of variances is not equal to 1 95 percent confidence interval: 0.3121757 0.9811254 sample estimates: ratio of variances 0.5534289 > (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 = -1.8748, df = 88.672, p-value = 0.06412 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -13637218.3 396644.8 sample estimates: mean of x mean of y 25614545 32234832 > (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 = 973, p-value = 0.1068 alternative hypothesis: true location shift is not equal to 0 Warning message: In wilcox.test.default(z[, par1], z[, par2], alternative = par4, : cannot compute exact p-value with ties > (ks.t <- ks.test(z[,par1],z[,par2],alternative=par4)) Two-sample Kolmogorov-Smirnov test data: z[, par1] and z[, par2] D = 0.3265, p-value = 0.01077 alternative hypothesis: two-sided Warning message: In ks.test(z[, par1], z[, par2], alternative = par4) : cannot compute exact p-value with 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.2041, p-value = 0.2593 alternative hypothesis: two-sided Warning message: In ks.test(newsam1, newsam2, alternative = par4) : cannot compute exact p-value with ties > mydf <- data.frame(cbind(z[,par1],z[,par2])) > colnames(mydf) <- c('Variable 1','Variable 2') > postscript(file="/var/wessaorg/rcomp/tmp/19sxc1383342747.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > boxplot(mydf, notch=TRUE, ylab='value',main=main) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/23knn1383342747.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/wessaorg/rcomp/tmp/3kgp51383342747.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/44b261383342747.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/wessaorg/rcomp/tmp/5wj401383342747.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/wessaorg/rcomp/tmp/68d2a1383342747.tab") > > try(system("convert tmp/19sxc1383342747.ps tmp/19sxc1383342747.png",intern=TRUE)) character(0) > try(system("convert tmp/23knn1383342747.ps tmp/23knn1383342747.png",intern=TRUE)) character(0) > try(system("convert tmp/3kgp51383342747.ps tmp/3kgp51383342747.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.601 0.609 3.231