R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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 <- c(20 + ,25 + ,15 + ,15 + ,25 + ,25 + ,25 + ,21 + ,30 + ,25 + ,20 + ,40 + ,13 + ,30 + ,25 + ,20 + ,25 + ,20 + ,25 + ,20 + ,20 + ,15 + ,15 + ,12 + ,20 + ,5 + ,20 + ,15 + ,25 + ,22 + ,20 + ,22 + ,25 + ,20 + ,20 + ,35 + ,30 + ,25 + ,20 + ,20 + ,20 + ,25 + ,25 + ,15 + ,20 + ,35 + ,25 + ,25 + ,30 + ,23 + ,10 + ,22 + ,25 + ,25 + ,22 + ,30 + ,20 + ,25 + ,25 + ,22 + ,25 + ,25 + ,25 + ,22 + ,25 + ,12 + ,18 + ,20 + ,20 + ,22 + ,30 + ,25 + ,22 + ,20 + ,50 + ,30 + ,25 + ,20 + ,30 + ,22 + ,25 + ,30 + ,22 + ,25 + ,22 + ,22 + ,25 + ,25 + ,25 + ,20 + ,22 + ,15 + ,20 + ,30 + ,20 + ,25 + ,30 + ,35 + ,22 + ,12 + ,30 + ,15 + ,10 + ,30 + ,9 + ,25 + ,20 + ,20 + ,35 + ,25 + ,35 + ,30 + ,12 + ,25 + ,15 + ,25 + ,25 + ,20 + ,20 + ,6 + ,15 + ,40 + ,20 + ,40 + ,25 + ,25 + ,20 + ,15 + ,15 + ,22 + ,24 + ,22 + ,20 + ,25 + ,25 + ,25 + ,35 + ,40 + ,20 + ,22 + ,22 + ,20 + ,25 + ,25 + ,18 + ,25 + ,20 + ,25 + ,30 + ,20 + ,22 + ,35 + ,22 + ,25 + ,25 + ,25 + ,25 + ,22 + ,23 + ,35 + ,15 + ,25 + ,18 + ,22 + ,25 + ,25 + ,28 + ,30 + ,20 + ,25 + ,25 + ,30 + ,22 + ,30 + ,10 + ,10 + ,25 + ,20 + ,22 + ,25 + ,25 + ,15 + ,22 + ,25 + ,25 + ,28 + ,22 + ,30 + ,25 + ,20 + ,25 + ,25 + ,20 + ,30 + ,20 + ,30 + ,50 + ,19 + ,20 + ,28 + ,20 + ,25 + ,35 + ,25 + ,25 + ,15 + ,16 + ,20 + ,20 + ,25 + ,30 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,25 + ,30 + ,22 + ,20 + ,25 + ,25 + ,18 + ,18 + ,20 + ,25 + ,25 + ,30 + ,25 + ,20 + ,25 + ,20 + ,20 + ,20 + ,22 + ,18 + ,22 + ,20 + ,15 + ,25 + ,25 + ,20 + ,25 + ,15 + ,22 + ,25 + ,25 + ,15 + ,12 + ,25 + ,30 + ,22 + ,15 + ,22 + ,25 + ,12 + ,18 + ,30 + ,25 + ,25 + ,40 + ,24 + ,25 + ,15 + ,25 + ,20 + ,25 + ,25 + ,25 + ,20 + ,30 + ,20 + ,25 + ,30 + ,22 + ,25 + ,25 + ,25 + ,50 + ,19 + ,50 + ,25 + ,35 + ,20 + ,20 + ,20 + ,20 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,20 + ,20 + ,25 + ,18 + ,25 + ,22 + ,22 + ,30 + ,30 + ,8 + ,20 + ,25 + ,30 + ,50 + ,22 + ,20 + ,10 + ,25 + ,25 + ,25 + ,25 + ,18 + ,25 + ,20 + ,25 + ,30 + ,18 + ,20 + ,25 + ,22 + ,22 + ,20 + ,20 + ,25 + ,20 + ,20 + ,20 + ,20 + ,25 + ,20 + ,10 + ,20 + ,25 + ,30 + ,25 + ,50 + ,30 + ,30 + ,50 + ,15 + ,25 + ,25 + ,22 + ,20 + ,22 + ,30 + ,25 + ,18 + ,22 + ,22 + ,30 + ,40 + ,25 + ,20 + ,10 + ,20 + ,9 + ,15 + ,20 + ,15 + ,20 + ,30 + ,12 + ,15 + ,12 + ,20 + ,15 + ,12 + ,25 + ,20 + ,25 + ,25 + ,25 + ,30 + ,20 + ,25 + ,15 + ,15 + ,22 + ,10 + ,15 + ,10 + ,20 + ,25 + ,20 + ,20 + ,38 + ,20 + ,20 + ,20 + ,40 + ,25 + ,25 + ,30 + ,25 + ,10 + ,20 + ,25 + ,12 + ,15 + ,25 + ,20 + ,22 + ,22 + ,20 + ,25 + ,25 + ,25 + ,15 + ,40 + ,20 + ,20 + ,16 + ,25 + ,15 + ,20 + ,25 + ,20 + ,30 + ,50 + ,20 + ,25 + ,20 + ,30 + ,30 + ,25 + ,25 + ,12 + ,25 + ,25 + ,25 + ,20 + ,20 + ,20 + ,15 + ,20 + ,25 + ,15 + ,25 + ,50 + ,30 + ,20 + ,20 + ,25 + ,12 + ,15 + ,20 + ,20 + ,35 + ,22 + ,15 + ,18 + ,30 + ,22 + ,12 + ,12 + ,20 + ,20 + ,15 + ,25 + ,15 + ,20 + ,20 + ,25 + ,18 + ,30 + ,20 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,20 + ,22 + ,15 + ,15 + ,22 + ,20 + ,10 + ,25 + ,20 + ,20 + ,15 + ,12 + ,20 + ,5 + ,20 + ,15 + ,15 + ,25 + ,25 + ,25 + ,15 + ,25 + ,22 + ,25 + ,20 + ,18 + ,22 + ,25 + ,35 + ,25 + ,25 + ,25 + ,35 + ,30 + ,22 + ,30 + ,50 + ,15 + ,25 + ,24 + ,20 + ,25 + ,25 + ,25 + ,12 + ,15 + ,22 + ,25 + ,25 + ,25 + ,25 + ,15 + ,20 + ,20 + ,15 + ,35 + ,30 + ,20 + ,22 + ,65 + ,20 + ,25 + ,22 + ,20 + ,25 + ,25 + ,20 + ,25 + ,15 + ,20 + ,12 + ,15 + ,10 + ,25 + ,15 + ,30 + ,35 + ,25 + ,25 + ,25 + ,25 + ,25 + ,40 + ,40 + ,25 + ,25 + ,20 + ,25 + ,25 + ,22 + ,25 + ,30 + ,25 + ,25 + ,30 + ,25 + ,25 + ,30 + ,25 + ,25 + ,20 + ,22 + ,22 + ,20 + ,25 + ,22 + ,25 + ,22 + ,40 + ,25 + ,25 + ,25 + ,22 + ,20 + ,35 + ,20 + ,35 + ,25 + ,22 + ,25 + ,25 + ,25 + ,25 + ,25 + ,40 + ,25 + ,30 + ,25 + ,20 + ,25 + ,25 + ,30 + ,22 + ,22 + ,20 + ,15 + ,15 + ,25 + ,25 + ,20 + ,20 + ,15 + ,25 + ,15 + ,20 + ,22 + ,25 + ,15 + ,15 + ,18 + ,5 + ,15 + ,25 + ,18 + ,40 + ,25 + ,25 + ,20 + ,30 + ,20 + ,25 + ,25 + ,25 + ,22 + ,22 + ,25 + ,25 + ,30 + ,25 + ,25 + ,25 + ,25 + ,20 + ,20 + ,25 + ,25 + ,25 + ,25 + ,20 + ,30 + ,25 + ,22 + ,30 + ,20 + ,20 + ,30 + ,25 + ,25 + ,30 + ,20 + ,25 + ,25 + ,24 + ,25 + ,30 + ,18 + ,15 + ,22 + ,22 + ,25 + ,22 + ,22 + ,25 + ,15 + ,20 + ,22 + ,18 + ,35 + ,20 + ,20 + ,20 + ,25 + ,25 + ,30 + ,15 + ,25 + ,22 + ,26 + ,25 + ,20 + ,25 + ,25 + ,25 + ,22 + ,25 + ,25 + ,20 + ,22 + ,30 + ,15 + ,30 + ,25 + ,20 + ,25 + ,25 + ,35 + ,22 + ,20 + ,25 + ,20 + ,20 + ,18 + ,20 + ,22 + ,25 + ,10 + ,20 + ,25 + ,20 + ,20 + ,30 + ,25 + ,20 + ,15 + ,20 + ,25 + ,10 + ,20 + ,25 + ,22 + ,22 + ,25 + ,25 + ,15 + ,25 + ,20 + ,10 + ,25 + ,16 + ,25 + ,35 + ,25 + ,15 + ,25 + ,25 + ,30 + ,25 + ,10 + ,22 + ,20 + ,25 + ,20 + ,20 + ,25 + ,22 + ,18 + ,30 + ,19 + ,25 + ,20 + ,25 + ,20 + ,25 + ,20 + ,22 + ,12 + ,30 + ,12 + ,22 + ,25 + ,25 + ,25 + ,25 + ,30 + ,30 + ,10 + ,22 + ,22 + ,25 + ,20 + ,22 + ,20 + ,25 + ,20 + ,15 + ,25 + ,20 + ,25 + ,20 + ,30 + ,15 + ,40 + ,25 + ,20 + ,22 + ,22 + ,30 + ,20 + ,40 + ,20 + ,25 + ,20 + ,25 + ,20 + ,50 + ,50 + ,25 + ,25 + ,40 + ,30 + ,22 + ,30 + ,20 + ,25 + ,25 + ,30 + ,25 + ,25 + ,20 + ,18 + ,18 + ,28 + ,25 + ,22 + ,15 + ,40 + ,40 + ,12 + ,12 + ,18 + ,12 + ,25 + ,26 + ,18 + ,25 + ,22 + ,15 + ,25 + ,15 + ,15 + ,15 + ,25 + ,15 + ,12 + ,22 + ,20 + ,20 + ,25 + ,20 + ,12 + ,9 + ,15 + ,12 + ,15 + ,25 + ,20 + ,20 + ,15 + ,15 + ,30 + ,21 + ,25 + ,22 + ,22 + ,50 + ,15 + ,25 + ,15 + ,25 + ,22 + ,18 + ,50 + ,20 + ,50 + ,20 + ,20 + ,30 + ,25 + ,20 + ,22 + ,25 + ,50 + ,40 + ,25 + ,25 + ,25 + ,25 + ,30 + ,40 + ,25 + ,30 + ,20) > par1 = '12' > #'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!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 900 > (np <- floor(n / par1)) [1] 75 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 75 75 75 75 75 75 75 75 75 75 75 75 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 20 13 20 30 30 25 22 22 30 35 15 20 18 [2,] 25 30 5 25 23 25 20 22 35 25 40 25 25 [3,] 15 25 20 20 10 25 50 25 22 35 20 25 20 [4,] 15 20 15 20 22 22 30 25 12 30 40 25 25 [5,] 25 25 25 20 25 25 25 25 30 12 25 35 30 [6,] 25 20 22 25 25 12 20 20 15 25 25 40 20 [7,] 25 25 20 25 22 18 30 22 10 15 20 20 22 [8,] 21 20 22 15 30 20 22 15 30 25 15 22 35 [9,] 30 20 25 20 20 20 25 20 9 25 15 22 22 [10,] 25 15 20 35 25 22 30 30 25 20 22 20 25 [11,] 20 15 20 25 25 30 22 20 20 20 24 25 25 [12,] 40 12 35 25 22 25 25 25 20 6 22 25 25 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 25 20 25 20 25 20 25 15 30 25 22 20 [2,] 22 25 15 30 15 25 30 25 22 15 25 20 [3,] 23 25 22 20 16 25 25 25 15 25 25 25 [4,] 35 30 25 30 20 30 20 20 22 20 25 25 [5,] 15 22 25 50 20 22 25 25 25 25 50 25 [6,] 25 30 28 19 25 20 20 15 12 25 19 20 [7,] 18 10 22 20 30 25 20 22 18 25 50 20 [8,] 22 10 30 28 20 25 20 25 30 20 25 20 [9,] 25 25 25 20 25 18 22 25 25 30 35 20 [10,] 25 20 20 25 25 18 18 15 25 20 20 25 [11,] 28 22 25 35 25 20 22 12 40 25 20 18 [12,] 30 25 25 25 20 25 20 25 24 30 20 25 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 22 25 22 20 20 20 12 10 25 22 20 12 [2,] 22 25 22 25 22 9 25 15 25 20 25 25 [3,] 30 25 20 30 30 15 20 10 30 25 20 25 [4,] 30 25 20 25 25 20 25 20 25 25 30 25 [5,] 8 18 25 50 18 15 25 25 10 25 50 20 [6,] 20 25 20 30 22 20 25 20 20 15 20 20 [7,] 25 20 20 30 22 30 30 20 25 40 25 20 [8,] 30 25 20 50 30 12 20 38 12 20 20 15 [9,] 50 30 20 15 40 15 25 20 15 20 30 20 [10,] 22 18 25 25 25 12 15 20 25 16 30 25 [11,] 20 20 20 25 20 20 15 20 20 25 25 15 [12,] 10 25 10 22 10 15 22 40 22 15 25 25 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 50 18 25 15 20 22 25 25 22 15 20 25 [2,] 30 30 18 15 15 25 24 15 20 30 25 25 [3,] 20 22 30 22 15 35 20 20 25 35 25 20 [4,] 20 12 20 20 25 25 25 20 25 25 22 22 [5,] 25 12 25 10 25 25 25 15 20 25 25 22 [6,] 12 20 25 25 25 25 25 35 25 25 30 20 [7,] 15 20 25 20 15 35 12 30 15 25 25 25 [8,] 20 15 20 20 25 30 15 20 20 25 25 22 [9,] 20 25 20 15 22 22 22 22 12 40 30 25 [10,] 35 15 25 12 25 30 25 65 15 40 25 22 [11,] 22 20 20 20 20 50 25 20 10 25 25 40 [12,] 15 20 22 5 18 15 25 25 25 25 30 25 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 25 25 20 25 30 25 30 30 18 25 30 20 [2,] 25 25 15 15 20 25 20 18 35 20 25 22 [3,] 22 40 15 15 25 20 20 15 20 25 20 25 [4,] 20 25 25 18 25 20 30 22 20 25 25 10 [5,] 35 30 25 5 25 25 25 22 20 25 25 20 [6,] 20 25 20 15 22 25 25 25 25 22 35 25 [7,] 35 20 20 25 22 25 30 22 25 25 22 20 [8,] 25 25 15 18 25 25 20 22 30 25 20 20 [9,] 22 25 25 40 25 20 25 25 15 20 25 30 [10,] 25 30 15 25 30 30 25 15 25 22 20 25 [11,] 25 22 20 25 25 25 24 20 22 30 20 20 [12,] 25 22 22 20 25 22 25 22 26 15 18 15 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 20 10 22 25 25 15 30 40 18 26 22 20 [2,] 25 25 20 20 30 25 20 30 28 18 20 15 [3,] 10 16 25 25 30 20 40 22 25 25 20 15 [4,] 20 25 20 20 10 25 20 30 22 22 25 30 [5,] 25 35 20 22 22 20 25 20 15 15 20 21 [6,] 22 25 25 12 22 30 20 25 40 25 12 25 [7,] 22 15 22 30 25 15 25 25 40 15 9 22 [8,] 25 25 18 12 20 40 20 30 12 15 15 22 [9,] 25 25 30 22 22 25 50 25 12 15 12 50 [10,] 15 30 19 25 20 20 50 25 18 25 15 15 [11,] 25 25 25 25 25 22 25 20 12 15 25 25 [12,] 20 10 20 25 20 22 25 18 25 12 20 15 [,74] [,75] [,76] [1,] 25 25 NA [2,] 22 50 NA [3,] 18 40 NA [4,] 50 25 NA [5,] 20 25 NA [6,] 50 25 NA [7,] 20 25 NA [8,] 20 30 NA [9,] 30 40 NA [10,] 25 25 NA [11,] 20 30 NA [12,] 22 20 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/wessaorg/rcomp/tmp/1p3bd1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2ydom1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/3c8fw1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/4ybls1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 13 15 15 15 15 15 15 15 15 15 15 15 [2,] 20 20 20 20 20 20 20 20 20 20 20 20 [3,] 22 25 22 25 25 25 22 21 25 25 22 22 [4,] 25 25 25 25 25 25 25 25 25 25 25 25 [5,] 30 30 30 30 30 30 30 30 30 30 30 30 $n [1] 75 75 75 75 75 75 75 75 75 75 75 75 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 21.08779 24.08779 21.08779 24.08779 24.08779 24.08779 21.08779 20.08779 [2,] 22.91221 25.91221 22.91221 25.91221 25.91221 25.91221 22.91221 21.91221 [,9] [,10] [,11] [,12] [1,] 24.08779 24.08779 21.08779 21.08779 [2,] 25.91221 25.91221 22.91221 22.91221 $out [1] 35 12 10 12 50 10 40 5 35 40 9 35 50 10 50 35 10 35 35 40 10 40 40 12 40 [26] 35 12 10 10 50 12 35 50 50 8 50 10 50 12 10 35 5 35 12 40 12 12 35 35 12 [51] 40 12 50 10 10 50 40 35 12 35 40 9 35 10 50 12 38 12 12 40 12 9 35 50 40 [76] 12 40 40 50 12 12 50 40 35 12 35 12 65 40 50 35 12 40 50 10 40 12 40 12 35 [101] 6 10 10 10 40 5 10 12 $group [1] 1 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 4 4 [26] 4 4 4 4 4 5 5 5 5 5 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 [51] 6 6 6 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 8 9 9 9 9 [76] 9 9 9 9 9 9 9 9 10 10 10 10 10 10 10 11 11 11 11 11 11 11 12 12 12 [101] 12 12 12 12 12 12 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(13, 20, 22, 25, 30, 15, 20, 25, 25, 30, 15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5h0ja1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 15 12 15.0 15 20 18.0 20 15 9.0 6.0 15.0 20 18 18.0 [2,] 20 15 20.0 20 22 20.0 22 20 13.5 17.5 17.5 21 21 22.0 [3,] 25 20 20.0 25 24 23.5 25 22 21.0 25.0 22.0 25 25 25.0 [4,] 25 25 23.5 25 25 25.0 30 25 30.0 27.5 25.0 25 25 26.5 [5,] 30 30 25.0 30 25 30.0 30 30 35.0 35.0 25.0 25 30 30.0 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 20.0 20 19 15.0 18.0 18 12.0 12.0 15 19 18 8 [2,] 20.0 22 20 20.0 20.0 20 15.0 20.0 20 20 20 20 [3,] 23.5 25 25 22.5 23.5 21 23.5 24.5 25 25 20 22 [4,] 25.0 25 30 25.0 25.0 25 25.0 27.5 25 30 25 30 [5,] 30.0 28 35 30.0 30.0 30 25.0 30.0 30 35 25 30 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 18 20 15.0 10.0 9.0 12.0 10.0 10.0 15 20 12.0 12.0 [2,] 20 20 23.5 20.0 13.5 17.5 17.5 17.5 18 20 17.5 17.5 [3,] 25 20 25.0 22.0 15.0 23.5 20.0 23.5 21 25 20.0 20.0 [4,] 25 22 30.0 27.5 20.0 25.0 22.5 25.0 25 30 25.0 27.5 [5,] 30 25 30.0 30.0 20.0 30.0 25.0 30.0 25 30 25.0 35.0 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 12 18.0 5.0 15.0 15.0 15 15.0 10 15.0 22.0 20.0 20 [2,] 15 20.0 13.5 16.5 23.5 21 20.0 15 25.0 25.0 22.0 22 [3,] 20 23.5 17.5 21.0 25.0 25 21.0 20 25.0 25.0 23.5 25 [4,] 21 25.0 20.0 25.0 32.5 25 27.5 25 32.5 27.5 25.0 25 [5,] 30 30.0 25.0 25.0 35.0 25 35.0 25 40.0 30.0 25.0 25 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 20.0 15.0 5 22.0 20 20.0 15.0 15.0 15 18.0 15 15 [2,] 23.5 15.0 15 23.5 21 22.0 19.0 20.0 21 20.0 20 20 [3,] 25.0 20.0 19 25.0 25 25.0 22.0 23.5 25 23.5 20 22 [4,] 27.5 23.5 25 25.0 25 27.5 23.5 25.5 25 25.0 25 25 [5,] 30.0 25.0 40 25.0 30 30.0 30.0 30.0 30 30.0 30 25 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 10.0 18 20.0 20 15 20 18 12.0 12.0 9.0 15.0 18.0 [2,] 15.5 20 20.0 20 20 20 21 13.5 15.0 13.5 15.0 20.0 [3,] 25.0 21 23.5 22 22 25 25 20.0 16.5 20.0 21.5 22.0 [4,] 25.0 25 25.0 25 25 35 30 26.5 25.0 21.0 25.0 27.5 [5,] 35.0 30 30.0 30 30 50 40 40.0 26.0 25.0 30.0 30.0 [,75] [,76] [1,] 20 NA [2,] 25 NA [3,] 25 NA [4,] 35 NA [5,] 50 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [51] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [76] 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 22.71947 15.43893 18.40363 22.71947 22.63168 21.21947 21.35115 19.71947 [2,] 27.28053 24.56107 21.59637 27.28053 25.36832 25.78053 28.64885 24.28053 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 13.47424 20.43893 18.5792 23.17557 23.17557 22.94752 21.21947 23.63168 [2,] 28.52576 29.56107 25.4208 26.82443 26.82443 27.05248 25.78053 26.36832 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 20.43893 20.21947 21.21947 18.71947 18.93893 21.0792 22.71947 20.43893 [2,] 29.56107 24.78053 25.78053 23.28053 28.06107 27.9208 27.28053 29.56107 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 17.71947 17.43893 22.71947 19.08779 22.03531 18.5792 12.03531 20.0792 [2,] 22.28053 26.56107 27.28053 20.91221 27.96469 25.4208 17.96469 26.9208 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 17.71947 20.0792 17.80725 20.43893 16.5792 15.43893 17.26336 21.21947 [2,] 22.28053 26.9208 24.19275 29.56107 23.4208 24.56107 22.73664 25.78053 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 14.53531 17.12309 20.89504 23.17557 17.5792 15.43893 21.5792 23.85973 [2,] 20.46469 24.87691 29.10496 26.82443 24.4208 24.56107 28.4208 26.14027 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [1,] 22.13168 23.63168 23.17557 16.12309 14.43893 24.31584 23.17557 22.49141 [2,] 24.86832 26.36832 26.82443 23.87691 23.56107 25.68416 26.82443 27.50859 [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 19.94752 20.99141 23.17557 21.21947 17.71947 19.71947 20.66699 18.71947 [2,] 24.05248 26.00859 26.82443 25.78053 22.28053 24.28053 29.33301 23.28053 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 21.21947 19.71947 19.71947 18.1584 20.89504 14.07061 11.93893 16.5792 [2,] 25.78053 24.28053 24.28053 31.8416 29.10496 25.92939 21.06107 23.4208 [,73] [,74] [,75] [,76] [1,] 16.93893 18.5792 20.43893 NA [2,] 26.06107 25.4208 29.56107 NA $out [1] 40 5 35 35 30 10 30 12 50 40 40 35 40 35 35 15 10 10 15 30 50 40 50 50 50 [26] 10 50 50 40 30 38 40 40 50 50 50 12 65 20 40 35 35 40 30 20 30 35 35 10 10 [51] 12 12 10 40 50 50 50 $group [1] 1 3 3 4 5 5 5 6 7 11 11 12 12 13 14 14 15 15 16 16 17 22 24 24 26 [26] 28 29 29 30 31 33 33 35 36 38 43 44 45 48 49 50 50 51 54 54 54 58 60 61 62 [51] 65 65 66 67 73 74 74 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" "73" "74" "75" [76] NA Warning message: In bxp(list(stats = c(15, 20, 25, 25, 30, 12, 15, 20, 25, 30, 15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6hg1m1322483116.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 5.982949 38 5 [2,] 6.342482 40 5 [3,] 6.660846 40 5 [4,] 7.525104 43 5 [5,] 8.374177 45 5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 6.121444 38.63168 5 [2,] 7.200247 41.36832 5 $out [1] 53 35 $group [1] 2 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(5.98294874419593, 6.34248174404548, 6.66084551359875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1p3bd1322483116.ps tmp/1p3bd1322483116.png",intern=TRUE)) character(0) > try(system("convert tmp/2ydom1322483116.ps tmp/2ydom1322483116.png",intern=TRUE)) character(0) > try(system("convert tmp/3c8fw1322483116.ps tmp/3c8fw1322483116.png",intern=TRUE)) character(0) > try(system("convert tmp/4ybls1322483116.ps tmp/4ybls1322483116.png",intern=TRUE)) character(0) > try(system("convert tmp/5h0ja1322483116.ps tmp/5h0ja1322483116.png",intern=TRUE)) character(0) > try(system("convert tmp/6hg1m1322483116.ps tmp/6hg1m1322483116.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.849 0.321 2.188