Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
95 81 84 87 90 97 91 91 88 84 87 87 92 89 92 84 85 91 84 81 80 69 80 73 84 75 73 84 84 77 85 91 86 89 87 82 83 90 82 86 72 73 91 96 99 90 88 90 77 73 78 88 86 100 88 67 73 88 89 95 102 94 93 103 88 91 84 77 81 87 89 98 88 104 88 93 89 90 101 107 85 88 96 79 113 111 107 95 84 82 99 83 98 85 82 84 105 105 99 105 100 89 96 111 105 88 90 94 93 87 92 93 90 81 69 67 NA 83 84 86 83 84 90 93 96 101 69 67 73 86 90 88 90 83 90 87 78 84 92 100 89 75 67 73 89 94 98 88 89 82 87 89 77 113 111 107 102 88 84 101 94 93 87 111 96 105 111 107 101 105 95 79 97 89 86 83 82 95 105 92 104 69 100 79 80 75 86 87 94 81 67 83 92 92 92 95 102 102 97 90 98 87 98 94 94 94 83 105 101 107 94 96 107 88 94 91 105 105 83 113 105 100 69 84 77 69 72 74 100 96 92 113 94 100 88 107 92 104 103 93 94 95 89 91 73 86 103 110 97 95 84 82 69 67 73 86 73 93 108 84 97 83 79 93 102 91 89 94 92 86 113 107 89 113 111 107 93 89 90 105 110 104 88 98 96 90 95 93 92 84 104 81 73 83 86 86 92 85 98 83 86 90 83 97 97 96 103 106 84 102 106 107
Names of X columns:
WJ10ARD WJ10AMA WJ10AKN
box colour
Notched Boxes
TRUE
TRUE
FALSE
% 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')
Compute
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
Click here to blog (archive) this computation