Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
99.2 96.7 101.0 99.0 98.1 100.1 100.0 100.0 100.0 111.6 104.9 90.6 122.2 104.9 86.5 117.6 109.5 89.7 121.1 110.8 90.6 136.0 112.3 82.8 154.2 109.3 70.1 153.6 105.3 65.4 158.5 101.7 61.3 140.6 95.4 62.5 136.2 96.4 63.6 168.0 97.6 52.6 154.3 102.4 59.7 149.0 101.6 59.5 165.5 103.8 61.3
Names of X columns:
X1 X2 X3
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation