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