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(6539 + ,2605 + ,11078 + ,7752 + ,6699 + ,2682 + ,11157 + ,8089 + ,6962 + ,2755 + ,11203 + ,8279 + ,6981 + ,2760 + ,11425 + ,8284 + ,7024 + ,2735 + ,11368 + ,8298 + ,6940 + ,2659 + ,11118 + ,8235 + ,6774 + ,2654 + ,11041 + ,8043 + ,6671 + ,2670 + ,10837 + ,7908 + ,6965 + ,2785 + ,11018 + ,8282 + ,6969 + ,2845 + ,11035 + ,8416 + ,6822 + ,2723 + ,11220 + ,8069 + ,6878 + ,2746 + ,11486 + ,7888 + ,6691 + ,2767 + ,11553 + ,7729 + ,6837 + ,2940 + ,11626 + ,8013 + ,7018 + ,2977 + ,11861 + ,8099 + ,7167 + ,2993 + ,11951 + ,8185 + ,7076 + ,2892 + ,11935 + ,8023 + ,7171 + ,2824 + ,11717 + ,8120 + ,7093 + ,2771 + ,11841 + ,8115 + ,6971 + ,2686 + ,11701 + ,8024 + ,7142 + ,2738 + ,11842 + ,8312 + ,7047 + ,2723 + ,11818 + ,8204 + ,6999 + ,2731 + ,12088 + ,7976 + ,6650 + ,2632 + ,11876 + ,7355 + ,6475 + ,2606 + ,11839 + ,7075 + ,6437 + ,2605 + ,11582 + ,6985 + ,6639 + ,2646 + ,11848 + ,7209 + ,6422 + ,2627 + ,11741 + ,6954 + ,6272 + ,2535 + ,11455 + ,6766 + ,6232 + ,2456 + ,11154 + ,6742 + ,6003 + ,2404 + ,10922 + ,6551 + ,5673 + ,2319 + ,10556 + ,6298 + ,6050 + ,2519 + ,10923 + ,6866 + ,5977 + ,2504 + ,11001 + ,6780 + ,5796 + ,2382 + ,10702 + ,6128 + ,5752 + ,2394 + ,10805 + ,5941 + ,5609 + ,2381 + ,10831 + ,5773 + ,5839 + ,2501 + ,10925 + ,5982 + ,6069 + ,2532 + ,11220 + ,6352 + ,6006 + ,2515 + ,11188 + ,6366 + ,5809 + ,2429 + ,10752 + ,6377 + ,5797 + ,2389 + ,10771 + ,6295 + ,5502 + ,2261 + ,10274 + ,6037 + ,5568 + ,2272 + ,10261 + ,5999 + ,5864 + ,2439 + ,10600 + ,6452 + ,5764 + ,2373 + ,10596 + ,6388 + ,5615 + ,2327 + ,10556 + ,5994 + ,5615 + ,2364 + ,10716 + ,5847 + ,5681 + ,2388 + ,10958 + ,5838 + ,5915 + ,2553 + ,11273 + ,6097 + ,6334 + ,2663 + ,11806 + ,6397 + ,6494 + ,2694 + ,11922 + ,6431 + ,6620 + ,2679 + ,12127 + ,6528 + ,6578 + ,2611 + ,11985 + ,6443 + ,6495 + ,2580 + ,11805 + ,6350 + ,6538 + ,2627 + ,11901 + ,6443 + ,6737 + ,2732 + ,12217 + ,6855 + ,6651 + ,2707 + ,12344 + ,6795 + ,6530 + ,2633 + ,12314 + ,6398 + ,6563 + ,2683 + ,12447 + ,6220) + ,dim=c(4 + ,60) + ,dimnames=list(c('Voeding-Mannen' + ,'Landbouw' + ,'Hotels&Restaurants' + ,'Voeding-Vrouwen') + ,1:60)) > y <- array(NA,dim=c(4,60),dimnames=list(c('Voeding-Mannen','Landbouw','Hotels&Restaurants','Voeding-Vrouwen'),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/1gw6t1257948815.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,] 5502.0 2261.0 10261.0 5773.0 [2,] 5946.0 2478.5 10924.0 6359.0 [3,] 6538.5 2632.5 11320.5 6825.0 [4,] 6909.0 2731.5 11841.5 8033.5 [5,] 7171.0 2993.0 12447.0 8416.0 $n [1] 60 60 60 60 $conf [,1] [,2] [,3] [,4] [1,] 6342.07 2580.894 11133.35 6483.44 [2,] 6734.93 2684.106 11507.65 7166.56 $out numeric(0) $group numeric(0) $names [1] "Voeding-Mannen" "Landbouw" "Hotels&Restaurants" [4] "Voeding-Vrouwen" > 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/26xav1257948815.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/3phjd1257948815.tab") > > system("convert tmp/1gw6t1257948815.ps tmp/1gw6t1257948815.png") > > > proc.time() user system elapsed 0.397 0.159 0.493