R version 2.8.0 (2008-10-20)
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(100.30
+ ,100.00
+ ,100.00
+ ,NA
+ ,101.90
+ ,102.83
+ ,97.07
+ ,NA
+ ,102.10
+ ,109.50
+ ,92.57
+ ,NA
+ ,103.20
+ ,115.91
+ ,87.69
+ ,NA
+ ,103.70
+ ,107.94
+ ,82.97
+ ,NA
+ ,106.20
+ ,110.86
+ ,78.80
+ ,100.00
+ ,107.70
+ ,118.89
+ ,75.16
+ ,101.42
+ ,109.90
+ ,123.38
+ ,72.07
+ ,97.84
+ ,111.70
+ ,113.33
+ ,68.06
+ ,100.00
+ ,114.90
+ ,116.38
+ ,65.13
+ ,95.64
+ ,116.00
+ ,122.04
+ ,63.94
+ ,98.94
+ ,118.30
+ ,125.47
+ ,63.29
+ ,98.39
+ ,120.40
+ ,115.62
+ ,62.91
+ ,99.98
+ ,126.00
+ ,117.91
+ ,62.96
+ ,96.14
+ ,128.10
+ ,122.40
+ ,65.24
+ ,99.33
+ ,130.10
+ ,125.05
+ ,69.96
+ ,97.82
+ ,130.80
+ ,114.18
+ ,76.14
+ ,102.85
+ ,133.60
+ ,114.74
+ ,84.06
+ ,103.93
+ ,134.20
+ ,120.63
+ ,90.56
+ ,107.88
+ ,135.50
+ ,123.68
+ ,94.69
+ ,110.82
+ ,136.20
+ ,112.84
+ ,98.54
+ ,115.53
+ ,139.10
+ ,115.64
+ ,99.57
+ ,117.09
+ ,139.00
+ ,122.32
+ ,101.52
+ ,123.40
+ ,139.60
+ ,124.59
+ ,106.07
+ ,126.57
+ ,138.70
+ ,116.33
+ ,109.00
+ ,130.34
+ ,140.90
+ ,117.45
+ ,104.83
+ ,81.79
+ ,141.30
+ ,125.64
+ ,103.96
+ ,83.26
+ ,141.80
+ ,128.38
+ ,103.04
+ ,83.65
+ ,142.00
+ ,119.87
+ ,99.19
+ ,83.88
+ ,144.50
+ ,121.22
+ ,97.89
+ ,84.73
+ ,144.60
+ ,128.98
+ ,93.82
+ ,86.27
+ ,145.50
+ ,131.35
+ ,89.48
+ ,87.71
+ ,146.80
+ ,121.35
+ ,87.74
+ ,88.52
+ ,149.50
+ ,123.72
+ ,85.79
+ ,89.32
+ ,149.90
+ ,131.06
+ ,85.47
+ ,90.61
+ ,150.10
+ ,134.55
+ ,82.00
+ ,90.95
+ ,150.90
+ ,125.93
+ ,83.51
+ ,90.90
+ ,152.80
+ ,128.90
+ ,82.38
+ ,91.02
+ ,153.10
+ ,136.19
+ ,79.88
+ ,91.66
+ ,154.00
+ ,140.34
+ ,76.03
+ ,91.23
+ ,154.90
+ ,130.48
+ ,70.44
+ ,90.81
+ ,156.90
+ ,134.68
+ ,68.71
+ ,90.90
+ ,158.40
+ ,141.05
+ ,65.08
+ ,91.27
+ ,159.70
+ ,145.44
+ ,62.09
+ ,90.45
+ ,160.20
+ ,136.21
+ ,59.82
+ ,89.60
+ ,163.20
+ ,139.85
+ ,58.35
+ ,88.49
+ ,163.70
+ ,147.13
+ ,59.38
+ ,88.31
+ ,164.40
+ ,151.44
+ ,58.57
+ ,87.83
+ ,163.70
+ ,143.62
+ ,58.68
+ ,87.62
+ ,165.50
+ ,148.55
+ ,57.48
+ ,87.60
+ ,165.60
+ ,153.54
+ ,55.69
+ ,88.08
+ ,166.80
+ ,159.79
+ ,54.39
+ ,88.08
+ ,167.50
+ ,152.55
+ ,53.52
+ ,87.41
+ ,170.60
+ ,155.84
+ ,52.01
+ ,87.53
+ ,170.90
+ ,160.38
+ ,49.51
+ ,88.31
+ ,172.00
+ ,164.22
+ ,49.46
+ ,89.04
+ ,171.80
+ ,156.40
+ ,48.21
+ ,89.34
+ ,173.90
+ ,160.05
+ ,47.83
+ ,90.22
+ ,174.00
+ ,165.60
+ ,48.86
+ ,91.02
+ ,173.80
+ ,171.15
+ ,49.51
+ ,91.69
+ ,173.90
+ ,161.90
+ ,49.95
+ ,91.62
+ ,176.00
+ ,167.21
+ ,49.57
+ ,91.32
+ ,176.60
+ ,171.34
+ ,51.52
+ ,91.82
+ ,178.20
+ ,176.83
+ ,48.48
+ ,91.75
+ ,179.20
+ ,166.27
+ ,50.43
+ ,91.57
+ ,181.30
+ ,172.30
+ ,48.26
+ ,92.19
+ ,181.80
+ ,176.71
+ ,48.70
+ ,92.81
+ ,182.90
+ ,182.99
+ ,47.72
+ ,50.99
+ ,183.80
+ ,172.07
+ ,45.55
+ ,50.92
+ ,186.30
+ ,178.17
+ ,45.82
+ ,50.37
+ ,187.40
+ ,182.20
+ ,44.41
+ ,50.32
+ ,189.20
+ ,188.49
+ ,45.44
+ ,49.77
+ ,189.70
+ ,176.88
+ ,45.17
+ ,49.45
+ ,191.90
+ ,182.13
+ ,45.39
+ ,49.15
+ ,192.60
+ ,185.32
+ ,46.10
+ ,49.10
+ ,193.70
+ ,192.86
+ ,50.27
+ ,48.85
+ ,194.20
+ ,180.27
+ ,50.54
+ ,48.76
+ ,197.60
+ ,184.92
+ ,53.04
+ ,48.58
+ ,199.30
+ ,187.82
+ ,53.04
+ ,48.90
+ ,201.40
+ ,194.94
+ ,52.82
+ ,48.97
+ ,203.00
+ ,184.36
+ ,52.98
+ ,48.78
+ ,206.30
+ ,188.80
+ ,51.52
+ ,48.64
+ ,207.10
+ ,193.42
+ ,50.70
+ ,48.64
+ ,209.80
+ ,199.76
+ ,49.51
+ ,48.16
+ ,211.10
+ ,188.78
+ ,51.08
+ ,48.12
+ ,215.30
+ ,191.49
+ ,53.69
+ ,48.05
+ ,217.40
+ ,194.87
+ ,58.30
+ ,48.25
+ ,215.50
+ ,198.28
+ ,64.05
+ ,47.77
+ ,210.90
+ ,183.24
+ ,72.45
+ ,46.72
+ ,212.60
+ ,204.87
+ ,80.91
+ ,NA)
+ ,dim=c(4
+ ,90)
+ ,dimnames=list(c('RP'
+ ,'FCE'
+ ,'UNEM'
+ ,'IS')
+ ,1:90))
>  y <- array(NA,dim=c(4,90),dimnames=list(c('RP','FCE','UNEM','IS'),1:90))
>  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/freestat/rcomp/tmp/1brtl1257971368.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,] 100.3 100.000  44.410  46.720
[2,] 139.0 122.320  50.540  50.955
[3,] 161.7 142.335  62.935  89.330
[4,] 182.9 176.830  83.510  92.500
[5,] 217.4 204.870 109.000 130.340

$n
[1] 90 90 90 84

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 154.3886 133.2565 57.44396 82.16797
[2,] 169.0114 151.4135 68.42604 96.49203

$out
numeric(0)

$group
numeric(0)

$names
[1] "RP"   "FCE"  "UNEM" "IS"  

Warning message:
In bxp(list(stats = c(100.3, 139, 161.7, 182.9, 217.4, 100, 122.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/2dg651257971368.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/freestat/rcomp/tmp/34kxy1257971368.tab") 
> 
> system("convert tmp/1brtl1257971368.ps tmp/1brtl1257971368.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.592   0.242   0.673