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.00
+ ,100.00
+ ,100.00
+ ,100.00
+ ,99.80
+ ,93.54
+ ,98.19
+ ,106.88
+ ,108.48
+ ,110.15
+ ,111.50
+ ,119.65
+ ,104.62
+ ,107.83
+ ,100.13
+ ,114.20
+ ,104.42
+ ,100.91
+ ,103.31
+ ,120.90
+ ,105.93
+ ,98.93
+ ,103.15
+ ,119.85
+ ,100.78
+ ,97.09
+ ,97.53
+ ,119.21
+ ,93.79
+ ,96.17
+ ,92.94
+ ,111.64
+ ,108.90
+ ,100.38
+ ,109.73
+ ,116.26
+ ,112.54
+ ,111.75
+ ,109.52
+ ,127.51
+ ,105.96
+ ,104.22
+ ,105.90
+ ,122.40
+ ,104.64
+ ,94.75
+ ,101.79
+ ,112.82
+ ,102.88
+ ,108.05
+ ,106.03
+ ,114.37
+ ,103.99
+ ,101.59
+ ,104.46
+ ,110.86
+ ,113.18
+ ,108.50
+ ,113.29
+ ,119.31
+ ,102.73
+ ,100.44
+ ,98.96
+ ,114.98
+ ,100.34
+ ,98.46
+ ,98.00
+ ,120.29
+ ,101.65
+ ,92.69
+ ,99.75
+ ,111.69
+ ,103.57
+ ,95.10
+ ,97.55
+ ,115.25
+ ,90.15
+ ,91.84
+ ,86.85
+ ,107.84
+ ,105.22
+ ,105.45
+ ,108.29
+ ,131.20
+ ,113.31
+ ,113.78
+ ,113.14
+ ,133.01
+ ,106.26
+ ,105.57
+ ,105.86
+ ,120.12
+ ,109.33
+ ,108.45
+ ,102.26
+ ,125.77
+ ,100.27
+ ,104.13
+ ,102.97
+ ,108.62
+ ,102.98
+ ,105.75
+ ,101.70
+ ,117.83
+ ,120.01
+ ,123.59
+ ,119.62
+ ,139.74
+ ,106.87
+ ,112.81
+ ,106.36
+ ,127.02
+ ,107.83
+ ,106.06
+ ,108.33
+ ,104.77
+ ,124.22
+ ,110.57
+ ,120.33
+ ,128.10
+ ,112.58
+ ,107.79
+ ,111.68
+ ,119.72
+ ,106.51
+ ,112.76
+ ,105.13
+ ,112.75
+ ,123.63
+ ,114.25
+ ,127.64
+ ,122.67
+ ,121.00
+ ,121.63
+ ,125.06
+ ,122.53
+ ,124.90
+ ,124.17
+ ,126.33
+ ,125.52
+ ,125.56
+ ,120.85
+ ,121.09
+ ,126.28
+ ,113.24
+ ,111.37
+ ,117.06
+ ,107.96
+ ,114.91
+ ,109.08
+ ,117.37
+ ,111.10
+ ,130.39
+ ,121.77
+ ,132.52
+ ,130.34
+ ,124.95
+ ,117.97
+ ,124.47
+ ,131.79
+ ,117.34
+ ,114.68
+ ,117.27
+ ,135.00
+ ,125.62
+ ,125.03
+ ,126.56
+ ,138.71
+ ,115.85
+ ,127.69
+ ,117.11
+ ,133.70
+ ,112.42
+ ,135.66
+ ,117.50
+ ,136.23
+ ,130.22
+ ,139.17
+ ,136.74
+ ,145.22
+ ,127.91
+ ,141.50
+ ,134.90
+ ,139.15
+ ,132.14
+ ,146.51
+ ,139.80
+ ,146.48
+ ,132.10
+ ,140.94
+ ,133.50
+ ,143.97
+ ,129.99
+ ,145.58
+ ,137.03
+ ,131.81
+ ,129.29
+ ,129.10
+ ,133.72
+ ,138.03
+ ,153.11
+ ,151.80
+ ,157.90
+ ,165.09
+ ,135.42
+ ,134.03
+ ,132.11
+ ,148.71
+ ,148.19
+ ,141.96
+ ,138.75
+ ,162.59
+ ,144.97
+ ,144.58
+ ,142.57
+ ,168.41
+ ,136.11
+ ,137.83
+ ,127.31
+ ,152.12
+ ,133.22
+ ,146.95
+ ,133.75
+ ,155.17
+ ,143.11
+ ,147.57
+ ,146.04
+ ,162.56
+ ,149.22
+ ,160.25
+ ,151.37
+ ,161.75
+ ,150.58
+ ,157.21
+ ,156.68
+ ,168.31
+ ,140.48
+ ,143.81
+ ,142.48
+ ,161.78
+ ,143.36
+ ,151.64
+ ,148.66
+ ,151.17
+ ,145.14
+ ,139.68
+ ,143.86
+ ,152.12
+ ,164.47
+ ,152.00
+ ,162.23
+ ,186.83
+ ,146.52
+ ,144.75
+ ,141.08
+ ,161.83
+ ,157.00
+ ,155.19
+ ,148.50
+ ,175.46
+ ,159.07
+ ,148.77
+ ,151.55
+ ,187.40
+ ,149.54
+ ,152.30
+ ,145.91
+ ,172.96
+ ,148.08
+ ,151.23
+ ,143.60
+ ,175.53
+ ,153.53
+ ,151.87
+ ,154.90
+ ,176.86
+ ,167.51
+ ,165.59
+ ,169.42
+ ,190.10
+ ,169.99
+ ,169.09
+ ,174.89
+ ,195.36
+ ,149.81
+ ,153.23
+ ,151.50
+ ,173.96)
+ ,dim=c(4
+ ,72)
+ ,dimnames=list(c('InvwEU'
+ ,'InvwnietEU'
+ ,'UitvwEU'
+ ,'UitvwnietEU')
+ ,1:72))
>  y <- array(NA,dim=c(4,72),dimnames=list(c('InvwEU','InvwnietEU','UitvwEU','UitvwnietEU'),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 = '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/10gm11258044188.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,]  90.150  91.840  86.850 100.000
[2,] 105.945 105.660 105.495 118.520
[3,] 120.505 119.410 119.975 130.770
[4,] 141.795 144.665 140.440 153.645
[5,] 169.990 169.090 174.890 195.360

$n
[1] 72 72 72 72

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 113.8296 112.1471 113.4681 124.2296
[2,] 127.1804 126.6729 126.4819 137.3104

$out
numeric(0)

$group
numeric(0)

$names
[1] "InvwEU"      "InvwnietEU"  "UitvwEU"     "UitvwnietEU"

> 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/29y9i1258044188.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/3ynou1258044188.tab") 
> 
> system("convert tmp/10gm11258044188.ps tmp/10gm11258044188.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.409   0.163   0.485