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(35.323 + ,1 + ,35.478 + ,1 + ,4.39 + ,1 + ,41.667 + ,2 + ,22.173 + ,1 + ,28.021 + ,1 + ,18.109 + ,2 + ,13.962 + ,1 + ,40.174 + ,1 + ,16.065 + ,1 + ,18.145 + ,2 + ,18.439 + ,2 + ,10.603 + ,1 + ,34.811 + ,1 + ,69.064 + ,1 + ,51.202 + ,2 + ,14.786 + ,1 + ,33.01 + ,2 + ,81.101 + ,2 + ,89.232 + ,1 + ,21.223 + ,2 + ,15.173 + ,1 + ,241.66 + ,1 + ,26.848 + ,2 + ,8.752 + ,1 + ,60.535 + ,1 + ,60.535 + ,1 + ,26.052 + ,1 + ,49.218 + ,2 + ,30.669 + ,1 + ,18.673 + ,2 + ,86 + ,1 + ,10.632 + ,1 + ,35.802 + ,2 + ,33.974 + ,2 + ,36.972 + ,2 + ,4.928 + ,1 + ,53.976 + ,2 + ,15.467 + ,2 + ,35.723 + ,1 + ,40.424 + ,1 + ,9.706 + ,1 + ,26.532 + ,1 + ,23.843 + ,2 + ,18.062 + ,2 + ,35.681 + ,1 + ,68.125 + ,2 + ,23.937 + ,2 + ,31.479 + ,1 + ,66.659 + ,2 + ,250.234 + ,1 + ,49.469 + ,1 + ,42.951 + ,1 + ,43.402 + ,1 + ,24.112 + ,2 + ,56.95 + ,1 + ,17.313 + ,1 + ,25.658 + ,1 + ,48.172 + ,2 + ,13.891 + ,2 + ,32.048 + ,1 + ,19.797 + ,1 + ,31.317 + ,1 + ,20.966 + ,2 + ,22.708 + ,1 + ,26.81 + ,2 + ,52.004 + ,2 + ,32.354 + ,2 + ,27.128 + ,1 + ,26.529 + ,1 + ,28.392 + ,1 + ,57.393 + ,2 + ,194.731 + ,2 + ,9.415 + ,1 + ,91.076 + ,1 + ,57.751 + ,1 + ,8.236 + ,1 + ,20.407 + ,2 + ,13.681 + ,2 + ,79.659 + ,2 + ,53.48 + ,2 + ,6.906 + ,1 + ,50.202 + ,37.877 + ,1 + ,85.903 + ,2 + ,35.351 + ,2 + ,283.801 + ,1 + ,5.974 + ,1 + ,3.441 + ,1 + ,51.987 + ,1 + ,13.22 + ,1 + ,1.455 + ,2 + ,18.187 + ,1 + ,21.29 + ,1 + ,5.686 + ,1 + ,4.944 + ,1 + ,32.789 + ,2 + ,50.494 + ,1 + ,35.162 + ,2 + ,38.095 + ,2 + ,19.172 + ,1 + ,24.5 + ,2 + ,20.573 + ,1 + ,42.042 + ,1 + ,302.912 + ,2 + ,25.027 + ,1 + ,16.488 + ,2 + ,32.36 + ,1 + ,6.193 + ,1 + ,37.7 + ,1 + ,6.343 + ,1 + ,23.025 + ,1 + ,48.578 + ,1 + ,21.564 + ,1 + ,33.697 + ,1 + ,10.831 + ,2 + ,19.172 + ,1 + ,21.075 + ,1 + ,33.189 + ,2 + ,60.5 + ,1 + ,33.686 + ,1 + ,40.838 + ,1 + ,13.491 + ,1 + ,106.637 + ,1 + ,35.897 + ,1 + ,7.314 + ,2 + ,49.094 + ,1 + ,14.667 + ,1 + ,54.179 + ,2 + ,145.846 + ,2 + ,18.56 + ,1 + ,23.525 + ,2 + ,21.804 + ,2 + ,26.301 + ,2 + ,41.33 + ,2 + ,10.5 + ,1 + ,13.338 + ,1 + ,60.31 + ,2 + ,34.256 + ,2 + ,48.267 + ,1 + ,41.559 + ,1 + ,32.45 + ,1 + ,10.951 + ,1 + ,22.561 + ,2 + ,57.095 + ,1 + ,19.105 + ,1 + ,13.151 + ,2 + ,27.426 + ,2 + ,15.355 + ,2 + ,13.82 + ,2 + ,47.21 + ,1 + ,110.349 + ,1 + ,34.985 + ,1 + ,27.257 + ,2 + ,23.556 + ,2 + ,50.108 + ,2 + ,18.158 + ,1 + ,87.357 + ,1 + ,18.187 + ,1 + ,28.33 + ,1 + ,13.474 + ,2 + ,26.244 + ,1) + ,dim=c(2 + ,162) + ,dimnames=list(c('A' + ,'B') + ,1:162)) > y <- array(NA,dim=c(2,162),dimnames=list(c('A','B'),1:162)) > 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 = 0.47, df = 322, p-value = 0.6387 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -6.208945 10.106587 sample estimates: mean of x mean of y 21.82813 19.87931 > (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.883, num df = 161, denom df = 161, p-value = 0.4307 alternative hypothesis: true ratio of variances is not equal to 1 95 percent confidence interval: 0.6476081 1.2039296 sample estimates: ratio of variances 0.8829918 > (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 = 0.47, df = 320.761, p-value = 0.6387 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -6.209064 10.106706 sample estimates: mean of x mean of y 21.82813 19.87931 > (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 = 13579.5, p-value = 0.5811 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.0679, p-value = 0.8492 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.4877, p-value < 2.2e-16 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/wessaorg/rcomp/tmp/1mt5n1355652439.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(1, 1, 5.917, 31.479, 69.064, 1, 1, 2, 26.301, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/21hen1355652439.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/3bzr41355652439.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/4ofon1355652439.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/5yhns1355652439.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/6knaq1355652439.tab") > > try(system("convert tmp/1mt5n1355652439.ps tmp/1mt5n1355652439.png",intern=TRUE)) character(0) > try(system("convert tmp/21hen1355652439.ps tmp/21hen1355652439.png",intern=TRUE)) character(0) > try(system("convert tmp/3bzr41355652439.ps tmp/3bzr41355652439.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.834 0.302 2.115