Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
74 71 69 67 73 74 71 69 67 73 74 71 69 67 73 74 71 69 67 73 74 71 69 67 73 74 71 69 67 73 75 72 75 69 73 77 72 77 72 74 78 73 79 72 75 79 77 79 73 77 79 77 81 73 77 80 78 81 73 77 80 78 83 73 78 81 79 83 75 79 81 79 83 77 80 82 79 83 78 81 82 80 84 79 81 83 80 84 80 82 83 80 84 80 82 84 80 84 81 82 85 81 84 81 82 85 81 84 82 82 85 81 85 83 82 85 82 85 83 83 86 82 85 83 83 86 82 86 84 83 86 82 86 84 83 87 82 86 84 83 88 82 86 84 83 88 84 86 84 84 88 84 86 84 84 88 85 86 84 84 88 86 87 84 84 88 86 87 85 84 88 86 87 86 85 89 86 87 86 86 89 86 87 87 86 89 86 87 87 86 89 86 88 87 86 90 86 88 87 87 90 87 88 88 88 91 88 88 88 88 91 88 88 88 88 91 88 88 89 89 91 88 88 89 89 91 89 88 89 89 91 89 88 89 89 91 89 88 89 89 91 90 89 89 89 91 90 89 89 90 92 90 90 90 90 92 91 90 90 90 92 91 90 90 90 93 91 91 90 90 93 91 91 90 91 94 91 91 90 91 94 92 92 90 91 94 92 92 91 92 94 92 92 91 92 95 92 92 91 92 95 92 93 92 92 95 92 93 92 92 95 92 93 94 92 96 92 93 94 92 96 93 93 94 93 97 94 94 94 93 97 94 94 94 93 97 95 94 94 93 97 95 94 95 93 97 95 95 95 93 98 95 95 96 94 98 95 95 96 94 98 96 95 96 94 98 96 95 96 95 98 96 96 96 95 98 97 97 97 96 98 97 97 97 96 98 97 99 98 96 100 98 100 98 97 100 98 101 98 97 100 99 101 100 97 102 100 101 100 98 102 100 102 101 98 102 100 102 102 98 102 102 102 103 98 102 103 102 104 99 103 103 103 105 99 103 103 103 105 100 104 104 103 105 100 104 105 104 105 100 106 105 104 105 100 106 105 105 106 101 106 105 105 106 102 107 105 105 107 104 108 108 105 107 104 112 109 105 107 105 112 111 105 110 107 113 112 108 110 107 114 112 113 111 107 114 112 113 111 107 114 112 113 111 107 114 112 113 111 107 114 112 113 111 107 114 112 113 111 NA
Names of X columns:
WJ10AFS WJ10AVA WJ10ARD WJ10AMA WJ10AKN
box colour
Notched Boxes
TRUE
FALSE
% to trim
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par3 <- '5' par2 <- 'TRUE' par1 <- '' 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