R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(100.00 + ,100.00 + ,100.00 + ,100.00 + ,100.00 + ,100.39 + ,100.37 + ,100.35 + ,100.33 + ,100.31 + ,100.15 + ,100.26 + ,100.38 + ,100.50 + ,100.61 + ,100.21 + ,100.37 + ,100.52 + ,100.68 + ,100.84 + ,100.03 + ,100.18 + ,100.34 + ,100.49 + ,100.64 + ,99.58 + ,99.78 + ,99.97 + ,100.17 + ,100.36 + ,99.40 + ,99.64 + ,99.88 + ,100.13 + ,100.37 + ,99.77 + ,100.01 + ,100.26 + ,100.50 + ,100.75 + ,100.41 + ,100.67 + ,100.93 + ,101.19 + ,101.45 + ,100.12 + ,100.50 + ,100.88 + ,101.25 + ,101.63 + ,99.83 + ,100.28 + ,100.73 + ,101.18 + ,101.63 + ,99.73 + ,100.24 + ,100.74 + ,101.25 + ,101.75 + ,98.74 + ,99.49 + ,100.25 + ,101.00 + ,101.76 + ,98.44 + ,99.36 + ,100.29 + ,101.22 + ,102.14 + ,98.79 + ,99.68 + ,100.57 + ,101.46 + ,102.35 + ,99.60 + ,100.42 + ,101.24 + ,102.05 + ,102.87 + ,99.82 + ,100.75 + ,101.69 + ,102.62 + ,103.55 + ,99.85 + ,100.87 + ,101.89 + ,102.90 + ,103.92 + ,100.01 + ,101.04 + ,102.07 + ,103.10 + ,104.13 + ,100.28 + ,101.36 + ,102.43 + ,103.51 + ,104.58 + ,100.63 + ,101.57 + ,102.51 + ,103.45 + ,104.39 + ,101.14 + ,101.93 + ,102.71 + ,103.50 + ,104.29 + ,101.51 + ,102.37 + ,103.22 + ,104.08 + ,104.93 + ,102.41 + ,103.10 + ,103.79 + ,104.48 + ,105.17 + ,102.46 + ,103.22 + ,103.99 + ,104.75 + ,105.52 + ,102.09 + ,102.96 + ,103.83 + ,104.70 + ,105.57 + ,101.99 + ,102.77 + ,103.55 + ,104.33 + ,105.11 + ,101.52 + ,102.38 + ,103.24 + ,104.11 + ,104.97 + ,102.44 + ,103.10 + ,103.77 + ,104.43 + ,105.09 + ,103.42 + ,103.90 + ,104.37 + ,104.85 + ,105.33 + ,103.63 + ,104.12 + ,104.61 + ,105.11 + ,105.60 + ,103.28 + ,103.75 + ,104.21 + ,104.68 + ,105.14 + ,103.98 + ,104.37 + ,104.77 + ,105.16 + ,105.56 + ,103.56 + ,103.94 + ,104.33 + ,104.71 + ,105.09 + ,103.42 + ,103.78 + ,104.14 + ,104.51 + ,104.87 + ,103.92 + ,104.15 + ,104.37 + ,104.59 + ,104.81 + ,103.81 + ,104.01 + ,104.20 + ,104.40 + ,104.60 + ,103.09 + ,103.33 + ,103.58 + ,103.83 + ,104.07 + ,102.60 + ,103.05 + ,103.51 + ,103.96 + ,104.41 + ,102.77 + ,103.08 + ,103.39 + ,103.71 + ,104.02 + ,102.60 + ,102.86 + ,103.11 + ,103.37 + ,103.62 + ,102.88 + ,103.08 + ,103.28 + ,103.48 + ,103.68 + ,102.17 + ,102.50 + ,102.83 + ,103.15 + ,103.48 + ,101.85 + ,102.20 + ,102.56 + ,102.91 + ,103.27 + ,101.66 + ,102.14 + ,102.62 + ,103.10 + ,103.58 + ,101.91 + ,102.28 + ,102.66 + ,103.03 + ,103.41 + ,102.13 + ,102.43 + ,102.72 + ,103.02 + ,103.31 + ,102.71 + ,102.82 + ,102.92 + ,103.02 + ,103.13 + ,103.17 + ,103.22 + ,103.26 + ,103.31 + ,103.36 + ,102.89 + ,102.95 + ,103.02 + ,103.08 + ,103.14 + ,102.94 + ,103.14 + ,103.33 + ,103.53 + ,103.73 + ,103.33 + ,103.45 + ,103.57 + ,103.68 + ,103.80 + ,103.75 + ,103.68 + ,103.61 + ,103.54 + ,103.46 + ,104.11 + ,103.98 + ,103.85 + ,103.72 + ,103.60 + ,104.77 + ,104.49 + ,104.22 + ,103.94 + ,103.67 + ,104.62 + ,104.39 + ,104.15 + ,103.92 + ,103.68 + ,105.00 + ,104.76 + ,104.52 + ,104.28 + ,104.04 + ,105.74 + ,105.51 + ,105.27 + ,105.03 + ,104.79 + ,105.94 + ,105.77 + ,105.60 + ,105.43 + ,105.26 + ,106.37 + ,106.18 + ,105.99 + ,105.80 + ,105.62 + ,106.65 + ,106.44 + ,106.23 + ,106.03 + ,105.82 + ,107.08 + ,106.74 + ,106.40 + ,106.05 + ,105.71 + ,106.77 + ,106.51 + ,106.25 + ,106.00 + ,105.74 + ,107.21 + ,106.97 + ,106.74 + ,106.50 + ,106.26 + ,107.34 + ,107.15 + ,106.96 + ,106.78 + ,106.59 + ,107.12 + ,106.93 + ,106.74 + ,106.55 + ,106.36 + ,106.86 + ,106.73 + ,106.59 + ,106.46 + ,106.33 + ,106.92 + ,106.78 + ,106.65 + ,106.51 + ,106.37 + ,106.95 + ,106.75 + ,106.56 + ,106.36 + ,106.17 + ,107.23 + ,106.96 + ,106.69 + ,106.42 + ,106.16 + ,106.94 + ,106.80 + ,106.66 + ,106.51 + ,106.37 + ,106.62 + ,106.51 + ,106.40 + ,106.29 + ,106.18 + ,105.94 + ,105.97 + ,105.99 + ,106.01 + ,106.03 + ,105.91 + ,105.95 + ,105.99 + ,106.03 + ,106.08 + ,106.52 + ,106.45 + ,106.38 + ,106.31 + ,106.24 + ,106.85 + ,106.63 + ,106.41 + ,106.19 + ,105.97 + ,107.22 + ,106.99 + ,106.75 + ,106.52 + ,106.28 + ,107.28 + ,107.09 + ,106.90 + ,106.71 + ,106.52 + ,107.86 + ,107.57 + ,107.29 + ,107.00 + ,106.72 + ,107.68 + ,107.46 + ,107.24 + ,107.02 + ,106.80 + ,108.07 + ,107.82 + ,107.56 + ,107.31 + ,107.06 + ,107.87 + ,107.66 + ,107.45 + ,107.23 + ,107.02 + ,107.65 + ,107.50 + ,107.35 + ,107.19 + ,107.04 + ,108.16 + ,107.89 + ,107.63 + ,107.36 + ,107.09 + ,108.60 + ,108.24 + ,107.88 + ,107.51 + ,107.15 + ,108.92 + ,108.57 + ,108.21 + ,107.86 + ,107.50 + ,109.66 + ,109.22 + ,108.78 + ,108.34 + ,107.90 + ,109.87 + ,109.40 + ,108.94 + ,108.48 + ,108.02 + ,109.54 + ,109.10 + ,108.66 + ,108.22 + ,107.78 + ,109.06 + ,108.72 + ,108.38 + ,108.04 + ,107.70) + ,dim=c(5 + ,90) + ,dimnames=list(c('IND90_UT10' + ,'IND70_UT30' + ,'IND50_UT50' + ,'IND30_UT70' + ,'IND10_UT90') + ,1:90)) > y <- array(NA,dim=c(5,90),dimnames=list(c('IND90_UT10','IND70_UT30','IND50_UT50','IND30_UT70','IND10_UT90'),1:90)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'grey' > ylab = 'value' > xlab = 'variables' > main = 'Notched Boxplots' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2007), Notched Boxplots (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_notchedbox1.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description > z <- as.data.frame(t(y)) > postscript(file="/var/www/html/rcomp/tmp/18joz1193403802.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > (r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=TRUE,col=par1)) $stats [,1] [,2] [,3] [,4] [,5] [1,] 98.440 99.360 99.88 100.000 100.00 [2,] 101.510 102.140 102.62 103.080 103.46 [3,] 103.375 103.715 103.92 104.365 104.84 [4,] 106.860 106.730 106.41 106.310 106.17 [5,] 109.870 109.400 108.94 108.480 108.02 $n [1] 90 90 90 90 90 $conf [,1] [,2] [,3] [,4] [,5] [1,] 102.4840 102.9506 103.2888 103.8271 104.3887 [2,] 104.2660 104.4794 104.5512 104.9029 105.2913 $out numeric(0) $group numeric(0) $names [1] "IND90_UT10" "IND70_UT30" "IND50_UT50" "IND30_UT70" "IND10_UT90" > dev.off() null device 1 > load(file='/var/www/html/rcomp/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,r$stats[j,i]) + } + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/2w4d01193403802.tab") > 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,r$conf[1,i]) + a<-table.element(a,r$stats[3,i]) + a<-table.element(a,r$conf[2,i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/3rqs61193403802.tab") > > system("convert tmp/18joz1193403802.ps tmp/18joz1193403802.png") > > > proc.time() user system elapsed 0.760 0.173 0.804