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



myrfcuid = 'r0625664'

x <- array(list(5315.1,5327.75,5349.45,5346.8,5346.6,5325.25,5340.35,5354.75,5382.85,5392.35,5400.35,5410.8,5444.35,5424,5441.85,5447.6,5454.45,5478.8,5490.5,5500.75,5504.25,5513.65,5523.75,5536.4,5547.65,5562.85,5570.4,5589.7,5621.7,5612.3,5631.7,5652.85,5645.45,5664.1,5675.25,5689.65,5700.8,5711.35,5701.85,5732.5,5714.6,5746.35,5753,5764.1,5767.8,5781.9,5805,5805.2,5835.4,5838.8,5851.1,5854.85,5854.95,5870.9,5873.6,5882.75,5867.7,5879.05,5895.6,5891.5,5954.05,5952.95,5960.15,5942.6,5957.55,5949.15,5940.5,5940.1,5926.2,5926.8,5915.3,5912.05,5897,5887.75,5882.6,5905.45,5872,5881.95,5878.4,5874.2,5896.4,5890,5888.5,5873.3,5898.9,5887.65,5907.2,5921.3,5918.75,5920.95,5935.65,5941.3,5936,5931.4,5943.8,5949.85,5953.75,5963.75,5977.1,5973.7,6005.75,6014.5,6023.35,6042.8,6027.7,6041.15,6058.45,6073.2,6096.1,6103.3,6101.55,6115.15,6146.25,6134.3,6136.65,6168.05,6182.8,6204.3,6220.85,6229.75),dim=c(1,120),dimnames=list(c('V1'),1:120))
 y <- array(NA,dim=c(1,120),dimnames=list(c('V1'),1:120))
 for (i in 1:dim(x)[1])
 {
 	for (j in 1:dim(x)[2])
 	{
 		y[i,j] <- as.numeric(x[i,j])
 	}
 }
par2 = 'no'
par1 = 'grey'
ylab = 'value'
xlab = 'variables'
main = 'Notched Boxplots'
par2 <- 'no'
par1 <- 'grey'
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Thu, 22 Sep 2016 14:56:16 +0200)
#Author: root
#To cite this work: Wessa P., (2016), Notched Boxplots (v1.0.9) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_notchedbox1.wasp/
#Source of accompanying publication: Office for Research, Development, and Education
#
if(par2=='yes') {
z <- na.omit(as.data.frame(t(y)))
} else {
z <- as.data.frame(t(y))
}
postscript(file="/home/pw/wessanet/rcomp/tmp/164s61481973405.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
(r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=TRUE,col=par1))
dev.off()

a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot statistics',6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,'lower whisker',1,TRUE)
a<-table.element(a,'lower hinge',1,TRUE)
a<-table.element(a,'median',1,TRUE)
a<-table.element(a,'upper hinge',1,TRUE)
a<-table.element(a,'upper whisker',1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
for (j in 1:5)
{
a<-table.element(a,r$stats[j,i])
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/2vije1481973405.tab") 
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Boxplot Notches',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Variable',1,TRUE)
a<-table.element(a,'lower bound',1,TRUE)
a<-table.element(a,'median',1,TRUE)
a<-table.element(a,'upper bound',1,TRUE)
a<-table.row.end(a)
for (i in 1:length(y[,1]))
{
a<-table.row.start(a)
a<-table.element(a,dimnames(t(x))[[2]][i],1,TRUE)
a<-table.element(a,r$conf[1,i])
a<-table.element(a,r$stats[3,i])
a<-table.element(a,r$conf[2,i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file="/home/pw/wessanet/rcomp/tmp/3ofwa1481973405.tab") 

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