R version 2.11.1 (2010-05-31)
Copyright (C) 2010 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(68
+ ,63
+ ,4
+ ,48
+ ,62
+ ,8
+ ,51
+ ,61
+ ,4
+ ,46
+ ,70
+ ,4
+ ,56
+ ,60
+ ,6
+ ,46
+ ,76
+ ,4
+ ,44
+ ,68
+ ,8
+ ,55
+ ,71
+ ,4
+ ,67
+ ,77
+ ,4
+ ,49
+ ,65
+ ,8
+ ,54
+ ,69
+ ,8
+ ,75
+ ,78
+ ,4
+ ,61
+ ,65
+ ,5
+ ,51
+ ,75
+ ,4
+ ,52
+ ,64
+ ,4
+ ,38
+ ,52
+ ,4
+ ,46
+ ,63
+ ,4
+ ,52
+ ,69
+ ,4
+ ,52
+ ,63
+ ,4
+ ,48
+ ,53
+ ,4
+ ,76
+ ,79
+ ,4
+ ,64
+ ,72
+ ,4
+ ,30
+ ,74
+ ,4
+ ,42
+ ,55
+ ,10
+ ,50
+ ,73
+ ,8
+ ,39
+ ,66
+ ,4
+ ,55
+ ,76
+ ,7
+ ,59
+ ,64
+ ,8
+ ,18
+ ,55
+ ,4
+ ,49
+ ,64
+ ,4
+ ,42
+ ,76
+ ,5
+ ,80
+ ,83
+ ,4
+ ,66
+ ,61
+ ,4
+ ,69
+ ,67
+ ,5
+ ,66
+ ,61
+ ,4
+ ,63
+ ,61
+ ,4
+ ,33
+ ,55
+ ,4
+ ,38
+ ,53
+ ,5
+ ,57
+ ,68
+ ,4
+ ,46
+ ,71
+ ,12
+ ,58
+ ,65
+ ,4
+ ,37
+ ,72
+ ,9
+ ,59
+ ,54
+ ,15
+ ,65
+ ,71
+ ,4
+ ,37
+ ,76
+ ,4
+ ,39
+ ,70
+ ,4
+ ,62
+ ,71
+ ,4
+ ,64
+ ,62
+ ,5
+ ,52
+ ,74
+ ,7
+ ,59
+ ,65
+ ,4
+ ,53
+ ,70
+ ,4
+ ,55
+ ,75
+ ,4
+ ,58
+ ,70
+ ,6
+ ,72
+ ,72
+ ,6
+ ,36
+ ,62
+ ,16
+ ,58
+ ,64
+ ,4
+ ,56
+ ,64
+ ,6
+ ,63
+ ,71
+ ,4
+ ,51
+ ,58
+ ,9
+ ,50
+ ,58
+ ,5
+ ,44
+ ,59
+ ,4
+ ,59
+ ,75
+ ,5
+ ,58
+ ,79
+ ,5
+ ,73
+ ,81
+ ,4
+ ,37
+ ,49
+ ,4
+ ,38
+ ,54
+ ,4
+ ,48
+ ,64
+ ,5
+ ,62
+ ,67
+ ,4
+ ,53
+ ,65
+ ,4
+ ,54
+ ,66
+ ,8
+ ,51
+ ,63
+ ,4
+ ,73
+ ,75
+ ,4
+ ,51
+ ,62
+ ,4
+ ,70
+ ,61
+ ,4
+ ,47
+ ,65
+ ,6
+ ,57
+ ,64
+ ,4
+ ,64
+ ,64
+ ,4
+ ,47
+ ,65
+ ,4
+ ,54
+ ,55
+ ,18
+ ,74
+ ,72
+ ,7
+ ,59
+ ,73
+ ,4
+ ,47
+ ,54
+ ,4
+ ,36
+ ,67
+ ,5
+ ,55
+ ,64
+ ,4
+ ,62
+ ,75
+ ,4
+ ,44
+ ,59
+ ,8
+ ,67
+ ,67
+ ,10
+ ,42
+ ,54
+ ,5
+ ,70
+ ,77
+ ,5
+ ,47
+ ,60
+ ,5
+ ,46
+ ,58
+ ,8
+ ,60
+ ,69
+ ,5
+ ,46
+ ,55
+ ,8
+ ,50
+ ,60
+ ,4
+ ,41
+ ,56
+ ,5
+ ,63
+ ,66
+ ,4
+ ,52
+ ,69
+ ,4
+ ,68
+ ,69
+ ,7
+ ,60
+ ,75
+ ,5
+ ,32
+ ,57
+ ,10
+ ,40
+ ,61
+ ,14
+ ,47
+ ,65
+ ,4
+ ,41
+ ,59
+ ,8
+ ,67
+ ,74
+ ,5
+ ,54
+ ,68
+ ,8
+ ,46
+ ,69
+ ,7
+ ,42
+ ,43
+ ,4
+ ,67
+ ,76
+ ,4
+ ,51
+ ,70
+ ,6
+ ,57
+ ,57
+ ,4
+ ,60
+ ,67
+ ,4
+ ,51
+ ,71
+ ,8
+ ,49
+ ,73
+ ,7
+ ,62
+ ,62
+ ,6
+ ,52
+ ,72
+ ,7
+ ,60
+ ,64
+ ,4
+ ,50
+ ,59
+ ,6
+ ,59
+ ,58
+ ,9
+ ,47
+ ,70
+ ,4
+ ,56
+ ,78
+ ,4
+ ,59
+ ,66
+ ,8
+ ,62
+ ,72
+ ,4
+ ,64
+ ,73
+ ,4
+ ,55
+ ,62
+ ,5
+ ,52
+ ,61
+ ,10
+ ,48
+ ,50
+ ,4
+ ,60
+ ,63
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,51
+ ,66
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,63
+ ,68
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,49
+ ,81
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,52
+ ,72
+ ,5
+ ,NA
+ ,NA
+ ,NA
+ ,48
+ ,53
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,50
+ ,61
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,67
+ ,77
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,44
+ ,75
+ ,7
+ ,NA
+ ,NA
+ ,NA
+ ,51
+ ,70
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,37
+ ,63
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,51
+ ,57
+ ,10
+ ,NA
+ ,NA
+ ,NA
+ ,60
+ ,70
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,38
+ ,67
+ ,5
+ ,NA
+ ,NA
+ ,NA
+ ,52
+ ,44
+ ,12
+ ,NA
+ ,NA
+ ,NA
+ ,65
+ ,81
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,70
+ ,71
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,44
+ ,67
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,50
+ ,61
+ ,7
+ ,NA
+ ,NA
+ ,NA
+ ,50
+ ,56
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,57
+ ,74
+ ,11
+ ,NA
+ ,NA
+ ,NA
+ ,63
+ ,68
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,36
+ ,60
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,54
+ ,72
+ ,7
+ ,NA
+ ,NA
+ ,NA
+ ,36
+ ,74
+ ,5
+ ,NA
+ ,NA
+ ,NA
+ ,70
+ ,73
+ ,8
+ ,NA
+ ,NA
+ ,NA
+ ,47
+ ,58
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,52
+ ,67
+ ,5
+ ,NA
+ ,NA
+ ,NA
+ ,52
+ ,76
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,69
+ ,67
+ ,4
+ ,NA
+ ,NA
+ ,NA
+ ,52
+ ,68
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,48
+ ,71
+ ,16
+ ,NA
+ ,NA
+ ,NA
+ ,51
+ ,70
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,53
+ ,61
+ ,6
+ ,NA
+ ,NA
+ ,NA
+ ,55
+ ,54
+ ,4
+ ,NA
+ ,NA
+ ,NA)
+ ,dim=c(6
+ ,98)
+ ,dimnames=list(c('I-M'
+ ,'E-M'
+ ,'A-M'
+ ,'I-V'
+ ,'E-V'
+ ,'A-V')
+ ,1:98))
>  y <- array(NA,dim=c(6,98),dimnames=list(c('I-M','E-M','A-M','I-V','E-V','A-V'),1:98))
>  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/rcomp/tmp/1vodp1287221740.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,height=8.3333333333333) 
> (r<-boxplot(z ,xlab=xlab,ylab=ylab,main=main,notch=TRUE,col=par1))
$stats
     [,1] [,2] [,3] [,4] [,5] [,6]
[1,]   30   49    4   32 50.0  4.0
[2,]   47   61    4   47 61.0  4.0
[3,]   52   67    5   55 66.0  4.0
[4,]   59   72    8   63 71.5  5.5
[5,]   76   81   14   80 83.0  7.0

$n
[1] 98 98 98 63 63 63

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 50.08475 65.24435 4.361584 51.81502 63.90986 3.701408
[2,] 53.91525 68.75565 5.638416 58.18498 68.09014 4.298592

$out
 [1] 18 43 44 15 16 18 16  8  8 10  8 12  9  8  8 10  8  9

$group
 [1] 1 2 2 3 3 3 3 6 6 6 6 6 6 6 6 6 6 6

$names
[1] "I-M" "E-M" "A-M" "I-V" "E-V" "A-V"

Warning message:
In bxp(list(stats = c(30, 47, 52, 59, 76, 49, 61, 67, 72, 81, 4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/2h6cv1287221740.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/rcomp/tmp/3k7aj1287221740.tab") 
> try(system("convert tmp/1vodp1287221740.ps tmp/1vodp1287221740.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.750   0.260   0.991