R version 2.7.2 (2008-08-25)
Copyright (C) 2008 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(110.40
+ ,109.20
+ ,99.90
+ ,72.50
+ ,96.40
+ ,88.60
+ ,99.80
+ ,59.40
+ ,101.90
+ ,94.30
+ ,99.80
+ ,85.70
+ ,106.20
+ ,98.30
+ ,100.30
+ ,88.20
+ ,81.00
+ ,86.40
+ ,99.90
+ ,62.80
+ ,94.70
+ ,80.60
+ ,99.90
+ ,87.00
+ ,101.00
+ ,104.10
+ ,100.00
+ ,79.20
+ ,109.40
+ ,108.20
+ ,100.10
+ ,112.00
+ ,102.30
+ ,93.40
+ ,100.10
+ ,79.20
+ ,90.70
+ ,71.90
+ ,100.20
+ ,132.10
+ ,96.20
+ ,94.10
+ ,100.30
+ ,40.10
+ ,96.10
+ ,94.90
+ ,100.60
+ ,69.00
+ ,106.00
+ ,96.40
+ ,100.00
+ ,59.40
+ ,103.10
+ ,91.10
+ ,100.10
+ ,73.80
+ ,102.00
+ ,84.40
+ ,100.20
+ ,57.40
+ ,104.70
+ ,86.40
+ ,100.00
+ ,81.10
+ ,86.00
+ ,88.00
+ ,100.10
+ ,46.60
+ ,92.10
+ ,75.10
+ ,100.10
+ ,41.40
+ ,106.90
+ ,109.70
+ ,100.10
+ ,71.20
+ ,112.60
+ ,103.00
+ ,100.50
+ ,67.90
+ ,101.70
+ ,82.10
+ ,100.50
+ ,72.00
+ ,92.00
+ ,68.00
+ ,100.50
+ ,145.50
+ ,97.40
+ ,96.40
+ ,96.30
+ ,39.70
+ ,97.00
+ ,94.30
+ ,96.30
+ ,51.90
+ ,105.40
+ ,90.00
+ ,96.80
+ ,73.70
+ ,102.70
+ ,88.00
+ ,96.80
+ ,70.90
+ ,98.10
+ ,76.10
+ ,96.90
+ ,60.80
+ ,104.50
+ ,82.50
+ ,96.80
+ ,61.00
+ ,87.40
+ ,81.40
+ ,96.80
+ ,54.50
+ ,89.90
+ ,66.50
+ ,96.80
+ ,39.10
+ ,109.80
+ ,97.20
+ ,96.80
+ ,66.60
+ ,111.70
+ ,94.10
+ ,97.00
+ ,58.50
+ ,98.60
+ ,80.70
+ ,97.00
+ ,59.80
+ ,96.90
+ ,70.50
+ ,97.00
+ ,80.90
+ ,95.10
+ ,87.80
+ ,96.80
+ ,37.30
+ ,97.00
+ ,89.50
+ ,96.90
+ ,44.60
+ ,112.70
+ ,99.60
+ ,97.20
+ ,48.70
+ ,102.90
+ ,84.20
+ ,97.30
+ ,54.00
+ ,97.40
+ ,75.10
+ ,97.30
+ ,49.50
+ ,111.40
+ ,92.00
+ ,97.20
+ ,61.60
+ ,87.40
+ ,80.80
+ ,97.30
+ ,35.00
+ ,96.80
+ ,73.10
+ ,97.30
+ ,35.70
+ ,114.10
+ ,99.80
+ ,97.30
+ ,51.30
+ ,110.30
+ ,90.00
+ ,97.30
+ ,49.00
+ ,103.90
+ ,83.10
+ ,97.30
+ ,41.50
+ ,101.60
+ ,72.40
+ ,97.30
+ ,72.50
+ ,94.60
+ ,78.80
+ ,98.10
+ ,42.10
+ ,95.90
+ ,87.30
+ ,96.80
+ ,44.10
+ ,104.70
+ ,91.00
+ ,96.80
+ ,45.10
+ ,102.80
+ ,80.10
+ ,96.80
+ ,50.30
+ ,98.10
+ ,73.60
+ ,96.80
+ ,40.90
+ ,113.90
+ ,86.40
+ ,96.80
+ ,47.20
+ ,80.90
+ ,74.50
+ ,96.80
+ ,36.90
+ ,95.70
+ ,71.20
+ ,96.80
+ ,40.90
+ ,113.20
+ ,92.40
+ ,96.80
+ ,38.30
+ ,105.90
+ ,81.50
+ ,96.80
+ ,46.30
+ ,108.80
+ ,85.30
+ ,96.80
+ ,28.40
+ ,102.30
+ ,69.90
+ ,96.80
+ ,78.40
+ ,99.00
+ ,84.20
+ ,96.90
+ ,36.80
+ ,100.70
+ ,90.70
+ ,97.10
+ ,50.70
+ ,115.50
+ ,100.30
+ ,97.10
+ ,42.80)
+ ,dim=c(4
+ ,61)
+ ,dimnames=list(c('X1'
+ ,'x2'
+ ,'X3'
+ ,'X4')
+ ,1:61))
>  y <- array(NA,dim=c(4,61),dimnames=list(c('X1','x2','X3','X4'),1:61))
>  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: Wessa P., (2007), Notched Boxplots (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_notchedbox1.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description
> z <- as.data.frame(t(y))
> z <- log(z)
> postscript(file="/var/www/html/rcomp/tmp/1nnny1225739912.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,] 4.454347 4.197202 4.567468 3.346389
[2,] 4.566429 4.389499 4.572647 3.756538
[3,] 4.622027 4.469350 4.577799 3.998201
[4,] 4.663439 4.544358 4.605170 4.276666
[5,] 4.749271 4.697749 4.611152 4.980176

$n
[1] 61 61 61 61

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 4.602402 4.438023 4.571220 3.892980
[2,] 4.641652 4.500678 4.584378 4.103422

$out
[1] 4.394449 4.393214

$group
[1] 1 1

$names
[1] "X1" "x2" "X3" "X4"

Warning message:
In bxp(list(stats = c(4.45434729625351, 4.56642935767166, 4.62202730305451,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> 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/28hg81225739912.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/3sic01225739912.tab") 
> z
         X1       x2       X3       X4
1  4.704110 4.693181 4.604170 4.283587
2  4.568506 4.484132 4.603168 4.084294
3  4.623992 4.546481 4.603168 4.450853
4  4.665324 4.588024 4.608166 4.479607
5  4.394449 4.458988 4.604170 4.139955
6  4.550714 4.389499 4.604170 4.465908
7  4.615121 4.645352 4.605170 4.371976
8  4.695011 4.683981 4.606170 4.718499
9  4.627910 4.536891 4.606170 4.371976
10 4.507557 4.275276 4.607168 4.883559
11 4.566429 4.544358 4.608166 3.691376
12 4.565389 4.552824 4.611152 4.234107
13 4.663439 4.568506 4.605170 4.084294
14 4.635699 4.511958 4.606170 4.301359
15 4.624973 4.435567 4.607168 4.050044
16 4.651099 4.458988 4.605170 4.395683
17 4.454347 4.477337 4.606170 3.841601
18 4.522875 4.318821 4.606170 3.723281
19 4.671894 4.697749 4.606170 4.265493
20 4.723842 4.634729 4.610158 4.218036
21 4.622027 4.407938 4.610158 4.276666
22 4.521789 4.219508 4.610158 4.980176
23 4.578826 4.568506 4.567468 3.681351
24 4.574711 4.546481 4.567468 3.949319
25 4.657763 4.499810 4.572647 4.300003
26 4.631812 4.477337 4.572647 4.261270
27 4.585987 4.332048 4.573680 4.107590
28 4.649187 4.412798 4.572647 4.110874
29 4.470495 4.399375 4.572647 3.998201
30 4.498698 4.197202 4.572647 3.666122
31 4.698661 4.576771 4.572647 4.198705
32 4.715817 4.544358 4.574711 4.069027
33 4.591071 4.390739 4.574711 4.091006
34 4.573680 4.255613 4.574711 4.393214
35 4.554929 4.475062 4.572647 3.618993
36 4.574711 4.494239 4.573680 3.797734
37 4.724729 4.601162 4.576771 3.885679
38 4.633758 4.433195 4.577799 3.988984
39 4.578826 4.318821 4.577799 3.901973
40 4.713127 4.521789 4.576771 4.120662
41 4.470495 4.391977 4.577799 3.555348
42 4.572647 4.291828 4.577799 3.575151
43 4.737075 4.603168 4.577799 3.937691
44 4.703204 4.499810 4.577799 3.891820
45 4.643429 4.420045 4.577799 3.725693
46 4.621044 4.282206 4.577799 4.283587
47 4.549657 4.366913 4.585987 3.740048
48 4.563306 4.469350 4.572647 3.786460
49 4.651099 4.510860 4.572647 3.808882
50 4.632785 4.383276 4.572647 3.918005
51 4.585987 4.298645 4.572647 3.711130
52 4.735321 4.458988 4.572647 3.854394
53 4.393214 4.310799 4.572647 3.608212
54 4.561218 4.265493 4.572647 3.711130
55 4.729156 4.526127 4.572647 3.645450
56 4.662495 4.400603 4.572647 3.835142
57 4.689511 4.446174 4.572647 3.346389
58 4.627910 4.247066 4.572647 4.361824
59 4.595120 4.433195 4.573680 3.605498
60 4.612146 4.507557 4.575741 3.925926
61 4.749271 4.608166 4.575741 3.756538
> system("convert tmp/1nnny1225739912.ps tmp/1nnny1225739912.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.598   0.160   0.657