Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
74 77 88 67 74 71 69 67 74 71 69 67 74 71 69 84 74 73 69 72 74 71 69 67 75 72 75 67 77 79 81 73 78 80 77 73 79 79 84 81 79 84 87 89 80 79 84 75 80 86 84 77 81 97 84 87 81 88 86 73 82 71 86 72 82 84 86 83 83 71 93 90 83 77 81 67 84 81 79 80 85 92 85 91 85 82 87 78 85 86 91 73 85 81 94 92 86 86 83 90 86 91 86 90 86 91 86 86 87 87 85 98 88 95 95 81 88 92 92 89 88 91 95 84 88 89 83 84 88 82 83 84 88 80 88 89 88 105 86 87 89 95 91 91 89 82 88 90 89 72 94 94 89 86 105 105 90 96 95 84 90 108 83 79 91 78 69 80 91 86 84 84 91 90 88 86 91 102 85 82 91 92 105 100 91 86 92 100 91 92 87 98 91 80 88 107 91 82 86 90 92 82 88 96 92 103 88 98 92 94 92 84 93 105 90 83 93 95 93 89 94 78 89 87 94 112 104 69 94 91 97 90 94 88 104 103 95 104 84 85 95 86 88 89 95 90 101 105 95 103 102 91 96 111 99 83 96 89 88 94 97 93 90 88 97 89 102 94 97 82 101 107 97 98 89 94 97 94 95 105 98 92 87 90 98 97 87 89 98 105 93 96 98 112 101 94 98 86 87 111 98 96 79 97 98 105 92 92 98 95 94 95 100 88 93 89 100 95 90 95 100 92 103 106 102 81 88 104 102 80 102 88 102 100 105 101 102 99 100 96 102 97 97 97 103 91 93 87 103 103 94 96 104 98 108 84 104 92 105 110 106 109 113 111 106 112 95 102 106 85 113 94 107 90 91 96 108 92 103 110 112 88 103 88 112 112 113 105 113 112 102 106 114 100 113 111 114 105 105 105 114 100 96 111 114 112 105 111 114 96 113 107 114 112 113 111
Names of X columns:
WJ10AFS WJ10AVA WJ10ARD WJ10AMA
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