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(109.04 + ,92.31 + ,100.00 + ,130.33 + ,109.49 + ,89.74 + ,100.00 + ,135.50 + ,107.23 + ,88.46 + ,94.00 + ,124.18 + ,109.00 + ,89.74 + ,93.00 + ,129.93 + ,109.12 + ,88.46 + ,92.00 + ,113.90 + ,109.24 + ,88.46 + ,89.00 + ,101.49 + ,108.92 + ,87.18 + ,87.00 + ,95.46 + ,109.53 + ,89.74 + ,87.00 + ,90.52 + ,107.06 + ,87.18 + ,91.00 + ,93.85 + ,109.11 + ,87.18 + ,86.00 + ,91.73 + ,109.26 + ,85.90 + ,89.00 + ,109.71 + ,109.99 + ,82.05 + ,94.00 + ,137.51 + ,110.17 + ,79.49 + ,102.00 + ,184.89 + ,110.28 + ,78.21 + ,100.00 + ,208.33 + ,109.13 + ,79.49 + ,98.00 + ,237.56 + ,110.15 + ,83.33 + ,100.00 + ,233.83 + ,109.39 + ,84.62 + ,102.00 + ,220.10 + ,108.45 + ,83.33 + ,104.00 + ,199.71 + ,108.23 + ,79.49 + ,108.00 + ,191.73 + ,107.44 + ,79.49 + ,108.00 + ,179.38 + ,104.86 + ,84.62 + ,105.00 + ,169.67 + ,106.23 + ,91.03 + ,107.00 + ,162.72 + ,105.85 + ,93.59 + ,103.00 + ,165.19 + ,104.95 + ,94.87 + ,110.00 + ,152.61 + ,104.46 + ,94.87 + ,109.00 + ,143.71 + ,104.66 + ,94.87 + ,109.00 + ,134.81 + ,103.05 + ,94.87 + ,110.00 + ,140.15 + ,104.16 + ,93.59 + ,112.00 + ,132.80 + ,104.08 + ,93.59 + ,113.00 + ,128.55 + ,104.20 + ,94.87 + ,113.00 + ,128.37 + ,103.68 + ,97.44 + ,110.00 + ,120.91 + ,103.69 + ,97.44 + ,112.00 + ,115.80 + ,101.29 + ,97.44 + ,110.00 + ,109.53 + ,103.03 + ,97.44 + ,103.00 + ,119.18 + ,102.90 + ,100.00 + ,112.00 + ,114.19 + ,102.68 + ,101.28 + ,113.00 + ,112.92 + ,102.98 + ,103.85 + ,109.00 + ,117.29 + ,103.47 + ,105.13 + ,109.00 + ,130.61 + ,101.72 + ,105.13 + ,109.00 + ,131.94 + ,102.82 + ,103.85 + ,109.00 + ,125.16 + ,102.74 + ,103.85 + ,105.00 + ,127.17 + ,102.38 + ,103.85 + ,107.00 + ,124.24 + ,101.81 + ,103.85 + ,106.00 + ,113.15 + ,101.88 + ,103.85 + ,109.00 + ,111.32 + ,99.60 + ,105.13 + ,110.00 + ,113.04 + ,100.93 + ,106.41 + ,106.00 + ,103.68 + ,100.85 + ,107.69 + ,102.00 + ,100.40 + ,100.93 + ,108.97 + ,103.00 + ,104.31 + ,101.10 + ,110.26 + ,97.00 + ,109.13 + ,101.10 + ,108.97 + ,101.00 + ,109.48 + ,99.31 + ,106.41 + ,100.00 + ,102.81 + ,100.33 + ,100.00 + ,100.00 + ,100.00 + ,99.99 + ,100.00 + ,100.00 + ,92.53 + ,99.82 + ,102.56 + ,106.00 + ,96.84 + ,99.65 + ,110.26 + ,109.00 + ,97.30 + ,99.06 + ,114.10 + ,108.00 + ,87.42 + ,96.92 + ,114.10 + ,105.00 + ,85.18 + ,98.20 + ,110.26 + ,105.00 + ,79.09 + ,98.54 + ,106.41 + ,104.00 + ,83.29 + ,98.71 + ,106.41 + ,105.00 + ,88.11 + ,98.20 + ,106.41 + ,109.00 + ,81.39 + ,98.29 + ,107.69 + ,109.00 + ,81.96 + ,96.67 + ,108.97 + ,107.00 + ,78.06 + ,97.69 + ,107.69 + ,111.00 + ,75.70 + ,97.78 + ,110.26 + ,105.00 + ,77.31 + ,97.44 + ,108.97 + ,107.00 + ,73.23 + ,96.92 + ,108.97 + ,108.00 + ,72.66 + ,96.84 + ,107.69 + ,110.00 + ,68.18 + ,95.05 + ,108.97 + ,108.00 + ,67.26 + ,96.33 + ,108.97 + ,105.00 + ,63.58 + ,96.33 + ,108.97 + ,105.00 + ,61.52 + ,96.16 + ,108.97 + ,96.00 + ,60.54 + ,96.50 + ,108.97 + ,106.00 + ,56.86 + ,96.33 + ,108.97 + ,100.00 + ,59.74 + ,94.71 + ,108.97 + ,98.00 + ,58.07 + ,95.82 + ,108.97 + ,101.00 + ,57.04 + ,95.47 + ,110.26 + ,102.00 + ,55.03 + ,95.82 + ,107.69 + ,100.00 + ,54.05 + ,95.99 + ,103.85 + ,93.00 + ,60.20 + ,95.73 + ,102.56 + ,98.00 + ,63.70 + ,93.77 + ,102.56 + ,102.00 + ,60.25 + ,94.71 + ,102.56 + ,103.00 + ,55.77 + ,94.62 + ,102.56 + ,107.00 + ,52.38 + ,94.79 + ,101.28 + ,108.00 + ,55.60 + ,94.88 + ,100.00 + ,108.00 + ,56.58 + ,94.79 + ,100.00 + ,108.00 + ,54.05 + ,93.43 + ,101.28 + ,110.00 + ,53.13 + ,94.37 + ,103.85 + ,111.00 + ,51.41 + ,94.62 + ,102.56 + ,112.00 + ,52.27 + ,94.45 + ,97.44 + ,111.00 + ,52.10 + ,94.37 + ,93.59 + ,113.00 + ,49.91 + ,94.20 + ,89.74 + ,112.00 + ,44.63 + ,92.66 + ,87.18 + ,110.00 + ,43.60 + ,93.51 + ,89.74 + ,103.00 + ,42.68 + ,93.60 + ,91.03 + ,93.00 + ,43.02 + ,93.60 + ,92.31 + ,97.00 + ,45.89 + ,93.77 + ,91.03 + ,107.00 + ,51.98 + ,93.60 + ,88.46 + ,111.00 + ,52.61 + ,92.41 + ,85.90 + ,115.00 + ,51.98 + ,93.60 + ,85.90 + ,115.00 + ,54.57 + ,93.34 + ,84.62 + ,114.00 + ,55.83 + ,92.92 + ,88.46 + ,114.00 + ,53.99 + ,92.07 + ,93.59 + ,118.00 + ,53.53 + ,91.89 + ,96.15 + ,119.00 + ,56.35 + ,90.27 + ,93.59 + ,124.00 + ,54.62 + ,91.72 + ,91.03 + ,126.00 + ,54.11 + ,91.98 + ,88.46 + ,125.00 + ,62.15 + ,91.81 + ,91.03 + ,125.00 + ,61.29 + ,91.98 + ,96.15 + ,121.00 + ,62.21 + ,91.30 + ,98.72 + ,127.00 + ,58.47 + ,89.93 + ,100.00 + ,124.00 + ,56.98 + ,90.87 + ,100.00 + ,126.00 + ,58.59 + ,90.53 + ,98.72 + ,124.00 + ,56.52 + ,90.27 + ,100.00 + ,123.00 + ,51.46 + ,90.10 + ,100.00 + ,124.00 + ,56.35 + ,89.68 + ,101.28 + ,122.00 + ,56.00 + ,87.89 + ,101.28 + ,120.00 + ,53.93) + ,dim=c(4 + ,117) + ,dimnames=list(c('icp' + ,'wtb' + ,'icv' + ,'icg ') + ,1:117)) > y <- array(NA,dim=c(4,117),dimnames=list(c('icp','wtb','icv','icg '),1:117)) > 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/1x66e1257601691.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,] 87.89 78.21 89 42.68 [2,] 94.37 91.03 102 56.35 [3,] 98.29 100.00 107 83.29 [4,] 103.69 106.41 111 119.18 [5,] 110.28 114.10 124 208.33 $n [1] 117 117 117 117 $conf [,1] [,2] [,3] [,4] [1,] 96.92862 97.75343 105.6854 74.11236 [2,] 99.65138 102.24657 108.3146 92.46764 $out [1] 87.00 87.00 86.00 126.00 125.00 125.00 127.00 126.00 237.56 233.83 [11] 220.10 $group [1] 3 3 3 3 3 3 3 3 4 4 4 $names [1] "icp" "wtb" "icv" "icg\r" Warning messages: 1: In axis(side = 1, at = 1:4, labels = c("icp", "wtb", "icv", "icg\r" : font width unknown for character 0xd 2: In axis(side = 1, at = 1:4, labels = c("icp", "wtb", "icv", "icg\r" : font width unknown for character 0xd > 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/2wrle1257601691.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/3m0en1257601691.tab") > > system("convert tmp/1x66e1257601691.ps tmp/1x66e1257601691.png") > > > proc.time() user system elapsed 0.426 0.151 0.508