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(13132.10 + ,12002.40 + ,4033.40 + ,5056.80 + ,17665.90 + ,15525.50 + ,4643.70 + ,5336.70 + ,16913.00 + ,14247.90 + ,4873.30 + ,5381.60 + ,17318.80 + ,15000.70 + ,4608.80 + ,5590.70 + ,16224.20 + ,14931.40 + ,4733.60 + ,5107.50 + ,15469.60 + ,13333.70 + ,4256.40 + ,5629.00 + ,16557.50 + ,14711.20 + ,4758.30 + ,5700.40 + ,19414.80 + ,17197.30 + ,5356.70 + ,6226.90 + ,17335.00 + ,14985.20 + ,5257.30 + ,5895.20 + ,16525.20 + ,14734.40 + ,5416.90 + ,5925.60 + ,18160.40 + ,15937.80 + ,5813.30 + ,6671.00 + ,15553.80 + ,13028.10 + ,5261.90 + ,5649.20 + ,15262.20 + ,13836.80 + ,4669.20 + ,5962.50 + ,18581.00 + ,16677.50 + ,5855.80 + ,6020.60 + ,17564.10 + ,15130.00 + ,5274.60 + ,5688.80 + ,18948.60 + ,17504.00 + ,5516.70 + ,6466.10 + ,17187.80 + ,16979.90 + ,5819.50 + ,6245.60 + ,17564.80 + ,16012.50 + ,5156.00 + ,6597.80 + ,17668.40 + ,16247.70 + ,5377.30 + ,6027.50 + ,20811.70 + ,19268.20 + ,6386.80 + ,7295.10 + ,17257.80 + ,15533.00 + ,5144.00 + ,5980.90 + ,18984.20 + ,16803.30 + ,6138.50 + ,6711.50 + ,20532.60 + ,17396.10 + ,5567.90 + ,6722.60 + ,17082.30 + ,15155.40 + ,5822.60 + ,6049.20 + ,16894.90 + ,15692.40 + ,5145.50 + ,6448.80 + ,20274.90 + ,18063.70 + ,5706.60 + ,6263.20 + ,20078.60 + ,17568.60 + ,6078.50 + ,7202.90 + ,19900.90 + ,18154.30 + ,6074.50 + ,7334.00 + ,17012.20 + ,15467.40 + ,5577.60 + ,6331.30 + ,19642.90 + ,16956.10 + ,5727.50 + ,7314.20 + ,19024.00 + ,16854.00 + ,6067.10 + ,6560.90 + ,21691.00 + ,19396.40 + ,7069.90 + ,7349.50 + ,18835.90 + ,16457.60 + ,5490.00 + ,6904.40 + ,19873.40 + ,17284.50 + ,5948.30 + ,7482.60 + ,21468.20 + ,18395.30 + ,6177.50 + ,7243.40 + ,19406.80 + ,16938.70 + ,6890.10 + ,7597.70 + ,18385.30 + ,16414.30 + ,5756.20 + ,7113.80 + ,20739.30 + ,18173.40 + ,6528.80 + ,7236.20 + ,22268.30 + ,19919.70 + ,6792.00 + ,8380.70 + ,21569.00 + ,19623.80 + ,6657.40 + ,7839.30 + ,17514.80 + ,17228.10 + ,5753.70 + ,6765.20 + ,21124.70 + ,18730.30 + ,5813.60 + ,8161.90 + ,21251.00 + ,19039.10 + ,5966.50 + ,6192.90 + ,21393.00 + ,19413.30 + ,6147.50 + ,7488.00 + ,22145.20 + ,20013.60 + ,7022.40 + ,8258.30 + ,20310.50 + ,17917.20 + ,6361.00 + ,7842.60 + ,23466.90 + ,21270.30 + ,6717.10 + ,8308.00 + ,21264.60 + ,18766.10 + ,7157.80 + ,8922.40 + ,18388.10 + ,16790.80 + ,5386.20 + ,8077.40 + ,22635.40 + ,19960.60 + ,6965.60 + ,9216.70 + ,22014.30 + ,19586.70 + ,6509.30 + ,8709.50 + ,18422.70 + ,17179.00 + ,5199.30 + ,7098.40 + ,16120.20 + ,14964.90 + ,5200.10 + ,7037.50 + ,16037.70 + ,13918.50 + ,4385.80 + ,6830.20 + ,16410.70 + ,14401.30 + ,4764.10 + ,6092.80 + ,17749.80 + ,15994.60 + ,5300.40 + ,6321.10 + ,16349.80 + ,14521.10 + ,4853.20 + ,5530.30 + ,15662.30 + ,13746.50 + ,4814.00 + ,5374.50 + ,17782.30 + ,15956.00 + ,5391.10 + ,5774.00 + ,16398.90 + ,14332.20 + ,5364.70 + ,5336.10) + ,dim=c(4 + ,60) + ,dimnames=list(c('UITintra-EU' + ,'INintra-EU' + ,'UITextra-EU' + ,'INextra-EU') + ,1:60)) > y <- array(NA,dim=c(4,60),dimnames=list(c('UITintra-EU','INintra-EU','UITextra-EU','INextra-EU'),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/1u1ep1260011875.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,] 13132.10 12002.40 4033.40 5056.80 [2,] 16962.60 14992.95 5177.65 5944.05 [3,] 18386.70 16734.15 5572.75 6579.35 [4,] 20421.55 18109.00 6108.50 7324.10 [5,] 23466.90 21270.30 7157.80 9216.70 $n [1] 60 60 60 60 $conf [,1] [,2] [,3] [,4] [1,] 17681.15 16098.55 5382.878 6297.851 [2,] 19092.25 17369.75 5762.622 6860.849 $out numeric(0) $group numeric(0) $names [1] "UITintra-EU" "INintra-EU" "UITextra-EU" "INextra-EU" > 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/2kcms1260011875.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/307f61260011875.tab") > > system("convert tmp/1u1ep1260011875.ps tmp/1u1ep1260011875.png") > > > proc.time() user system elapsed 0.368 0.155 0.469