Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
86 36 88 86 56 94 103 48 90 74 32 73 63 44 68 82 39 80 93 34 86 77 41 86 111 50 91 71 39 79 103 62 96 89 52 92 75 37 72 88 50 96 84 41 70 85 55 86 70 41 87 104 56 88 88 39 79 77 52 90 77 46 95 72 44 85 83 41 90 110 50 115 91 50 84 80 44 79 91 52 94 86 54 97 85 44 86 107 52 111 93 37 87 87 52 98 84 50 87 73 36 68 84 50 88 86 52 82 99 55 111 75 31 75 87 36 94 79 49 95 82 42 80 95 37 95 84 41 68 85 30 94 95 52 88 63 30 84 85 44 101 86 66 98 75 48 78 98 43 109 71 57 102 63 46 81 71 54 97 84 48 75 81 48 97 79 62 101 63 58 101 93 58 95 92 62 95 83 46 95 80 34 90 111 66 107 92 52 92 79 55 86 69 55 70 83 57 95 80 56 96 91 55 91 97 56 87 85 54 92 85 55 97 99 46 102 67 52 91 87 32 68 68 44 88 81 46 97 80 59 90 93 46 101 93 46 94 102 54 101 104 66 109 90 56 100 85 59 103 92 57 94 82 52 97 85 48 85 89 44 75 77 41 77 79 50 87 76 48 78 101 48 108 81 59 97 92 34 105 89 46 106 81 54 107 77 55 95 95 54 107 85 59 115 81 44 101 76 54 85 93 52 90 104 66 115 89 44 95 76 57 97 77 39 112 71 60 97 79 45 77 89 41 90 81 50 94 99 39 103 81 43 77 84 48 98 85 37 90 111 58 111 78 46 77 111 43 88 78 44 75 87 34 92 92 30 78 93 50 106 70 39 80 84 37 87 75 55 92 96 41 111 85 39 86 87 36 85 75 43 90 103 50 101 86 55 94 77 43 86 74 60 86 74 48 90 76 30 75 83 43 86 101 39 91 83 52 97 92 39 91 74 39 70 87 56 98 71 59 96 79 46 95 83 57 100 80 50 95 90 54 97 80 50 97 96 60 92 109 59 115 98 41 88 85 48 87 83 59 100 86 60 98 72 56 102 75 51 96
Names of X columns:
MWARM30 MC30VRB WISCRY7V
box colour
Notched Boxes
3
TRUE
FALSE
% to trim
Chart options
Title:
Label y-axis:
Label x-axis:
R Code
par3 <- '0' par2 <- 'TRUE' par1 <- '3' 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