R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(60747 + ,57619 + ,60194 + ,57077 + ,58401 + ,56161 + ,58280 + ,55614 + ,59246 + ,57029 + ,60363 + ,57667 + ,59486 + ,57275 + ,60288 + ,57514 + ,61064 + ,58519 + ,61128 + ,58759 + ,60793 + ,58170 + ,61681 + ,58277 + ,63587 + ,61050 + ,65952 + ,63191 + ,66995 + ,64035 + ,66366 + ,63358 + ,65038 + ,62149 + ,64644 + ,61789 + ,63404 + ,60973 + ,63654 + ,61229 + ,61635 + ,60566 + ,62560 + ,61922 + ,63514 + ,62945 + ,65319 + ,64810 + ,67256 + ,66641 + ,67334 + ,67675 + ,67931 + ,66585 + ,66922 + ,66573 + ,66512 + ,66193 + ,66495 + ,65545 + ,65767 + ,64835 + ,68454 + ,67184 + ,71116 + ,69421 + ,73792 + ,72359 + ,75737 + ,74730 + ,77046 + ,75067 + ,76729 + ,74927 + ,76567 + ,74845 + ,78036 + ,75335 + ,80147 + ,78121 + ,82512 + ,79944 + ,84873 + ,82779 + ,83387 + ,81484 + ,81732 + ,79939 + ,81449 + ,80611 + ,80647 + ,79088 + ,80680 + ,79992 + ,79230 + ,77800 + ,77543 + ,76277 + ,75689 + ,75425 + ,74454 + ,74524 + ,73085 + ,72844 + ,71606 + ,70768 + ,71043 + ,70172 + ,67799 + ,67726 + ,67936 + ,68032 + ,67596 + ,67096 + ,67730 + ,68261 + ,66928 + ,67363 + ,66809 + ,67322 + ,55559 + ,57465 + ,55167 + ,57031 + ,51834 + ,54046 + ,45337 + ,47435 + ,40839 + ,43623 + ,44921 + ,48866 + ,47745 + ,53075 + ,48088 + ,53778 + ,45971 + ,51237 + ,43752 + ,50393 + ,42504 + ,49947 + ,42367 + ,50660 + ,41333 + ,50307 + ,41358 + ,52235 + ,40150 + ,51783 + ,38294 + ,51602 + ,34220 + ,47498 + ,31910 + ,45981 + ,29516 + ,43588 + ,27291 + ,42791 + ,25732 + ,41325 + ,23261 + ,38917 + ,20671 + ,36971 + ,18158 + ,33628 + ,16021 + ,31445 + ,13587 + ,28478 + ,8757 + ,19247 + ,4977 + ,12209 + ,4025 + ,10467 + ,3711 + ,10127 + ,3613 + ,10344 + ,3241 + ,10117 + ,2300 + ,8142 + ,1659 + ,6404 + ,1085 + ,4519 + ,807 + ,3482 + ,470 + ,2373 + ,296 + ,1772 + ,190 + ,1178 + ,102 + ,746 + ,89 + ,503 + ,37 + ,293 + ,16 + ,161 + ,4 + ,91 + ,0 + ,55 + ,2 + ,20 + ,2 + ,14 + ,0 + ,8 + ,0 + ,2 + ,0 + ,0 + ,0 + ,1) + ,dim=c(2 + ,111) + ,dimnames=list(c('mannen' + ,'vrouwen') + ,1:111)) > y <- array(NA,dim=c(2,111),dimnames=list(c('mannen','vrouwen'),1:111)) > 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 = 'paired' > par4 = 'two.sided' > par3 = '0.95' > par2 = '2' > par1 = '1' > main = 'Two Samples' > #'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)) Paired t-test data: z[, par1] and z[, par2] t = -3.9039, df = 110, p-value = 0.0001636 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -3038.1445 -992.1979 sample estimates: mean of the differences -2015.171 > (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.1799, num df = 110, denom df = 110, p-value = 0.3870 alternative hypothesis: true ratio of variances is not equal to 1 95 percent confidence interval: 0.8103828 1.7179647 sample estimates: ratio of variances 1.179919 > (r.w <- t.test(z[,par1],z[,par2],var.equal=FALSE,alternative=par4,paired=paired,mu=par6,conf.level=par3)) Paired t-test data: z[, par1] and z[, par2] t = -3.9039, df = 110, p-value = 0.0001636 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -3038.1445 -992.1979 sample estimates: mean of the differences -2015.171 > (w.t <- wilcox.test(z[,par1],z[,par2],alternative=par4,paired=paired,mu=par6,conf.level=par3)) Wilcoxon signed rank test with continuity correction data: z[, par1] and z[, par2] V = 2566.5, p-value = 0.1476 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.1081, p-value = 0.5354 alternative hypothesis: two-sided Warning message: In ks.test(z[, par1], z[, par2], alternative = par4) : cannot compute correct p-values 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.1351, p-value = 0.2629 alternative hypothesis: two-sided Warning message: In ks.test(newsam1, newsam2, alternative = par4) : cannot compute correct p-values with ties > mydf <- data.frame(cbind(z[,par1],z[,par2])) > colnames(mydf) <- c('Variable 1','Variable 2') > postscript(file="/var/www/html/rcomp/tmp/17dea1291737780.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/www/html/rcomp/tmp/27dea1291737780.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/www/html/rcomp/tmp/37dea1291737780.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/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/www/html/rcomp/tmp/4soqu1291737780.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/www/html/rcomp/tmp/5h7661291737780.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/www/html/rcomp/tmp/6vg3e1291737780.tab") > > try(system("convert tmp/17dea1291737780.ps tmp/17dea1291737780.png",intern=TRUE)) character(0) > try(system("convert tmp/27dea1291737780.ps tmp/27dea1291737780.png",intern=TRUE)) character(0) > try(system("convert tmp/37dea1291737780.ps tmp/37dea1291737780.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.687 0.436 1.498