R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-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(8 + ,10 + ,10 + ,13 + ,14 + ,12 + ,11 + ,8 + ,8 + ,10 + ,10 + ,12 + ,12 + ,12 + ,11 + ,12 + ,12 + ,12 + ,12 + ,12 + ,12 + ,10 + ,12 + ,10 + ,11 + ,10 + ,10 + ,12 + ,10 + ,12 + ,7 + ,12 + ,18 + ,12 + ,11 + ,13 + ,10 + ,10 + ,10 + ,8 + ,12 + ,10 + ,10 + ,8 + ,14 + ,9 + ,8 + ,12 + ,15 + ,14 + ,1 + ,9 + ,7 + ,8 + ,12 + ,57 + ,12 + ,10 + ,10 + ,8 + ,8 + ,16 + ,14 + ,13 + ,10 + ,12 + ,9 + ,12 + ,11 + ,10 + ,8 + ,8 + ,9 + ,12 + ,8 + ,12 + ,10 + ,12 + ,9 + ,8 + ,12 + ,8 + ,12 + ,10 + ,12 + ,9 + ,28 + ,10 + ,12 + ,9 + ,14 + ,12 + ,12 + ,99 + ,13 + ,13 + ,14 + ,12 + ,12 + ,10 + ,11 + ,12 + ,14 + ,10 + ,12 + ,12 + ,6 + ,12 + ,10 + ,12 + ,12 + ,12 + ,9 + ,12 + ,12 + ,13 + ,8 + ,12 + ,10 + ,10 + ,10 + ,9 + ,12 + ,9 + ,10 + ,8 + ,12 + ,10 + ,8 + ,8 + ,9 + ,12 + ,12 + ,10 + ,10 + ,9 + ,11 + ,10 + ,9 + ,15 + ,10 + ,8 + ,10 + ,8 + ,9 + ,9 + ,6 + ,16 + ,12 + ,12 + ,12 + ,12 + ,10 + ,12 + ,8 + ,9 + ,12 + ,12 + ,8 + ,14 + ,10 + ,12 + ,8 + ,11 + ,10 + ,12 + ,12 + ,12 + ,12 + ,8 + ,10 + ,7 + ,10 + ,10 + ,12 + ,11 + ,9 + ,10 + ,12 + ,14 + ,13 + ,10 + ,11 + ,10 + ,10 + ,8 + ,10 + ,10 + ,10 + ,8 + ,8 + ,4 + ,14 + ,8 + ,12 + ,12 + ,10 + ,8 + ,12 + ,12 + ,10 + ,10 + ,12 + ,12 + ,9 + ,11 + ,14 + ,10 + ,8 + ,12 + ,8 + ,10 + ,11 + ,12 + ,10 + ,10 + ,12 + ,8 + ,9 + ,12 + ,8 + ,8 + ,10 + ,10 + ,10 + ,14 + ,10 + ,12 + ,12 + ,13 + ,9 + ,12 + ,12 + ,10 + ,12 + ,6 + ,8 + ,12 + ,10 + ,9 + ,11 + ,11 + ,9 + ,10 + ,15 + ,12 + ,7 + ,7 + ,10 + ,9 + ,10 + ,10 + ,9 + ,12 + ,10 + ,9 + ,12 + ,10 + ,7 + ,12 + ,10 + ,10 + ,12 + ,8 + ,12 + ,10 + ,10 + ,9 + ,8 + ,8 + ,12 + ,12 + ,10 + ,12 + ,10 + ,9 + ,10 + ,10 + ,8 + ,10 + ,12 + ,12 + ,16 + ,10 + ,9 + ,12 + ,12 + ,10 + ,7 + ,12 + ,10 + ,10 + ,6 + ,9 + ,6 + ,18 + ,13 + ,10 + ,12 + ,15 + ,12 + ,12 + ,9 + ,7 + ,12 + ,13 + ,14 + ,13 + ,12 + ,8 + ,8 + ,10 + ,10 + ,8 + ,12 + ,10 + ,12 + ,12 + ,12 + ,9 + ,12 + ,7 + ,12 + ,8 + ,8 + ,12 + ,14 + ,10 + ,5 + ,9 + ,8 + ,13 + ,10 + ,10 + ,14 + ,10 + ,99 + ,10 + ,12 + ,17 + ,14 + ,8 + ,14 + ,12 + ,12 + ,10 + ,10 + ,8 + ,12 + ,12 + ,12 + ,10 + ,12 + ,10 + ,10 + ,12 + ,12 + ,12 + ,12 + ,13 + ,12 + ,8 + ,10 + ,12 + ,8 + ,10 + ,10 + ,12 + ,12 + ,12 + ,12 + ,12 + ,12 + ,14 + ,10 + ,12 + ,14 + ,12 + ,14 + ,12 + ,13 + ,8 + ,12 + ,14 + ,10 + ,10 + ,11 + ,16 + ,12 + ,10 + ,10 + ,99 + ,8 + ,11 + ,12 + ,12 + ,11 + ,10 + ,20 + ,9 + ,14 + ,12 + ,10 + ,12 + ,10 + ,12 + ,12 + ,8 + ,12 + ,12 + ,10 + ,99 + ,12 + ,2 + ,10 + ,10 + ,10 + ,12 + ,12 + ,12 + ,12 + ,88 + ,9 + ,12 + ,14 + ,8 + ,12 + ,10 + ,10 + ,10 + ,7 + ,8 + ,10 + ,1 + ,10 + ,10 + ,9 + ,15 + ,10 + ,12 + ,12 + ,12 + ,11 + ,12 + ,12 + ,14 + ,8 + ,12 + ,12 + ,10 + ,14 + ,8 + ,10 + ,12 + ,10 + ,10 + ,10 + ,12 + ,9 + ,12 + ,11 + ,8 + ,14 + ,12 + ,10 + ,12 + ,10 + ,8 + ,14 + ,12 + ,12 + ,12 + ,8 + ,12 + ,12 + ,10 + ,12 + ,12 + ,12 + ,9 + ,11 + ,10 + ,15 + ,10 + ,9 + ,9 + ,10 + ,7 + ,10 + ,9 + ,10 + ,10 + ,10 + ,15 + ,12 + ,12 + ,10 + ,12 + ,8 + ,12 + ,11 + ,8 + ,14 + ,8 + ,12 + ,10 + ,15 + ,9 + ,13 + ,12 + ,14 + ,12 + ,12 + ,17 + ,10 + ,13 + ,12 + ,12 + ,10 + ,12 + ,10 + ,12 + ,10 + ,10 + ,10 + ,1 + ,8 + ,12 + ,10 + ,10 + ,10 + ,12 + ,12 + ,11 + ,12 + ,8 + ,8 + ,12 + ,12 + ,10 + ,12 + ,9 + ,10 + ,12 + ,12 + ,12 + ,12 + ,10 + ,10 + ,9 + ,12 + ,10 + ,9 + ,12 + ,7 + ,14 + ,10 + ,10 + ,9 + ,10 + ,8 + ,10 + ,12 + ,12 + ,10 + ,9 + ,10 + ,9 + ,12 + ,10 + ,12 + ,10 + ,9 + ,7 + ,12 + ,11 + ,12 + ,9 + ,13 + ,12 + ,12 + ,7 + ,8 + ,12 + ,12 + ,12 + ,11 + ,12 + ,13 + ,10 + ,12 + ,10 + ,12 + ,12 + ,15 + ,12 + ,12 + ,13 + ,10 + ,10 + ,8 + ,11 + ,12 + ,12 + ,12 + ,12 + ,10 + ,10 + ,12 + ,15 + ,12 + ,10 + ,10 + ,7 + ,12 + ,10 + ,11 + ,10 + ,10 + ,10 + ,10 + ,11 + ,7 + ,15 + ,8 + ,10 + ,6 + ,8 + ,9 + ,8 + ,7 + ,10 + ,12 + ,14 + ,11 + ,8 + ,10 + ,8 + ,8 + ,14 + ,12 + ,15 + ,12 + ,12 + ,9 + ,12 + ,12 + ,9 + ,11 + ,15 + ,11 + ,12 + ,7 + ,15 + ,9 + ,10 + ,15 + ,15 + ,8 + ,11 + ,12 + ,10 + ,10 + ,12 + ,7 + ,12 + ,10 + ,11 + ,12 + ,10 + ,10 + ,8 + ,9 + ,8 + ,10 + ,10 + ,14 + ,10 + ,10 + ,12 + ,12 + ,7 + ,12 + ,10 + ,12 + ,9 + ,9 + ,13 + ,14 + ,10 + ,12 + ,12 + ,12 + ,12 + ,12 + ,10 + ,10 + ,8 + ,12 + ,8 + ,14 + ,10 + ,70 + ,12 + ,10 + ,8 + ,8 + ,11 + ,10 + ,8 + ,7 + ,8 + ,50 + ,9 + ,12 + ,12 + ,7 + ,10 + ,8 + ,10 + ,10 + ,10 + ,8 + ,12 + ,7 + ,13 + ,13 + ,8 + ,8 + ,11 + ,6 + ,12 + ,9 + ,12 + ,13 + ,13 + ,12 + ,12 + ,10 + ,8 + ,12 + ,10 + ,10 + ,15 + ,12 + ,10 + ,12 + ,8 + ,8 + ,12 + ,12 + ,10 + ,12 + ,9 + ,12 + ,12 + ,10 + ,9 + ,10 + ,10 + ,10 + ,12 + ,12 + ,12 + ,12 + ,8 + ,10 + ,12 + ,15 + ,10 + ,8 + ,15 + ,10 + ,9 + ,12 + ,99 + ,10 + ,10 + ,11 + ,11 + ,12 + ,12 + ,14 + ,12 + ,14 + ,9 + ,10 + ,12 + ,13 + ,10 + ,11 + ,10 + ,12 + ,12 + ,12 + ,13 + ,14 + ,9 + ,10 + ,10 + ,12 + ,12 + ,10 + ,99 + ,8 + ,99 + ,12 + ,99 + ,12 + ,10 + ,12 + ,10 + ,12 + ,12 + ,12 + ,10 + ,12 + ,12 + ,10 + ,10 + ,12 + ,99 + ,12 + ,12 + ,9 + ,99 + ,12 + ,12 + ,9 + ,12 + ,12 + ,12 + ,15 + ,12 + ,12 + ,12 + ,8 + ,8 + ,12 + ,8 + ,12 + ,12 + ,10 + ,10 + ,12 + ,99 + ,8 + ,8 + ,99 + ,10 + ,10 + ,5 + ,9 + ,9 + ,99 + ,9 + ,10 + ,12) > par1 = '12' > par1 <- '12' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Standard Deviation Plot (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_sdplot.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > par1 <- as.numeric(par1) > (n <- length(x)) [1] 858 > (np <- floor(n / par1)) [1] 71 > 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] 72 72 72 72 72 72 71 71 71 71 71 71 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 8 12 11 10 15 8 9 12 14 10 10 12 9 [2,] 10 12 10 10 14 16 12 9 12 12 9 10 9 [3,] 10 11 10 10 1 14 8 28 12 12 12 10 6 [4,] 13 12 12 8 9 13 12 10 10 12 9 9 16 [5,] 14 12 10 12 7 10 10 12 11 9 10 11 12 [6,] 12 12 12 10 8 12 12 9 12 12 8 10 12 [7,] 11 12 7 10 12 9 9 14 14 12 12 9 12 [8,] 8 12 12 8 57 12 8 12 10 13 10 15 12 [9,] 8 12 18 14 12 11 12 12 12 8 8 10 10 [10,] 10 10 12 9 10 10 8 99 12 12 8 8 12 [11,] 10 12 11 8 10 8 12 13 6 10 9 10 8 [12,] 12 10 13 12 8 8 10 13 12 10 12 8 9 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 12 12 13 14 9 12 12 11 9 12 10 7 [2,] 12 8 10 8 11 8 13 11 12 10 10 12 [3,] 8 10 11 12 14 9 9 9 10 10 8 10 [4,] 14 7 10 12 10 12 12 10 9 9 10 10 [5,] 10 10 10 10 8 8 12 15 12 8 12 6 [6,] 12 10 8 8 12 8 10 12 10 8 12 9 [7,] 8 12 10 12 8 10 12 7 7 12 16 6 [8,] 11 11 10 12 10 10 6 7 12 12 10 18 [9,] 10 9 10 10 11 10 8 10 10 10 9 13 [10,] 12 10 8 10 12 14 12 9 10 12 12 10 [11,] 12 12 8 12 10 10 10 10 12 10 12 12 [12,] 12 14 4 12 10 12 9 10 8 9 10 15 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 12 10 8 99 12 12 12 10 11 12 12 10 [2,] 12 8 12 10 12 8 14 10 10 12 88 1 [3,] 9 12 14 12 12 10 10 11 20 10 9 10 [4,] 7 10 10 17 10 12 12 16 9 99 12 10 [5,] 12 12 5 14 12 8 14 12 14 12 14 9 [6,] 13 12 9 8 10 10 12 10 12 2 8 15 [7,] 14 12 8 14 10 10 14 10 10 10 12 10 [8,] 13 9 13 12 12 12 12 99 12 10 10 12 [9,] 12 12 10 12 12 12 13 8 10 10 10 12 [10,] 8 7 10 10 12 12 8 11 12 12 10 12 [11,] 8 12 14 10 12 12 12 12 12 12 7 11 [12,] 10 8 10 8 13 12 14 12 8 12 8 12 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 12 10 14 11 10 12 12 10 12 9 10 12 [2,] 14 12 12 10 15 10 12 10 9 12 9 9 [3,] 8 9 12 15 12 15 10 10 10 7 10 13 [4,] 12 12 12 10 12 9 12 12 12 14 9 12 [5,] 12 11 8 9 10 13 10 12 12 10 12 12 [6,] 10 8 12 9 12 12 12 11 12 10 10 7 [7,] 14 14 12 10 8 14 10 12 12 9 12 8 [8,] 8 12 10 7 12 12 10 8 10 10 10 12 [9,] 10 10 12 10 11 12 10 8 10 8 9 12 [10,] 12 12 12 9 8 17 1 12 9 10 7 12 [11,] 10 10 12 10 14 10 8 12 12 12 12 11 [12,] 10 8 9 10 8 13 12 10 10 12 11 12 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 13 8 10 8 10 11 12 9 12 10 10 10 [2,] 10 11 7 10 8 15 10 8 9 8 8 10 [3,] 12 12 12 6 8 11 10 10 9 12 7 8 [4,] 10 12 10 8 14 12 12 10 13 8 8 12 [5,] 12 12 11 9 12 7 7 14 14 14 50 7 [6,] 12 12 10 8 15 15 12 10 10 10 9 13 [7,] 15 10 10 7 12 9 10 10 12 70 12 13 [8,] 12 10 10 10 12 10 11 12 12 12 12 8 [9,] 12 12 10 12 9 15 12 12 12 10 7 8 [10,] 13 15 11 14 12 15 10 7 12 8 10 11 [11,] 10 12 7 11 12 8 10 12 12 8 8 6 [12,] 10 10 15 8 9 11 8 10 10 11 10 12 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 9 12 10 15 12 12 8 12 9 12 9 [2,] 12 10 9 10 12 12 99 12 12 12 9 [3,] 13 12 10 8 14 12 12 10 12 10 99 [4,] 13 8 10 15 12 13 99 10 12 10 9 [5,] 12 8 10 10 14 14 12 12 15 12 10 [6,] 12 12 12 9 9 9 10 99 12 99 12 [7,] 10 12 12 12 10 10 12 12 12 8 NA [8,] 8 10 12 99 12 10 10 12 12 8 NA [9,] 12 12 12 10 13 12 12 9 8 99 NA [10,] 10 9 8 10 10 12 12 99 8 10 NA [11,] 10 12 10 11 11 10 12 12 12 10 NA [12,] 15 12 12 11 10 99 10 12 8 5 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/1yr781386253282.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/2kink1386253282.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/3w6rc1386253282.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/4y8th1386253282.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,] 7 7 6 7 7 7.0 7 7 7 7 7 5 [2,] 10 9 9 10 10 9.0 10 10 10 9 10 9 [3,] 11 10 10 12 12 10.5 12 12 10 10 11 10 [4,] 12 12 12 12 12 12.0 12 12 12 12 12 12 [5,] 15 16 15 15 15 15.0 15 15 15 15 14 15 $n [1] 72 72 72 72 72 72 71 71 71 71 71 71 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 10.62759 9.441386 9.441386 11.62759 11.62759 9.941386 11.62498 11.62498 [2,] 11.37241 10.558614 10.558614 12.37241 12.37241 11.058614 12.37502 12.37502 [,9] [,10] [,11] [,12] [1,] 9.624977 9.437465 10.62498 9.437465 [2,] 10.375023 10.562535 11.37502 10.562535 $out [1] 99 88 1 99 1 28 20 99 16 17 16 99 99 6 5 50 2 99 99 16 6 70 57 6 18 [26] 99 99 18 99 99 17 1 99 6 6 4 99 $group [1] 1 2 2 2 3 3 3 3 4 4 4 4 4 5 5 5 6 6 6 7 7 7 8 8 8 [26] 8 8 9 9 10 10 10 10 11 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(7, 10, 11, 12, 15, 7, 9, 10, 12, 16, 6, 9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5uhmv1386253282.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,] 8 11.0 7.0 8.0 1 8.0 8.0 9.0 10.0 8 8.0 8.0 6 8 [2,] 9 11.5 10.0 8.5 8 8.5 8.5 11.0 10.5 10 8.5 9.0 9 10 [3,] 10 12.0 11.5 10.0 10 10.5 10.0 12.0 12.0 12 9.5 10.0 11 12 [4,] 12 12.0 12.0 11.0 13 12.5 12.0 13.5 12.0 12 11.0 10.5 12 12 [5,] 14 12.0 13.0 14.0 15 16.0 12.0 14.0 14.0 13 12.0 12.0 16 14 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 7.0 8 8 8.0 8.0 6 7 7 8 8 6.0 7.0 [2,] 9.5 8 10 9.5 8.5 9 9 9 9 10 8.0 8.5 [3,] 10.0 10 12 10.0 10.0 11 10 10 10 10 10.0 12.0 [4,] 12.0 10 12 11.5 12.0 12 11 12 12 12 12.5 12.5 [5,] 14.0 13 14 14.0 14.0 13 12 12 12 12 18.0 14.0 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 7.0 5.0 8 10 8 10 8 8.0 10 7.0 9.0 8 [2,] 8.5 8.5 10 11 10 12 10 10.0 10 8.5 10.0 10 [3,] 11.0 10.0 12 12 12 12 11 11.5 12 10.0 10.5 11 [4,] 12.0 12.5 14 12 12 14 12 12.0 12 12.0 12.0 12 [5,] 12.0 14.0 17 13 12 14 12 14.0 12 14.0 15.0 14 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 8.0 10 9 8.0 9.0 8 8.0 9 7 7.0 7 10.0 [2,] 9.5 11 9 9.0 11.0 10 10.0 10 9 9.0 10 10.0 [3,] 10.5 12 10 11.5 12.0 10 10.5 11 10 10.0 12 12.0 [4,] 12.0 12 10 12.0 13.5 12 12.0 12 12 11.5 12 12.5 [5,] 14.0 12 11 15.0 17.0 12 12.0 12 14 12.0 13 15.0 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] 8 10 6.0 8 7.0 7 7.0 9 8 7.0 6 8.0 [2,] 10 10 8.0 9 9.5 10 9.5 10 8 8.0 8 10.0 [3,] 12 10 8.5 12 11.0 10 10.0 12 10 9.5 10 12.0 [4,] 12 11 10.5 12 15.0 12 12.0 12 12 11.0 12 12.5 [5,] 15 12 14.0 15 15.0 12 14.0 14 14 12.0 13 15.0 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 8.0 8 8.0 9.0 9.0 8 10 8.0 5 9.0 [2,] 9.5 10 10.0 10.0 10.0 10 11 8.5 9 9.0 [3,] 12.0 10 10.5 12.0 12.0 12 12 12.0 10 9.5 [4,] 12.0 12 13.5 12.5 12.5 12 12 12.0 12 12.0 [5,] 12.0 12 15.0 14.0 14.0 12 12 15.0 12 12.0 $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 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 8.63168 11.77195 10.58779 8.859733 7.719466 8.675573 8.403627 [2,] 11.36832 12.22805 12.41221 11.140267 12.280534 12.324427 11.596373 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 10.85973 11.31584 11.08779 8.359733 9.31584 9.63168 11.08779 8.859733 [2,] 13.14027 12.68416 12.91221 10.640267 10.68416 12.36832 12.91221 11.140267 [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 9.087787 11.08779 9.087787 8.403627 9.63168 9.087787 8.63168 [2,] 10.912213 12.91221 10.912213 11.596373 12.36832 10.912213 11.36832 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [1,] 8.63168 9.087787 7.94752 10.17557 9.403627 8.175573 10.17557 11.54389 [2,] 11.36832 10.912213 12.05248 13.82443 12.596373 11.824427 13.82443 12.45611 [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 11.08779 11.08779 10.08779 10.58779 11.08779 8.403627 9.587787 10.08779 [2,] 12.91221 12.91221 11.91221 12.41221 12.91221 11.596373 11.412213 11.91221 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [1,] 9.359733 11.54389 9.543893 10.13168 10.85973 9.087787 9.587787 [2,] 11.640267 12.45611 10.456107 12.86832 13.14027 10.912213 11.412213 [,46] [,47] [,48] [,49] [,50] [,51] [,52] [,53] [1,] 10.08779 8.63168 8.859733 11.08779 10.85973 11.08779 9.543893 7.359733 [2,] 11.91221 11.36832 11.140267 12.91221 13.14027 12.91221 10.456107 9.640267 [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 10.63168 8.491413 9.087787 8.859733 11.08779 8.175573 8.13168 [2,] 13.36832 13.508587 10.912213 11.140267 12.91221 11.824427 10.86832 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [1,] 8.175573 10.85973 10.85973 9.087787 8.903627 10.85973 10.85973 11.08779 [2,] 11.824427 13.14027 13.14027 10.912213 12.096373 13.14027 13.14027 12.91221 [,69] [,70] [,71] [,72] [1,] 11.54389 10.40363 8.63168 7.564903 [2,] 12.45611 13.59637 11.36832 11.435097 $out [1] 10 10 18 57 28 99 6 15 4 15 16 99 8 16 99 20 99 2 88 1 14 8 9 15 7 [26] 1 7 7 15 70 50 99 99 99 99 99 9 99 99 99 99 $group [1] 2 2 3 5 8 8 9 12 16 21 24 29 32 33 33 34 35 35 36 37 40 40 40 41 41 [26] 44 52 52 52 59 60 65 67 68 68 69 69 69 71 71 72 $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" NA Warning message: In bxp(list(stats = c(8, 9, 10, 12, 14, 11, 11.5, 12, 12, 12, 7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6ry1m1386253282.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,] 1.786716 64.0 2 [2,] 8.933591 78.0 2 [3,] 10.809897 92.5 2 [4,] 14.741133 97.5 3 [5,] 15.699248 98.0 3 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8.161038 83.60592 1.543893 [2,] 13.458756 101.39408 2.456107 $out [1] 45 8 $group [1] 2 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(1.78671601681646, 8.93359086331201, 10.8098969130596, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1yr781386253282.ps tmp/1yr781386253282.png",intern=TRUE)) character(0) > try(system("convert tmp/2kink1386253282.ps tmp/2kink1386253282.png",intern=TRUE)) character(0) > try(system("convert tmp/3w6rc1386253282.ps tmp/3w6rc1386253282.png",intern=TRUE)) character(0) > try(system("convert tmp/4y8th1386253282.ps tmp/4y8th1386253282.png",intern=TRUE)) character(0) > try(system("convert tmp/5uhmv1386253282.ps tmp/5uhmv1386253282.png",intern=TRUE)) character(0) > try(system("convert tmp/6ry1m1386253282.ps tmp/6ry1m1386253282.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 7.164 1.979 9.140