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.1
+ ,12002.4
+ ,4033.4
+ ,5056.8
+ ,17665.9
+ ,15525.5
+ ,4643.7
+ ,5336.7
+ ,16913
+ ,14247.9
+ ,4873.3
+ ,5381.6
+ ,17318.8
+ ,15000.7
+ ,4608.8
+ ,5590.7
+ ,16224.2
+ ,14931.4
+ ,4733.6
+ ,5107.5
+ ,15469.6
+ ,13333.7
+ ,4256.4
+ ,5629
+ ,16557.5
+ ,14711.2
+ ,4758.3
+ ,5700.4
+ ,19414.8
+ ,17197.3
+ ,5356.7
+ ,6226.9
+ ,17335
+ ,14985.2
+ ,5257.3
+ ,5895.2
+ ,16525.2
+ ,14734.4
+ ,5416.9
+ ,5925.6
+ ,18160.4
+ ,15937.8
+ ,5813.3
+ ,6671
+ ,15553.8
+ ,13028.1
+ ,5261.9
+ ,5649.2
+ ,15262.2
+ ,13836.8
+ ,4669.2
+ ,5962.5
+ ,18581
+ ,16677.5
+ ,5855.8
+ ,6020.6
+ ,17564.1
+ ,15130
+ ,5274.6
+ ,5688.8
+ ,18948.6
+ ,17504
+ ,5516.7
+ ,6466.1
+ ,17187.8
+ ,16979.9
+ ,5819.5
+ ,6245.6
+ ,17564.8
+ ,16012.5
+ ,5156
+ ,6597.8
+ ,17668.4
+ ,16247.7
+ ,5377.3
+ ,6027.5
+ ,20811.7
+ ,19268.2
+ ,6386.8
+ ,7295.1
+ ,17257.8
+ ,15533
+ ,5144
+ ,5980.9
+ ,18984.2
+ ,16803.3
+ ,6138.5
+ ,6711.5
+ ,20532.6
+ ,17396.1
+ ,5567.9
+ ,6722.6
+ ,17082.3
+ ,15155.4
+ ,5822.6
+ ,6049.2
+ ,16894.9
+ ,15692.4
+ ,5145.5
+ ,6448.8
+ ,20274.9
+ ,18063.7
+ ,5706.6
+ ,6263.2
+ ,20078.6
+ ,17568.6
+ ,6078.5
+ ,7202.9
+ ,19900.9
+ ,18154.3
+ ,6074.5
+ ,7334
+ ,17012.2
+ ,15467.4
+ ,5577.6
+ ,6331.3
+ ,19642.9
+ ,16956.1
+ ,5727.5
+ ,7314.2
+ ,19024
+ ,16854
+ ,6067.1
+ ,6560.9
+ ,21691
+ ,19396.4
+ ,7069.9
+ ,7349.5
+ ,18835.9
+ ,16457.6
+ ,5490
+ ,6904.4
+ ,19873.4
+ ,17284.5
+ ,5948.3
+ ,7482.6
+ ,21468.2
+ ,18395.3
+ ,6177.5
+ ,7243.4
+ ,19406.8
+ ,16938.7
+ ,6890.1
+ ,7597.7
+ ,18385.3
+ ,16414.3
+ ,5756.2
+ ,7113.8
+ ,20739.3
+ ,18173.4
+ ,6528.8
+ ,7236.2
+ ,22268.3
+ ,19919.7
+ ,6792
+ ,8380.7
+ ,21569
+ ,19623.8
+ ,6657.4
+ ,7839.3
+ ,17514.8
+ ,17228.1
+ ,5753.7
+ ,6765.2
+ ,21124.7
+ ,18730.3
+ ,5813.6
+ ,8161.9
+ ,21251
+ ,19039.1
+ ,5966.5
+ ,6192.9
+ ,21393
+ ,19413.3
+ ,6147.5
+ ,7488
+ ,22145.2
+ ,20013.6
+ ,7022.4
+ ,8258.3
+ ,20310.5
+ ,17917.2
+ ,6361
+ ,7842.6
+ ,23466.9
+ ,21270.3
+ ,6717.1
+ ,8308
+ ,21264.6
+ ,18766.1
+ ,7157.8
+ ,8922.4
+ ,18388.1
+ ,16790.8
+ ,5386.2
+ ,8077.4
+ ,22635.4
+ ,19960.6
+ ,6965.6
+ ,9216.7
+ ,22014.3
+ ,19586.7
+ ,6509.3
+ ,8709.5
+ ,18422.7
+ ,17179
+ ,5199.3
+ ,7098.4
+ ,16120.2
+ ,14964.9
+ ,5200.1
+ ,7037.5
+ ,16037.7
+ ,13918.5
+ ,4385.8
+ ,6830.2
+ ,16410.7
+ ,14401.3
+ ,4764.1
+ ,6092.8
+ ,17749.8
+ ,15994.6
+ ,5300.4
+ ,6321.1
+ ,16349.8
+ ,14521.1
+ ,4853.2
+ ,5530.3
+ ,15662.3
+ ,13746.5
+ ,4814
+ ,5374.5
+ ,17782.3
+ ,15956
+ ,5391.1
+ ,5774
+ ,16398.9
+ ,14332.2
+ ,5364.7
+ ,5336.1)
+ ,dim=c(4
+ ,60)
+ ,dimnames=list(c('uitintraeu'
+ ,'inintraeu'
+ ,'uitextraeu'
+ ,'inextraeu')
+ ,1:60))
>  y <- array(NA,dim=c(4,60),dimnames=list(c('uitintraeu','inintraeu','uitextraeu','inextraeu'),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/11wh21257352934.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] "uitintraeu" "inintraeu"  "uitextraeu" "inextraeu" 

> 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/2ik2v1257352934.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/3cv7v1257352934.tab") 
> 
> system("convert tmp/11wh21257352934.ps tmp/11wh21257352934.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.390   0.169   0.492