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(1258 + ,476049 + ,280 + ,47 + ,1199 + ,474605 + ,557 + ,55 + ,1158 + ,470439 + ,831 + ,55 + ,1427 + ,461251 + ,1081 + ,51 + ,934 + ,454724 + ,1318 + ,55 + ,709 + ,455626 + ,1578 + ,58 + ,1186 + ,516847 + ,1859 + ,49 + ,986 + ,525192 + ,2141 + ,52 + ,1033 + ,522975 + ,2428 + ,52 + ,1257 + ,518585 + ,2715 + ,57 + ,1105 + ,509239 + ,3004 + ,61 + ,1179 + ,512238 + ,3309 + ,61 + ,1092 + ,519164 + ,269 + ,32 + ,1092 + ,517009 + ,537 + ,43 + ,1087 + ,509933 + ,813 + ,54 + ,2028 + ,509127 + ,1068 + ,81 + ,2039 + ,500857 + ,1411 + ,72 + ,2010 + ,506971 + ,1675 + ,69 + ,754 + ,569323 + ,1958 + ,30 + ,760 + ,579714 + ,2242 + ,53 + ,715 + ,577992 + ,2524 + ,57 + ,855 + ,565464 + ,2836 + ,64 + ,971 + ,547344 + ,3143 + ,71 + ,815 + ,554788 + ,3522 + ,65 + ,915 + ,562325 + ,285 + ,59 + ,843 + ,560854 + ,574 + ,68 + ,761 + ,555332 + ,865 + ,70 + ,1858 + ,543599 + ,1147 + ,75 + ,2968 + ,536662 + ,1516 + ,74 + ,4061 + ,542722 + ,1789 + ,68 + ,3661 + ,593530 + ,2087 + ,62 + ,3269 + ,610763 + ,2372 + ,56 + ,2857 + ,612613 + ,2669 + ,55 + ,2568 + ,611324 + ,2966 + ,74 + ,2274 + ,594167 + ,3270 + ,58 + ,1987 + ,595454 + ,3652 + ,72 + ,683 + ,590865 + ,329 + ,63 + ,381 + ,589379 + ,658 + ,68 + ,71 + ,584428 + ,988 + ,55 + ,1772 + ,573100 + ,1303 + ,71 + ,3485 + ,567456 + ,1603 + ,79 + ,5181 + ,569028 + ,1929 + ,74 + ,4479 + ,620735 + ,2235 + ,75 + ,3782 + ,628884 + ,2544 + ,69 + ,3067 + ,628232 + ,2872 + ,63 + ,2489 + ,612117 + ,3198 + ,88 + ,1903 + ,595404 + ,3544 + ,68 + ,1330 + ,597141 + ,3903 + ,89 + ,736 + ,593408 + ,332 + ,63 + ,483 + ,590072 + ,665 + ,74 + ,242 + ,579799 + ,1001 + ,63 + ,1334 + ,574205 + ,1329 + ,79 + ,2423 + ,572775 + ,1639 + ,80 + ,3523 + ,572942 + ,1975 + ,77 + ,2986 + ,619567 + ,2304 + ,96 + ,2462 + ,625809 + ,2640 + ,88 + ,1908 + ,619916 + ,2992 + ,56 + ,1575 + ,587625 + ,3330 + ,74 + ,1237 + ,565742 + ,3690 + ,82 + ,904 + ,557274 + ,4063 + ,69) + ,dim=c(4 + ,60) + ,dimnames=list(c('bel.' + ,'werkl.' + ,'ink.' + ,'verz.') + ,1:60)) > y <- array(NA,dim=c(4,60),dimnames=list(c('bel.','werkl.','ink.','verz.'),1:60)) > 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/1mysa1258127023.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] [1,] 71.0 454724.0 269.0 30.0 [2,] 909.5 518874.5 1034.5 55.5 [3,] 1257.5 568242.0 1943.5 64.5 [4,] 2442.5 593469.0 2854.0 74.0 [5,] 4479.0 628884.0 4063.0 96.0 $n [1] 60 60 60 60 $conf [,1] [,2] [,3] [,4] [1,] 944.803 553026.4 1572.364 60.72642 [2,] 1570.197 583457.6 2314.636 68.27358 $out [1] 5181 $group [1] 1 $names [1] "bel." "werkl." "ink." "verz." > 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/2o80q1258127023.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/3edpx1258127023.tab") > > system("convert tmp/1mysa1258127023.ps tmp/1mysa1258127023.png") > > > proc.time() user system elapsed 0.375 0.165 0.445