Wessa.net

Send output to:

Data X:
88 95 95 98 92 87 89 95 91 81 97 84 88 92 92 95 104 84 79 79 84 91 78 69 80 79 84 91 86 84 85 92 85 94 78 89 86 86 83 82 71 86 107 90 91 97 93 90 78 80 77 91 90 88 74 77 88 95 86 88 97 89 102 112 88 103 80 86 84 98 97 87 102 81 88 100 88 93 97 82 101 92 82 88 114 100 113 88 91 95 96 111 99 91 102 85 114 105 105 91 92 105 114 100 96 89 82 88 103 91 93 83 71 93 74 71 69 88 89 83 88 82 83 98 105 93 74 71 69 86 91 86 93 105 90 85 82 87 91 86 92 75 72 75 97 98 89 88 80 88 79 84 87 106 109 113 102 80 102 98 112 101 98 86 87 114 112 105 95 90 101 98 96 79 82 84 86 97 94 95 94 112 104 84 81 79 88 105 86 83 77 81 98 105 92 106 112 95 94 91 97 91 92 87 89 72 94 102 100 105 103 103 94 96 89 88 89 86 105 112 112 113 74 71 69 74 73 69 102 99 100 106 85 113 91 80 88 94 88 104 98 95 94 85 86 91 108 92 103 90 96 95 74 71 69 81 88 86 104 98 108 90 108 83 95 103 102 85 81 94 114 96 113 114 112 113 93 95 93 104 92 105 92 103 88 100 95 90 92 94 92 77 79 81 86 91 86 87 87 85 91 82 86 102 97 97 100 92 103 113 112 102
Names of X columns:
WJ7ARD WJ7AMA WJ7AWL
box colour


Notched Boxes


% to trim


Chart options


Title:
Label y-axis:
Label x-axis:

R Code
par1 <- as.numeric(par1) #colour par2<- as.logical(par2) # Notches ? par3<-as.numeric(par3) # % trim if(par3>45){par3<-45;warning('trim limited to 45%')} if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')} lotrm<-as.integer(length(y[1,])par3/100)+1 hitrm<-as.integer(length(y[1,])(100-par3)/100) y1<-array(dim=c(dim(y)[1], hitrm-lotrm+1), dimnames=list(dimnames(y)[[1]], 1:(hitrm-lotrm+1) )) for(i in 1:dim(y)[1]){ tmp<-order(y[i,]) y1[i,]<- y[i, tmp[lotrm:hitrm] ] } bitmap(file='test2.png') pairs(t(y)) dev.off() y<-y1 z <- as.data.frame(t(y)) bitmap(file='test1.png') (r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=par2,col=par1)) dev.off() load(file='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,round(r$stats[j,i], digits=2)) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab') if (par2){ 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,round(r$conf[1,i], digits=2)) a<-table.element(a, round(r$stats[3,i], digits=2)) a<-table.element(a,round(r$conf[2,i], digits=2)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable1.tab') } a<-table.start() a<-table.row.start(a) a<-table.element(a,'Boxplot Means',4,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/trimmed_mean.htm','trimmed mean','definition of trimmed mean'),1,TRUE) a<-table.element(a,hyperlink('http://www.xycoon.com/unbiased1.htm','unbiased SD','definition of unbiased SD'),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,round(mean(unlist(z[i]), trim=par3/100, na.rm=TRUE), digits=2)) a<-table.element(a,round(sd(unlist(z[i]), na.rm=TRUE), digits=2)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable2.tab')

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	0 seconds
R Server	Big Analytics Cloud Computing Center
R Framework error message	Warning: there are blank lines in the 'Data X' field. Please, use NA for missing data - blank lines are simply deleted and are NOT treated as missing values.
R Engine error message	Error in list(88, 95, 95, , 98, 92, 87, 89, 95, 91, 81, 97, 84, 88, 92, : argument 4 is empty Calls: array Execution halted