source('/home/pw/wessanet/cretab') myrfcuid = 'q0893053-2019' x <- array(list(10,10,15,9,14,12,14,14,8,6,19,13,17,12,18,13,10,6,15,12,16,10,12,9,13,12,10,7,14,10,15,11,20,15,9,10,12,12,13,10,16,12,12,11,14,11,15,12,19,15,16,12,16,11,14,9,14,11,14,11,13,9,18,15,15,12,15,9,15,12,13,12,14,9,15,9,14,11,19,12,16,12,16,12,12,12,10,6,11,11,13,12,14,9,11,11,11,9,16,10,9,10,16,9,19,12,13,11,15,9,14,9,15,12,11,6,14,10,15,12,17,11,16,14,13,8,15,9,14,10,15,10,14,10,12,11,12,10,15,12,17,14,13,10,5,8,7,8,10,7,15,11,9,6,9,9,15,12,14,12,11,12,18,9,20,15,20,15,16,13,15,9,14,12,13,9,18,15,14,11,12,11,9,6,19,14,13,11,12,8,14,10,6,10,14,9,11,8,11,9,14,10,12,11,19,14,13,12,14,9,17,13,12,8,16,12,15,14,15,9,15,10,16,12,15,12,12,9,13,9,14,12,17,15,14,12,14,11,14,8,15,11,11,11,11,10,16,12,12,9,12,11,19,15,18,14,16,6,16,9,13,9,11,8,10,7,14,10,14,6,14,9,16,9,10,7,16,11,7,9,16,12,15,9,17,10,11,11,11,7,10,12,13,8,14,13,13,11,13,11,12,12,10,11,15,12,6,3,15,10,15,13,11,10,14,6,14,11,16,12,12,9,15,10,20,15,12,9,9,6,13,9,15,15,19,15,11,9,11,11,17,9,15,11,14,10,15,9,11,6,12,12,15,13,16,12,16,12),dim=c(2,179),dimnames=list(c('Perceived_Ease_of_Use','Perceived_Usefulness'),1:179)) y <- array(NA,dim=c(2,179),dimnames=list(c('Perceived_Ease_of_Use','Perceived_Usefulness'),1:179)) 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 (Mon, 10 Jul 2017 10:23:18 +0200) #Author: root #To cite this work: Wessa P., (2017), Notched Boxplots (v1.0.10) 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/1w6lw1580201459.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/2w5ju1580201459.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/30yck1580201459.tab") try(system("convert /home/pw/wessanet/rcomp/tmp/1w6lw1580201459.ps /home/pw/wessanet/rcomp/tmp/1w6lw1580201459.png",intern=TRUE))