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(5 + ,5 + ,5 + ,6 + ,8 + ,9 + ,9 + ,9 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,10 + ,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 + ,13 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,15 + ,16 + ,16 + ,16 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,18 + ,19 + ,19 + ,19 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,20 + ,21 + ,21 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,22 + ,23 + ,23 + ,24 + ,24 + ,24 + ,24 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,25 + ,26 + ,26 + ,28 + ,28 + ,28 + ,28 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,30 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,35 + ,38 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,40 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,50 + ,65) > 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,] 5 10 10 12 12 15 15 15 15 15 15 18 18 [2,] 5 10 12 12 12 15 15 15 15 15 15 18 18 [3,] 5 10 12 12 13 15 15 15 15 15 15 18 18 [4,] 6 10 12 12 15 15 15 15 15 15 15 18 18 [5,] 8 10 12 12 15 15 15 15 15 15 15 18 18 [6,] 9 10 12 12 15 15 15 15 15 15 15 18 18 [7,] 9 10 12 12 15 15 15 15 15 15 15 18 18 [8,] 9 10 12 12 15 15 15 15 15 15 15 18 18 [9,] 10 10 12 12 15 15 15 15 15 15 15 18 18 [10,] 10 10 12 12 15 15 15 15 15 15 16 18 18 [11,] 10 10 12 12 15 15 15 15 15 15 16 18 18 [12,] 10 10 12 12 15 15 15 15 15 15 16 18 18 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 18 20 20 20 20 20 20 20 20 20 20 20 [2,] 19 20 20 20 20 20 20 20 20 20 20 20 [3,] 19 20 20 20 20 20 20 20 20 20 20 20 [4,] 19 20 20 20 20 20 20 20 20 20 20 20 [5,] 20 20 20 20 20 20 20 20 20 20 20 20 [6,] 20 20 20 20 20 20 20 20 20 20 20 20 [7,] 20 20 20 20 20 20 20 20 20 20 20 20 [8,] 20 20 20 20 20 20 20 20 20 20 20 20 [9,] 20 20 20 20 20 20 20 20 20 20 20 20 [10,] 20 20 20 20 20 20 20 20 20 20 20 20 [11,] 20 20 20 20 20 20 20 20 20 20 20 20 [12,] 20 20 20 20 20 20 20 20 20 20 20 20 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 20 20 20 20 20 22 22 22 22 22 22 22 [2,] 20 20 20 20 20 22 22 22 22 22 22 22 [3,] 20 20 20 20 20 22 22 22 22 22 22 22 [4,] 20 20 20 20 20 22 22 22 22 22 22 22 [5,] 20 20 20 20 20 22 22 22 22 22 22 22 [6,] 20 20 20 20 20 22 22 22 22 22 22 22 [7,] 20 20 20 20 20 22 22 22 22 22 22 22 [8,] 20 20 20 20 21 22 22 22 22 22 22 22 [9,] 20 20 20 20 21 22 22 22 22 22 22 22 [10,] 20 20 20 20 22 22 22 22 22 22 22 22 [11,] 20 20 20 20 22 22 22 22 22 22 22 22 [12,] 20 20 20 20 22 22 22 22 22 22 22 22 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 22 24 25 25 25 25 25 25 25 25 25 25 [2,] 22 24 25 25 25 25 25 25 25 25 25 25 [3,] 22 24 25 25 25 25 25 25 25 25 25 25 [4,] 22 24 25 25 25 25 25 25 25 25 25 25 [5,] 22 25 25 25 25 25 25 25 25 25 25 25 [6,] 22 25 25 25 25 25 25 25 25 25 25 25 [7,] 22 25 25 25 25 25 25 25 25 25 25 25 [8,] 22 25 25 25 25 25 25 25 25 25 25 25 [9,] 22 25 25 25 25 25 25 25 25 25 25 25 [10,] 22 25 25 25 25 25 25 25 25 25 25 25 [11,] 23 25 25 25 25 25 25 25 25 25 25 25 [12,] 23 25 25 25 25 25 25 25 25 25 25 25 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 25 25 25 25 25 25 25 25 25 25 25 25 [2,] 25 25 25 25 25 25 25 25 25 25 25 25 [3,] 25 25 25 25 25 25 25 25 25 25 25 25 [4,] 25 25 25 25 25 25 25 25 25 25 25 25 [5,] 25 25 25 25 25 25 25 25 25 25 25 25 [6,] 25 25 25 25 25 25 25 25 25 25 25 25 [7,] 25 25 25 25 25 25 25 25 25 25 25 25 [8,] 25 25 25 25 25 25 25 25 25 25 25 25 [9,] 25 25 25 25 25 25 25 25 25 25 25 25 [10,] 25 25 25 25 25 25 25 25 25 25 25 25 [11,] 25 25 25 25 25 25 25 25 25 25 25 25 [12,] 25 25 25 25 25 25 25 25 25 25 25 25 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 25 25 26 30 30 30 30 30 30 30 35 40 [2,] 25 25 26 30 30 30 30 30 30 30 35 40 [3,] 25 25 28 30 30 30 30 30 30 35 35 40 [4,] 25 25 28 30 30 30 30 30 30 35 35 40 [5,] 25 25 28 30 30 30 30 30 30 35 35 40 [6,] 25 25 28 30 30 30 30 30 30 35 35 40 [7,] 25 25 30 30 30 30 30 30 30 35 35 40 [8,] 25 25 30 30 30 30 30 30 30 35 35 40 [9,] 25 25 30 30 30 30 30 30 30 35 35 40 [10,] 25 25 30 30 30 30 30 30 30 35 35 40 [11,] 25 25 30 30 30 30 30 30 30 35 38 40 [12,] 25 25 30 30 30 30 30 30 30 35 40 40 [,74] [,75] [,76] [1,] 40 50 NA [2,] 40 50 NA [3,] 40 50 NA [4,] 40 50 NA [5,] 40 50 NA [6,] 40 50 NA [7,] 40 50 NA [8,] 50 50 NA [9,] 50 50 NA [10,] 50 50 NA [11,] 50 50 NA [12,] 50 65 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/1vra31322496505.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/2dqpv1322496505.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/34jl01322496505.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/4ow5r1322496505.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,] 15 15 13 15 15 15 15 15 15 15 15 15 [2,] 20 20 20 20 20 20 20 20 20 20 20 20 [3,] 22 22 22 22 22 22 22 22 22 22 23 23 [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 21.08779 21.08779 21.08779 21.08779 21.08779 21.08779 21.08779 [2,] 22.91221 22.91221 22.91221 22.91221 22.91221 22.91221 22.91221 22.91221 [,9] [,10] [,11] [,12] [1,] 21.08779 21.08779 22.08779 22.08779 [2,] 22.91221 22.91221 23.91221 23.91221 $out [1] 5 10 10 12 12 35 40 40 50 5 10 12 12 12 35 40 40 50 5 10 12 12 35 35 40 [26] 40 50 6 10 12 12 35 35 40 40 50 8 10 12 12 35 35 40 40 50 9 10 12 12 35 [51] 35 40 40 50 9 10 12 12 35 35 40 40 50 9 10 12 12 35 35 40 50 50 10 10 12 [76] 12 35 35 40 50 50 10 10 12 12 35 35 40 50 50 10 10 12 12 35 38 40 50 50 10 [101] 10 12 12 35 40 40 50 65 $group [1] 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 [26] 3 3 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 5 6 6 6 6 6 [51] 6 6 6 6 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 8 9 9 9 [76] 9 9 9 9 9 9 10 10 10 10 10 10 10 10 10 11 11 11 11 11 11 11 11 11 12 [101] 12 12 12 12 12 12 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5ipju1322496505.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,] 5.0 10 12 12 13 15 15 15 15 15 15.0 18 18 18 [2,] 5.5 10 12 12 14 15 15 15 15 15 15.0 18 18 19 [3,] 9.0 10 12 12 15 15 15 15 15 15 15.0 18 18 20 [4,] 10.0 10 12 12 15 15 15 15 15 15 15.5 18 18 20 [5,] 10.0 10 12 12 15 15 15 15 15 15 16.0 18 18 20 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 20 20 20 20 20 20 20 20 20 20 20 20 [2,] 20 20 20 20 20 20 20 20 20 20 20 20 [3,] 20 20 20 20 20 20 20 20 20 20 20 20 [4,] 20 20 20 20 20 20 20 20 20 20 20 20 [5,] 20 20 20 20 20 20 20 20 20 20 20 20 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 20 20 20 20.0 22 22 22 22 22 22 22 22 [2,] 20 20 20 20.0 22 22 22 22 22 22 22 22 [3,] 20 20 20 20.0 22 22 22 22 22 22 22 22 [4,] 20 20 20 21.5 22 22 22 22 22 22 22 22 [5,] 20 20 20 22.0 22 22 22 22 22 22 22 22 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 24 25 25 25 25 25 25 25 25 25 25 25 [2,] 24 25 25 25 25 25 25 25 25 25 25 25 [3,] 25 25 25 25 25 25 25 25 25 25 25 25 [4,] 25 25 25 25 25 25 25 25 25 25 25 25 [5,] 25 25 25 25 25 25 25 25 25 25 25 25 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 25 25 25 25 25 25 25 25 25 25 25 25 [2,] 25 25 25 25 25 25 25 25 25 25 25 25 [3,] 25 25 25 25 25 25 25 25 25 25 25 25 [4,] 25 25 25 25 25 25 25 25 25 25 25 25 [5,] 25 25 25 25 25 25 25 25 25 25 25 25 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] 25 26 30 30 30 30 30 30 35 35 40 40 [2,] 25 28 30 30 30 30 30 30 35 35 40 40 [3,] 25 29 30 30 30 30 30 30 35 35 40 40 [4,] 25 30 30 30 30 30 30 30 35 35 40 50 [5,] 25 30 30 30 30 30 30 30 35 35 40 50 [,75] [,76] [1,] 50 NA [2,] 50 NA [3,] 50 NA [4,] 50 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] [,9] [,10] [,11] [,12] [1,] 6.94752 10 12 12 14.54389 15 15 15 15 15 14.77195 18 [2,] 11.05248 10 12 12 15.45611 15 15 15 15 15 15.22805 18 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 18 19.54389 20 20 20 20 20 20 20 20 20 20 [2,] 18 20.45611 20 20 20 20 20 20 20 20 20 20 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 20 20 20 20 20 19.31584 22 22 22 22 22 22 [2,] 20 20 20 20 20 20.68416 22 22 22 22 22 22 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 22 22 24.54389 25 25 25 25 25 25 25 25 25 [2,] 22 22 25.45611 25 25 25 25 25 25 25 25 25 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 25 25 25 25 25 25 25 25 25 25 25 25 [2,] 25 25 25 25 25 25 25 25 25 25 25 25 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 25 25 25 28.08779 30 30 30 30 30 30 35 35 [2,] 25 25 25 29.91221 30 30 30 30 30 30 35 35 [,73] [,74] [,75] [,76] [1,] 40 35.43893 50 NA [2,] 40 44.56107 50 NA $out [1] 10 12 12 23 23 30 30 38 40 65 $group [1] 3 5 5 38 38 71 71 72 72 75 $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(5, 5.5, 9, 10, 10, 10, 10, 10, 10, 10, 12, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6qem81322496505.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,] 6.656785 40.0 5 [2,] 6.710925 40.5 5 [3,] 6.808849 41.5 5 [4,] 7.089506 45.0 5 [5,] 7.131531 45.0 5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 6.636176 39.44752 5 [2,] 6.981522 43.55248 5 $out [1] 8.092372 55.000000 $group [1] 1 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(6.65678546993913, 6.71092494350862, 6.808849077895, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1vra31322496505.ps tmp/1vra31322496505.png",intern=TRUE)) character(0) > try(system("convert tmp/2dqpv1322496505.ps tmp/2dqpv1322496505.png",intern=TRUE)) character(0) > try(system("convert tmp/34jl01322496505.ps tmp/34jl01322496505.png",intern=TRUE)) character(0) > try(system("convert tmp/4ow5r1322496505.ps tmp/4ow5r1322496505.png",intern=TRUE)) character(0) > try(system("convert tmp/5ipju1322496505.ps tmp/5ipju1322496505.png",intern=TRUE)) character(0) > try(system("convert tmp/6qem81322496505.ps tmp/6qem81322496505.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.652 0.310 1.981