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. Natural language support but running in an English locale 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/www/html/freestat/rcomp/tmp/1y6u61291560928.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/www/html/freestat/rcomp/tmp/2y6u61291560928.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/www/html/freestat/rcomp/tmp/39xur1291560928.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/www/html/freestat/rcomp/tmp/49xur1291560928.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/www/html/freestat/rcomp/tmp/5jotu1291560928.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/www/html/freestat/rcomp/tmp/6jotu1291560928.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/1y6u61291560928.ps tmp/1y6u61291560928.png",intern=TRUE)) character(0) > try(system("convert tmp/2y6u61291560928.ps tmp/2y6u61291560928.png",intern=TRUE)) character(0) > try(system("convert tmp/39xur1291560928.ps tmp/39xur1291560928.png",intern=TRUE)) character(0) > try(system("convert tmp/49xur1291560928.ps tmp/49xur1291560928.png",intern=TRUE)) character(0) > try(system("convert tmp/5jotu1291560928.ps tmp/5jotu1291560928.png",intern=TRUE)) character(0) > try(system("convert tmp/6jotu1291560928.ps tmp/6jotu1291560928.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.062 1.308 2.281