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(100 + ,100 + ,100 + ,100 + ,108.5365854 + ,108.8888889 + ,105.4082715 + ,95.65217391 + ,135.9756098 + ,129.4444444 + ,122.6935313 + ,100 + ,139.0243902 + ,131.6666667 + ,123.8600212 + ,108.6956522 + ,111.5853659 + ,112.7777778 + ,105.8324496 + ,114.4927536 + ,136.5853659 + ,126.6666667 + ,101.8027572 + ,111.5942029 + ,145.7317073 + ,135 + ,122.9056204 + ,94.20289855 + ,129.8780488 + ,119.4444444 + ,115.6945917 + ,88.4057971 + ,140.2439024 + ,130.5555556 + ,124.3902439 + ,92.75362319 + ,130.4878049 + ,123.3333333 + ,116.4369035 + ,98.55072464 + ,129.2682927 + ,116.1111111 + ,119.6182397 + ,102.8985507 + ,127.4390244 + ,123.3333333 + ,117.3913043 + ,105.7971014 + ,109.1463415 + ,108.3333333 + ,106.0445387 + ,104.3478261 + ,126.2195122 + ,117.2222222 + ,120.1484624 + ,101.4492754 + ,135.3658537 + ,122.2222222 + ,129.7985154 + ,101.4492754 + ,120.7317073 + ,106.6666667 + ,119.300106 + ,101.4492754 + ,107.9268293 + ,98.88888889 + ,110.4984093 + ,105.7971014 + ,119.5121951 + ,106.6666667 + ,98.0911983 + ,108.6956522 + ,126.8292683 + ,110.5555556 + ,124.2841994 + ,104.3478261 + ,120.7317073 + ,108.8888889 + ,115.9066808 + ,111.5942029 + ,113.4146341 + ,100.5555556 + ,112.5132556 + ,115.942029 + ,128.0487805 + ,113.3333333 + ,125.980912 + ,114.4927536 + ,113.4146341 + ,100.5555556 + ,111.6648993 + ,115.942029 + ,115.2439024 + ,103.3333333 + ,112.407211 + ,115.942029 + ,105.4878049 + ,97.77777778 + ,108.1654295 + ,114.4927536 + ,121.9512195 + ,107.7777778 + ,119.7242842 + ,114.4927536 + ,121.3414634 + ,107.2222222 + ,123.5418876 + ,127.5362319 + ,118.902439 + ,103.3333333 + ,121.7391304 + ,117.3913043 + ,98.7804878 + ,93.88888889 + ,106.5747614 + ,117.3913043 + ,107.3170732 + ,91.11111111 + ,90.56203606 + ,118.8405797 + ,120.7317073 + ,105.5555556 + ,91.83457052 + ,115.942029 + ,118.2926829 + ,103.8888889 + ,116.542948 + ,120.2898551 + ,104.8780488 + ,95 + ,106.7868505 + ,123.1884058 + ,128.6585366 + ,119.4444444 + ,122.4814422 + ,124.6376812 + ,108.5365854 + ,98.88888889 + ,106.7868505 + ,126.0869565 + ,106.7073171 + ,100.5555556 + ,104.9840933 + ,126.0869565 + ,109.7560976 + ,105.5555556 + ,108.4835631 + ,123.1884058 + ,116.4634146 + ,105 + ,115.3764581 + ,121.7391304 + ,107.9268293 + ,93.33333333 + ,112.3011665 + ,123.1884058 + ,117.0731707 + ,100.5555556 + ,120.0424178 + ,126.0869565 + ,92.07317073 + ,87.22222222 + ,101.4846235 + ,126.0869565 + ,99.3902439 + ,83.88888889 + ,85.79003181 + ,124.6376812 + ,113.4146341 + ,101.6666667 + ,120.7847296 + ,114.4927536 + ,104.8780488 + ,91.66666667 + ,104.0296925 + ,117.3913043 + ,108.5365854 + ,93.88888889 + ,109.0137858 + ,118.8405797 + ,116.4634146 + ,102.2222222 + ,111.028632 + ,123.1884058 + ,101.2195122 + ,91.11111111 + ,101.6967126 + ,124.6376812 + ,97.56097561 + ,87.22222222 + ,100.3181336 + ,123.1884058 + ,101.8292683 + ,93.88888889 + ,107.7412513 + ,120.2898551 + ,106.097561 + ,92.22222222 + ,110.1802757 + ,118.8405797 + ,109.1463415 + ,92.77777778 + ,116.9671262 + ,126.0869565 + ,108.5365854 + ,92.22222222 + ,120.9968187 + ,134.7826087 + ,84.75609756 + ,80 + ,102.6511135 + ,134.7826087 + ,96.95121951 + ,80.55555556 + ,92.68292683 + ,127.5362319 + ,109.1463415 + ,97.22222222 + ,118.1336161 + ,107.2463768 + ,93.90243902 + ,79.44444444 + ,103.2873807 + ,104.3478261 + ,100 + ,85.55555556 + ,109.1198303 + ,108.6956522 + ,109.1463415 + ,95.55555556 + ,119.5121951 + ,120.2898551 + ,93.29268293 + ,81.11111111 + ,102.8632025 + ,127.5362319 + ,89.02439024 + ,78.88888889 + ,100.8483563 + ,128.9855072 + ,90.85365854 + ,82.77777778 + ,102.757158 + ,124.6376812 + ,91.46341463 + ,78.33333333 + ,104.5599152 + ,121.7391304 + ,101.8292683 + ,86.66666667 + ,118.4517497 + ,121.7391304 + ,99.3902439 + ,81.11111111 + ,116.4369035 + ,121.7391304 + ,71.34146341 + ,66.11111111 + ,95.3340403 + ,121.7391304 + ,92.07317073 + ,75 + ,92.68292683 + ,120.2898551 + ,94.51219512 + ,78.88888889 + ,110.8165429 + ,110.1449275 + ,91.46341463 + ,76.11111111 + ,104.0296925 + ,110.1449275 + ,93.90243902 + ,80 + ,108.9077413 + ,114.4927536 + ,97.56097561 + ,85 + ,111.7709438 + ,115.942029 + ,89.63414634 + ,79.44444444 + ,102.8632025 + ,118.8405797 + ,90.24390244 + ,80.55555556 + ,103.2873807 + ,120.2898551) + ,dim=c(4 + ,72) + ,dimnames=list(c('Import' + ,'Export' + ,'Industriele_productie' + ,'Werkloosheid') + ,1:72)) > y <- array(NA,dim=c(4,72),dimnames=list(c('Import','Export','Industriele_productie','Werkloosheid'),1:72)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par1 = 'blue' > 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/13axo1258396897.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.34146 66.11111 85.79003 88.4058 [2,] 99.08537 86.11111 103.28738 108.6957 [3,] 108.53659 98.88889 110.33934 117.3913 [4,] 120.73171 108.05556 118.87593 123.1884 [5,] 145.73171 135.00000 129.79852 134.7826 $n [1] 72 72 72 72 $conf [,1] [,2] [,3] [,4] [1,] 104.5059 94.80273 107.4367 114.6927 [2,] 112.5672 102.97505 113.2420 120.0899 $out numeric(0) $group numeric(0) $names [1] "Import" "Export" "Industriele_productie" [4] "Werkloosheid" > 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/2gwyj1258396898.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/3jxc61258396898.tab") > > system("convert tmp/13axo1258396897.ps tmp/13axo1258396897.png") > > > proc.time() user system elapsed 0.411 0.190 0.482