source('/home/pw/wessanet/cretab')



myrfcuid = ''

x <- array(list(70.8,66.67,69.6,66.33,69.87,64.33,67.47,64,67.6,63.33,67.13,61.33,66.27,64.67,66.73,63,68.07,60.67,67.8,63.67,64.8,60.67,64.6,61.67,64.2,62.33,64.2,60.33,63.67,59.67,61,60.33,59.67,59.33,59.67,58.67,59.8,58.67,60.73,59.33,59.4,57.33,58.07,59.33,57.47,56,70.73,53.67,72.87,58.67,66,49.33,66.07,71.33,66,70.33,66.27,69,64,66,63.67,66,63.73,63.33,63.33,65.33,63.53,64.33,63.53,64,62.87,61.67,59.53,63.67,62.8,64.67,60.8,61.67,59.8,62,56.67,61.33,57.67,63.67,58.4,61.33,55.47,62.33,56.2,59.67,71.53,59.33,68.67,61.67,65.67,58.67,66.73,58,67.33,56.67,66.73,59.67,66.87,58,65.8,57,64.73,57.67,65.47,58.67,63.6,55.33,64.07,56,64.67,55.67,63.73,53.33,62.53,53.67,61.93,51,62.67,47,62.8,4.33,61.33,70,62.6,68.67,59.13,67.67,61.27,66,59.47,65.67,57.87,65.67,59.73,63.67,61.4,63.67,58.8,64,58.33,62,57.47,62,57.13,61.67,55,61.67,51.53,63.33,72.73,61,73,62.33,70.8,60.33,70.07,60.33,71.67,60.67,71.07,57.67,70.67,58.33,70.73,58,70.73,57.33,68.6,56.67,69.6,58,66.47,55.33,67.07,55.67,68.67,54.67,66.93,56.33,65.93,55,68.87,55,66.53,54.67,65.8,54.33,66.6,49,66,48.33,65,49.67,66.8,43.67,65.6,6.33,66,3,65.67,66.67,64.67,67.33,65.07,65.33,64.67,66,65.07,65.67,65.2,66.67,64.87,65.67,63.47,65,62.6,64.67,64.07,66.67,63.73,63.67,64.67,63.33,61.6,63.67,61.6,63.33,60.47,63.67,61.27,63,63,61.67,61.47,61.33,60.87,60.67,61.67,60,62.87,61.67,62.4,61.33,59.73,58.67,60.13,60.33,58.8,59.67,59.6,59.33,58.93,59.67,60.13,61,58.2,61,58.27,60,58.27,60,55.07,58.67,53.87,58.33,52.33,58,47.2,56.33,37.93,54.67,72.73,55.33,70.07,54,70.67,52.67,72.07,44,68.8,65.67,68.8,65,67.47,66.33,66.73,64,66.53,62.33,66,61.33,67.6,63,66,63.67,66,62,66.53,61.33,65.8,64.67,64.27,62.67,64.67,64,64.6,61,64.13,60.67,65.47,59.67,62.93,60.33,63.53,56.67,62.13,56.67,63.87,54.33,64.67,51,63.33,51,63.13,47,62.8,68,62.4,65,62.4,64,62.6,64,61.47,64,62.2,62,63,61,61.8,60,59.73,60,60.33,62,60.13,60,59.53,59,59,61,55.93,60,41.87,60,36.33,58,71.67,58,71.47,60,70.47,58,69.53,59,70.73,56,69.93,54,68.73,51,67.53,47),dim=c(2,189),dimnames=list(c('X1','X2'),1:189))
 y <- array(NA,dim=c(2,189),dimnames=list(c('X1','X2'),1:189))
 for (i in 1:dim(x)[1])
 {
 	for (j in 1:dim(x)[2])
 	{
 		y[i,j] <- as.numeric(x[i,j])
 	}
 }
par8 = '1'
par7 = '0.95'
par6 = '0.0'
par5 = '1'
par4 = '1'
par3 = '2'
par2 = '1'
par1 = '-1" OR 3+229-229-1=0+0+0+1 -- '
ylab = 'Y-Axis Title'
xlab = 'X-Axis Title'
main = 'Boxplots Title'
par8 <- '1'
par7 <- '0.95'
par6 <- '0.0'
par5 <- '1'
par4 <- '1'
par3 <- '2'
par2 <- '1'
par1 <- '-1' OR 3+229-229-1=0+0+0+1 -- '
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Sun, 16 Jul 2017 09:48:03 +0200)
#Author: root
#To cite this work: Ian E. Holliday (2017), T-Tests (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL https://www.wessa.net/rwasp_Tests%20to%20Compare%20Two%20Means.wasp/
#Source of accompanying publication: 
#
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.character(par4)
par5 <- as.character(par5)
par6 <- as.numeric(par6)
par7 <- as.numeric(par7)
par8 <- as.logical(par8)
if ( par5 == 'unpaired') paired <- FALSE else paired <- TRUE
x <- t(y)
if(par8){
postscript(file="/home/pw/wessanet/rcomp/tmp/1dqly1638888518.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
(r<-boxplot(x ,xlab=xlab,ylab=ylab,main=main,notch=FALSE,col=2))
dev.off()
}

if( par4 == 'Wilcoxon-Mann_Whitney'){
a<-table.start()
a <- table.row.start(a)
a <- table.element(a,'Wilcoxon Test',3,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'',1,TRUE)
a <- table.element(a,'Statistic',1,TRUE)
a <- table.element(a,'P-value',1,TRUE)
a <- table.row.end(a)
W <- wilcox.test(x[,par2],x[,par3],alternative=par1, paired = paired)
a<-table.row.start(a)
a<-table.element(a,'Wilcoxon Test',1,TRUE)
a<-table.element(a,W$statistic[[1]])
a<-table.element(a,round(W$p.value, digits=5) )
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/2xx7q1638888518.tab") 
}
if( par4 == 'T-Test')
{
T <- t.test(x[,par2],x[,par3],alternative=par1, paired=paired, mu=par6, conf.level=par7)
a<-table.start()
a <- table.row.start(a)
a <- table.element(a,'T-Test',3,TRUE)
a <- table.row.end(a)
if(paired){
a <- table.row.start(a)
a <- table.element(a,'Difference: Mean1 - Mean2',1,TRUE)
a<-table.element(a,round(T$estimate, digits=5) )
a <- table.row.end(a)
}
if(!paired){
a <- table.row.start(a)
a <- table.element(a,'Mean1',1,TRUE)
a<-table.element(a,round(T$estimate[1], digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Mean2',1,TRUE)
a<-table.element(a,round(T$estimate[2], digits=5) )
a <- table.row.end(a)
}
a <- table.row.start(a)
a <- table.element(a,'T Statistic',1,TRUE)
a<-table.element(a,round(T$statistic, digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'P-value',1,TRUE)
a<-table.element(a,round(T$p.value, digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Lower Confidence Limit',1,TRUE)
a<-table.element(a,round(T$conf.int[1], digits=5) )
a <- table.row.end(a)
a<-table.row.start(a)
a <- table.element(a,'Upper Confidence Limit',1,TRUE)
a<-table.element(a,round(T$conf.int[2], digits=5) )
a <- table.row.end(a)
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/39r271638888518.tab") 
}
a<-table.start()
a <- table.row.start(a)
a <- table.element(a,'Standard Deviations',3,TRUE)
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Variable 1',1,TRUE)
a<-table.element(a,round(sd(x[,par2], na.rm=TRUE), digits=5) )
a <- table.row.end(a)
a <- table.row.start(a)
a <- table.element(a,'Variable 2',1,TRUE)
a<-table.element(a,round(sd(x[,par3], na.rm=TRUE), digits=5) )
a <- table.row.end(a)
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/4ya8x1638888518.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/1dqly1638888518.ps /home/pw/wessanet/rcomp/tmp/1dqly1638888518.png",intern=TRUE))