R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(3016.7 + ,5585 + ,2756.76 + ,3052.4 + ,5710 + ,2849.27 + ,3099.6 + ,5511 + ,2921.44 + ,3103.3 + ,5403 + ,2981.85 + ,3119.8 + ,5826 + ,3080.58 + ,3093.7 + ,5884 + ,3106.22 + ,3164.9 + ,5965 + ,3119.31 + ,3311.5 + ,5960 + ,3061.26 + ,3410.6 + ,6064 + ,3097.31 + ,3392.6 + ,6046 + ,3161.69 + ,3338.2 + ,5954 + ,3257.16 + ,3285.1 + ,5952 + ,3277.01 + ,3294.8 + ,5960 + ,3295.32 + ,3611.2 + ,5983 + ,3363.99 + ,3611.3 + ,5996 + ,3494.17 + ,3521 + ,6021 + ,3667.03 + ,3519.3 + ,6094 + ,3813.06 + ,3438.3 + ,6202 + ,3917.96 + ,3534.9 + ,6276 + ,3895.51 + ,3705.8 + ,6306 + ,3801.06 + ,3807.6 + ,6342 + ,3570.12 + ,3663 + ,6345 + ,3701.61 + ,3604.5 + ,6328 + ,3862.27 + ,3563.8 + ,6191 + ,3970.1 + ,3511.4 + ,6261 + ,4138.52 + ,3546.5 + ,6253 + ,4199.75 + ,3525.4 + ,6198 + ,4290.89 + ,3529.9 + ,6247 + ,4443.91 + ,3591.6 + ,6293 + ,4502.64 + ,3668.3 + ,6381 + ,4356.98 + ,3728.8 + ,6448 + ,4591.27 + ,3853.6 + ,6470 + ,4696.96 + ,3897.7 + ,6516 + ,4621.4 + ,3640.7 + ,6532 + ,4562.84 + ,3495.5 + ,6526 + ,4202.52 + ,3495.1 + ,6533 + ,4296.49 + ,3268 + ,6498 + ,4435.23 + ,3479.1 + ,6507 + ,4105.18 + ,3417.8 + ,6464 + ,4116.68 + ,3521.3 + ,6453 + ,3844.49 + ,3487.1 + ,6468 + ,3720.98 + ,3529.9 + ,6497 + ,3674.4 + ,3544.3 + ,6808 + ,3857.62 + ,3710.8 + ,6793 + ,3801.06 + ,3641.9 + ,6907 + ,3504.37 + ,3447.1 + ,6792 + ,3032.6 + ,3386.8 + ,6757 + ,3047.03 + ,3438.5 + ,6734 + ,2962.34 + ,3364.3 + ,6654 + ,2197.82 + ,3462.7 + ,6589 + ,2014.45 + ,3291.9 + ,6469 + ,1862.83 + ,3550 + ,6521 + ,1905.41 + ,3611 + ,6448 + ,1810.99 + ,3708.6 + ,6410 + ,1670.07 + ,3771.1 + ,6528 + ,1864.44 + ,4042.7 + ,6445 + ,2052.02 + ,3988.4 + ,6458 + ,2029.6 + ,3851.2 + ,6215 + ,2070.83 + ,3876.7 + ,6167 + ,2293.41) + ,dim=c(3 + ,59) + ,dimnames=list(c('zichtrekeningen' + ,'cunsumptiekrediet' + ,'bel20 ') + ,1:59)) > y <- array(NA,dim=c(3,59),dimnames=list(c('zichtrekeningen','cunsumptiekrediet','bel20 '),1:59)) > 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: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > z <- as.data.frame(t(y)) > postscript(file="/var/www/html/rcomp/tmp/1pcjo1258151982.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] [1,] 3016.7 5403.0 1670.070 [2,] 3389.7 6055.0 2972.095 [3,] 3521.3 6342.0 3504.370 [4,] 3641.3 6502.5 4037.640 [5,] 3988.4 6907.0 4696.960 $n [1] 59 59 59 $conf [,1] [,2] [,3] [1,] 3469.546 6249.95 3285.189 [2,] 3573.054 6434.05 3723.551 $out [1] 4042.7 $group [1] 1 $names [1] "zichtrekeningen" "cunsumptiekrediet" "bel20\r" Warning messages: 1: In axis(side = 1, at = 1:3, labels = c("zichtrekeningen", "cunsumptiekrediet", : font width unknown for character 0xd 2: In axis(side = 1, at = 1:3, labels = c("zichtrekeningen", "cunsumptiekrediet", : font width unknown for character 0xd > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > 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/2lstg1258151982.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/3pyxe1258151982.tab") > > system("convert tmp/1pcjo1258151982.ps tmp/1pcjo1258151982.png") > > > proc.time() user system elapsed 0.391 0.172 0.454