x <- array(list(23 ,NA ,17 ,NA ,23 ,NA ,24 ,NA ,17 ,NA ,20 ,NA ,22 ,NA ,18 ,NA ,20 ,NA ,NA ,20 ,NA ,21 ,NA ,21 ,24 ,NA ,20 ,NA ,24 ,NA ,27 ,NA ,28 ,NA ,22 ,NA ,NA ,28 ,NA ,19 ,NA ,23 ,27 ,NA ,22 ,NA ,20 ,NA ,24 ,NA ,16 ,NA ,25 ,NA ,23 ,NA ,18 ,NA ,23 ,NA ,NA ,24 ,NA ,25 ,NA ,27 ,NA ,27 ,NA ,17 ,NA ,27 ,27 ,NA ,14 ,NA ,22 ,NA ,28 ,NA ,11 ,NA ,24 ,NA ,27 ,NA ,27 ,NA ,25 ,NA ,NA ,23 ,NA ,20 ,NA ,22 ,24 ,NA ,22 ,NA ,28 ,NA ,NA ,28 ,NA ,22 ,NA ,28 ,NA ,27 ,NA ,21 ,NA ,27 ,25 ,NA ,23 ,NA ,25 ,NA ,NA ,19 ,NA ,17 ,NA ,16 ,24 ,NA ,24 ,NA ,28 ,NA ,20 ,NA ,14 ,NA ,21 ,NA ,NA ,28 ,NA ,17 ,NA ,24 ,26 ,NA ,23 ,NA ,27 ,NA ,23 ,NA ,24 ,NA ,14 ,NA ,23 ,NA ,24 ,NA ,14 ,NA ,20 ,NA ,8 ,NA ,27 ,NA ,NA ,11 ,NA ,22 ,NA ,20 ,24 ,NA ,23 ,NA ,21 ,NA ,NA ,25 ,NA ,25 ,NA ,22 ,23 ,NA ,21 ,NA ,21 ,NA ,18 ,NA ,24 ,NA ,12 ,NA ,NA ,20 ,NA ,15 ,NA ,20 ,NA ,20 ,NA ,22 ,NA ,24 ,NA ,24 ,NA ,21 ,NA ,19 ,23 ,NA ,25 ,NA ,28 ,NA ,NA ,25 ,NA ,16 ,NA ,23 ,NA ,28 ,NA ,28 ,NA ,27 ,26 ,NA ,23 ,NA ,22 ,NA ,26 ,NA ,21 ,NA ,27 ,NA ,23 ,NA ,21 ,NA ,26 ,NA ,22 ,NA ,26 ,NA ,22 ,NA ,NA ,24 ,NA ,22 ,NA ,21 ,NA ,21 ,NA ,21 ,NA ,19 ,20 ,NA ,18 ,NA ,24 ,NA ,NA ,22 ,NA ,12 ,NA ,19 ,NA ,20 ,NA ,25 ,NA ,26 ,25 ,NA ,17 ,NA ,22 ,NA ,NA ,20 ,NA ,24 ,NA ,28 ,22 ,NA ,15 ,NA ,21 ,NA ,23 ,NA ,13 ,NA ,23 ,NA ,25 ,NA ,26 ,NA ,28 ,NA ,23 ,NA ,16 ,NA ,10 ,NA ,NA ,23 ,NA ,24 ,NA ,24 ,22 ,NA ,21 ,NA ,21 ,NA ,24 ,NA ,20 ,NA ,21 ,NA ,25 ,NA ,14 ,NA ,24 ,NA ,NA ,21 ,NA ,25 ,NA ,24 ,NA ,12 ,NA ,25 ,NA ,25 ,17 ,NA ,20 ,NA ,25 ,NA ,20 ,NA ,22 ,NA ,23 ,NA ,23 ,NA ,20 ,NA ,21 ,NA ,NA ,23 ,NA ,26 ,NA ,16 ,20 ,NA ,18 ,NA ,17 ,NA ,NA ,28 ,NA ,22 ,NA ,25 ,NA ,24 ,NA ,24 ,NA ,24 ,NA ,24 ,NA ,17 ,NA ,23 ,24 ,NA ,24 ,NA ,25 ,NA ,24 ,NA ,20 ,NA ,23 ,NA ,28 ,NA ,19 ,NA ,28 ,NA ,NA ,25 ,NA ,20 ,NA ,26 ,NA ,21 ,NA ,15 ,NA ,22 ,25 ,NA ,23 ,NA ,19 ,NA ,25 ,NA ,26 ,NA ,26 ,NA ,18 ,NA ,22 ,NA ,18 ,NA ,17 ,NA ,20 ,NA ,18 ,NA ,NA ,26 ,NA ,24 ,NA ,25 ,NA ,28 ,NA ,26 ,NA ,27 ,NA ,21 ,NA ,21 ,NA ,12 ,NA ,27 ,NA ,25 ,NA ,15 ,22 ,NA ,13 ,NA ,21 ,NA ,NA ,NA ,NA ,NA ,NA ,NA ,25 ,NA ,22 ,NA ,22 ,NA ,NA ,22 ,NA ,23 ,NA ,21 ,NA ,23 ,NA ,28 ,NA ,24 ,26 ,NA ,22 ,NA ,27 ,NA ,19 ,NA ,20 ,NA ,22 ,NA ,25 ,NA ,6 ,NA ,28 ,NA ,21 ,NA ,21 ,NA ,26 ,NA ,13 ,NA ,20 ,NA ,10 ,NA ,NA ,24 ,NA ,18 ,NA ,19 ,25 ,NA ,23 ,NA ,22 ,NA ,26 ,NA ,20 ,NA ,21 ,NA ,25 ,NA ,24 ,NA ,24 ,NA ,25 ,NA ,22 ,NA ,25 ,NA ,NA ,22 ,NA ,21 ,NA ,21 ,21 ,NA ,18 ,NA ,20 ,NA ,NA ,23 ,NA ,21 ,NA ,21 ,NA ,25 ,NA ,23 ,NA ,24 ,24 ,NA ,23 ,NA ,23 ,NA ,NA ,21 ,NA ,15 ,NA ,18 ,21 ,NA ,21 ,NA ,24 ,NA ,25 ,NA ,24 ,NA ,24 ,NA ,NA ,22 ,NA ,23 ,NA ,19 ,20 ,NA ,21 ,NA ,20 ,NA ,NA ,20 ,NA ,21 ,NA ,18 ,23 ,NA ,20 ,NA ,20 ,NA ,28 ,NA ,11 ,NA ,27 ,NA ,23 ,NA ,22 ,NA ,23 ,NA ,28 ,NA ,27 ,NA ,26 ,NA ,24 ,NA ,25 ,NA ,23 ,NA ,18 ,NA ,18 ,NA ,17 ,NA ,20 ,NA ,20 ,NA ,21 ,NA ,28 ,NA ,24 ,NA ,25 ,NA ,NA ,21 ,NA ,10 ,NA ,23 ,21 ,NA ,27 ,NA ,27 ,NA ,25 ,NA ,21 ,NA ,24 ,NA ,NA ,19 ,NA ,21 ,NA ,20 ,18 ,NA ,18 ,NA ,27 ,NA ,21 ,NA ,15 ,NA ,21 ,NA ,22 ,NA ,24 ,NA ,24 ,NA ,24 ,NA ,22 ,NA ,21 ,NA ,15 ,NA ,14 ,NA ,15 ,NA ,28 ,NA ,28 ,NA ,25 ,NA ,NA ,26 ,NA ,18 ,NA ,25 ,23 ,NA ,26 ,NA ,22 ,NA ,26 ,NA ,17 ,NA ,24 ,NA ,NA ,20 ,NA ,19 ,NA ,21 ,NA ,22 ,NA ,22 ,NA ,22 ,20 ,NA ,18 ,NA ,23 ,NA ,NA ,23 ,NA ,24 ,NA ,22 ,NA ,22 ,NA ,15 ,NA ,20 ,NA ,24 ,NA ,18 ,NA ,23 ,NA ,23 ,NA ,26 ,NA ,25 ,22 ,NA ,11 ,NA ,23 ,NA ,26 ,NA ,26 ,NA ,22 ,NA ,NA ,23 ,NA ,21 ,NA ,25 ,NA ,27 ,NA ,23 ,NA ,26 ,23 ,NA ,23 ,NA ,22 ,NA ,21 ,NA ,15 ,NA ,24 ,NA ,26 ,NA ,22 ,NA ,24 ,NA ,23 ,NA ,26 ,NA ,25 ,NA ,NA ,21 ,NA ,16 ,NA ,20 ,27 ,NA ,20 ,NA ,26 ,NA ,19 ,NA ,18 ,NA ,21 ,NA ,NA ,23 ,NA ,22 ,NA ,26 ,NA ,25 ,NA ,16 ,NA ,21 ,NA ,23 ,NA ,19 ,NA ,22 ,NA ,22 ,NA ,20 ,NA ,16 ,22 ,NA ,19 ,NA ,26 ,NA ,25 ,NA ,23 ,NA ,28 ,NA ,25 ,NA ,24 ,NA ,18 ,NA ,NA ,28 ,NA ,25 ,NA ,25 ,NA ,28 ,NA ,21 ,NA ,23 ,NA ,20 ,NA ,21 ,NA ,21 ,25 ,NA ,23 ,NA ,20 ,NA ,19 ,NA ,27 ,NA ,25 ,NA ,25 ,NA ,23 ,NA ,22 ,NA ,22 ,NA ,18 ,NA ,21 ,NA ,NA ,18 ,NA ,16 ,NA ,16 ,20 ,NA ,16 ,NA ,18 ,NA) ,dim=c(6 ,162) ,dimnames=list(c('I1_man' ,'I1_vrouw' ,'I2_man' ,'I2_vrouw' ,'I3_man' ,'I3_vrouw') ,1:162)) y <- array(NA,dim=c(6,162),dimnames=list(c('I1_man','I1_vrouw','I2_man','I2_vrouw','I3_man','I3_vrouw'),1:162)) for (i in 1:dim(x)[1]) { for (j in 1:dim(x)[2]) { y[i,j] <- as.numeric(x[i,j]) } } par1 = 'grey' ylab = 'value' xlab = 'variables' main = 'Notched Boxplots' z <- as.data.frame(t(y)) postscript(file="/var/wessaorg/rcomp/tmp/1gs681318932219.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() #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,hyperlink('http://www.xycoon.com/overview.htm','Boxplot statistics','Boxplot overview'),6,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Variable',1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/lower_whisker.htm','lower whisker','definition of lower whisker'),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/lower_hinge.htm','lower hinge','definition of lower hinge'),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/central_tendency.htm','median','definitions about measures of central tendency'),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/upper_hinge.htm','upper hinge','definition of upper hinge'),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/upper_whisker.htm','upper whisker','definition of 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="/var/wessaorg/rcomp/tmp/2d09m1318932219.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="/var/wessaorg/rcomp/tmp/3sriu1318932219.tab") try(system("convert tmp/1gs681318932219.ps tmp/1gs681318932219.png",intern=TRUE))