R version 2.7.0 (2008-04-22) 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. 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 + ,2541086 + ,39690 + ,43129 + ,37863 + ,35953 + ,29133 + ,24693 + ,22205 + ,21725 + ,27192 + ,21790 + ,13253 + ,37702 + ,30364 + ,32609 + ,30212 + ,29965 + ,28352 + ,25814 + ,22414 + ,20506 + ,28806 + ,22228 + ,13971 + ,36845 + ,35338 + ,35022 + ,34777 + ,26887 + ,23970 + ,22780 + ,17351 + ,21382 + ,24561 + ,17409 + ,11514 + ,31514 + ,27071 + ,29462 + ,26105 + ,22397 + ,23843 + ,21705 + ,18089 + ,20764 + ,25316 + ,17704 + ,15548 + ,28029 + ,29383 + ,36438 + ,32034 + ,22679 + ,24319 + ,18004 + ,17537 + ,20366 + ,22782 + ,19169 + ,13807 + ,29743 + ,25591 + ,29096 + ,26482 + ,22405 + ,27044 + ,17970 + ,18730 + ,19684 + ,19785 + ,18479 + ,10698 + ,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] 971 > (np <- floor(n / par1)) [1] 80 > 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] 81 81 81 81 81 81 81 81 81 81 81 80 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 20 13253 13971 11514 15548 13807 10698 30 5 25 23 25 [2,] 2541086 37702 36845 31514 28029 29743 15 25 20 20 10 25 [3,] 39690 30364 35338 27071 29383 25591 15 20 15 20 22 22 [4,] 43129 32609 35022 29462 36438 29096 25 25 25 20 25 25 [5,] 37863 30212 34777 26105 32034 26482 25 20 22 25 25 12 [6,] 35953 29965 26887 22397 22679 22405 25 25 20 25 22 18 [7,] 29133 28352 23970 23843 24319 27044 21 20 22 15 30 20 [8,] 24693 25814 22780 21705 18004 17970 30 20 25 20 20 20 [9,] 22205 22414 17351 18089 17537 18730 25 15 20 35 25 22 [10,] 21725 20506 21382 20764 20366 19684 20 15 20 25 25 30 [11,] 27192 28806 24561 25316 22782 19785 40 12 35 25 22 25 [12,] 21790 22228 17409 17704 19169 18479 13 20 30 30 25 22 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 20 22 35 25 40 25 25 22 25 15 30 15 [2,] 50 25 22 35 20 25 20 23 25 22 20 16 [3,] 30 25 12 30 40 25 25 35 30 25 30 20 [4,] 25 25 30 12 25 35 30 15 22 25 50 20 [5,] 20 20 15 25 25 40 20 25 30 28 19 25 [6,] 30 22 10 15 20 20 22 18 10 22 20 30 [7,] 22 15 30 25 15 22 35 22 10 30 28 20 [8,] 25 20 9 25 15 22 22 25 25 25 20 25 [9,] 30 30 25 20 22 20 25 25 20 20 25 25 [10,] 22 20 20 20 24 25 25 28 22 25 35 25 [11,] 25 25 20 6 22 25 25 30 25 25 25 20 [12,] 22 30 35 15 20 18 25 20 25 20 25 20 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 25 30 25 22 15 25 20 22 25 22 25 22 [2,] 25 25 25 15 25 25 25 30 25 20 30 30 [3,] 30 20 20 22 20 25 25 30 25 20 25 25 [4,] 22 25 25 25 25 50 25 8 18 25 50 18 [5,] 20 20 15 12 25 19 20 20 25 20 30 22 [6,] 25 20 22 18 25 50 20 25 20 20 30 22 [7,] 25 20 25 30 20 25 20 30 25 20 50 30 [8,] 18 22 25 25 30 35 20 50 30 20 15 40 [9,] 18 18 15 25 20 20 25 22 18 25 25 25 [10,] 20 22 12 40 25 20 18 20 20 20 25 20 [11,] 25 20 25 24 30 20 25 10 25 10 22 10 [12,] 25 15 30 25 22 20 22 25 22 20 20 20 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 9 25 15 25 20 25 25 30 30 18 15 15 [2,] 15 20 10 30 25 20 25 20 22 30 22 15 [3,] 20 25 20 25 25 30 25 20 12 20 20 25 [4,] 15 25 25 10 25 50 20 25 12 25 10 25 [5,] 20 25 20 20 15 20 20 12 20 25 25 25 [6,] 30 30 20 25 40 25 20 15 20 25 20 15 [7,] 12 20 38 12 20 20 15 20 15 20 20 25 [8,] 15 25 20 15 20 30 20 20 25 20 15 22 [9,] 12 15 20 25 16 30 25 35 15 25 12 25 [10,] 20 15 20 20 25 25 15 22 20 20 20 20 [11,] 15 22 40 22 15 25 25 15 20 22 5 18 [12,] 12 10 25 22 20 12 50 18 25 15 20 22 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 25 24 15 20 30 25 25 25 25 15 15 20 [2,] 35 20 20 25 35 25 20 22 40 15 15 25 [3,] 25 25 20 25 25 22 22 20 25 25 18 25 [4,] 25 25 15 20 25 25 22 35 30 25 5 25 [5,] 25 25 35 25 25 30 20 20 25 20 15 22 [6,] 35 12 30 15 25 25 25 35 20 20 25 22 [7,] 30 15 20 20 25 25 22 25 25 15 18 25 [8,] 22 22 22 12 40 30 25 22 25 25 40 25 [9,] 30 25 65 15 40 25 22 25 30 15 25 30 [10,] 50 25 20 10 25 25 40 25 22 20 25 25 [11,] 15 25 25 25 25 30 25 25 22 22 20 25 [12,] 25 25 22 15 20 25 25 25 20 25 30 25 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 25 20 18 35 20 25 22 25 25 20 20 30 [2,] 20 20 15 20 25 20 25 10 16 25 25 30 [3,] 20 30 22 20 25 25 10 20 25 20 20 10 [4,] 25 25 22 20 25 25 20 25 35 20 22 22 [5,] 25 25 25 25 22 35 25 22 25 25 12 22 [6,] 25 30 22 25 25 22 20 22 15 22 30 25 [7,] 25 20 22 30 25 20 20 25 25 18 12 20 [8,] 20 25 25 15 20 25 30 25 25 30 22 22 [9,] 30 25 15 25 22 20 25 15 30 19 25 20 [10,] 25 24 20 22 30 20 20 25 25 25 25 25 [11,] 22 25 22 26 15 18 15 20 10 20 25 20 [12,] 30 30 18 25 30 20 20 10 22 25 25 15 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [1,] 25 20 30 28 18 20 15 22 50 [2,] 20 40 22 25 25 20 15 18 40 [3,] 25 20 30 22 22 25 30 50 25 [4,] 20 25 20 15 15 20 21 20 25 [5,] 30 20 25 40 25 12 25 50 25 [6,] 15 25 25 40 15 9 22 20 25 [7,] 40 20 30 12 15 15 22 20 30 [8,] 25 50 25 12 15 12 50 30 40 [9,] 20 50 25 18 25 15 15 25 25 [10,] 22 25 20 12 15 25 25 20 30 [11,] 22 25 18 25 12 20 15 22 20 [12,] 30 40 18 26 22 20 25 25 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/rcomp/tmp/1t6i51231599468.ps",horizontal=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/rcomp/tmp/2filj1231599468.ps",horizontal=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/rcomp/tmp/3twes1231599468.ps",horizontal=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/rcomp/tmp/49h8g1231599468.ps",horizontal=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 15 15 15 15 10 9 15 15 15 13.0 [2,] 20 20 20 20 20 20 20 20 20 20 20 20.0 [3,] 25 25 25 25 25 22 22 25 25 24 22 23.5 [4,] 25 25 25 25 25 25 28 30 25 25 25 25.0 [5,] 30 30 30 30 30 30 40 40 30 30 30 30.0 $n [1] 81 81 81 81 81 81 81 81 81 81 81 80 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 24.12222 24.12222 24.12222 24.12222 24.12222 21.12222 20.59556 23.24444 [2,] 25.87778 25.87778 25.87778 25.87778 25.87778 22.87778 23.40444 26.75556 [,9] [,10] [,11] [,12] [1,] 24.12222 23.12222 21.12222 22.61675 [2,] 25.87778 24.87778 22.87778 24.38325 $out [1] 13253 13971 11514 15548 13807 10698 5 35 40 [10] 9 35 50 2541086 37702 36845 31514 28029 29743 [19] 10 50 35 10 35 35 40 10 40 [28] 40 39690 30364 35338 27071 29383 25591 12 40 [37] 35 12 10 10 50 43129 32609 35022 29462 [46] 36438 29096 12 35 50 50 8 50 10 [55] 50 12 10 35 5 35 37863 30212 34777 [64] 26105 32034 26482 12 40 12 12 35 35 [73] 12 40 12 50 35953 29965 26887 22397 22679 [82] 22405 10 10 50 40 35 12 35 40 [91] 9 29133 28352 23970 23843 24319 27044 50 24693 [100] 25814 22780 21705 18004 17970 50 50 50 22205 [109] 22414 17351 18089 17537 18730 35 12 35 12 [118] 65 40 50 21725 20506 21382 20764 20366 19684 [127] 35 12 40 50 10 40 12 27192 28806 [136] 24561 25316 22782 19785 40 12 35 6 10 [145] 10 10 40 5 10 12 21790 22228 17409 [154] 17704 19169 18479 35 12 10 12 50 10 [163] 40 $group [1] 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 [26] 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4 4 [51] 4 4 4 4 4 4 4 4 4 4 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 [76] 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 7 7 7 7 7 7 7 8 8 [101] 8 8 8 8 8 8 8 9 9 9 9 9 9 9 9 9 9 9 9 9 10 10 10 10 10 [126] 10 10 10 10 10 10 10 10 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 [151] 12 12 12 12 12 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(15, 20, 25, 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/rcomp/tmp/5bb8o1231599468.ps",horizontal=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] [1,] 20.0 13253 13971.0 11514.0 15548.0 13807.0 13.0 12.0 15 15 20 [2,] 21997.5 22321 19395.5 19426.5 18586.5 18604.5 17.5 17.5 20 20 22 [3,] 28162.5 28579 24265.5 23120.0 22730.5 21095.0 25.0 20.0 21 25 24 [4,] 38776.5 30288 34899.5 26588.0 28706.0 26763.0 27.5 25.0 25 25 25 [5,] 43129.0 37702 36845.0 31514.0 36438.0 29743.0 40.0 30.0 30 30 25 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 18 20 15.0 9.0 6.0 15 18 20 15.0 20 15 19 [2,] 20 22 20.0 13.5 15.0 20 21 22 21.0 21 21 20 [3,] 22 25 23.5 21.0 22.5 22 25 25 24.0 25 25 25 [4,] 25 30 25.0 30.0 25.0 25 25 25 26.5 25 25 30 [5,] 30 30 30.0 35.0 35.0 25 25 25 30.0 30 30 35 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] 15 18 15.0 12.0 15.0 15 19 18 8.0 18 20 15.0 [2,] 20 20 20.0 17.5 20.0 20 20 20 20.0 20 20 23.5 [3,] 20 25 20.0 25.0 24.5 25 25 21 23.5 25 20 25.0 [4,] 25 25 23.5 25.0 25.0 25 30 25 30.0 25 21 30.0 [5,] 30 30 25.0 30.0 30.0 30 35 25 30.0 30 22 30.0 [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [1,] 10.0 9 10.0 15 10.0 15 12 15.0 12.0 12.0 15 5.0 [2,] 20.0 12 17.5 20 17.5 18 20 20.0 16.5 15.0 20 13.5 [3,] 22.0 15 23.5 20 22.0 20 25 22.5 20.0 20.0 21 20.0 [4,] 27.5 20 25.0 25 25.0 25 30 25.0 23.5 23.5 25 20.0 [5,] 30.0 30 30.0 25 30.0 25 30 25.0 30.0 30.0 30 25.0 [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [1,] 15.0 15.0 15 15.0 10 20.0 22.0 20.0 20 20.0 15 5 [2,] 16.5 25.0 21 20.0 15 25.0 25.0 22.0 22 22.0 15 15 [3,] 22.0 25.0 25 21.0 20 25.0 25.0 23.5 25 25.0 20 19 [4,] 25.0 32.5 25 27.5 25 32.5 27.5 25.0 25 27.5 25 25 [5,] 25.0 35.0 25 35.0 25 40.0 30.0 25.0 25 30.0 25 40 [,60] [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [1,] 22.0 20 20.0 15 15.0 15 18 15 10.0 10 18 20.0 [2,] 23.5 21 22.0 18 20.0 21 20 20 17.5 19 20 20.0 [3,] 25.0 25 25.0 22 25.0 25 21 20 22.0 25 21 23.5 [4,] 25.0 25 27.5 22 25.5 25 25 25 25.0 25 25 25.0 [5,] 25.0 30 30.0 25 30.0 30 25 30 25.0 30 30 30.0 [,72] [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [1,] 15 15.0 20 18.0 12.0 12.0 9.0 15 18.0 20 [2,] 20 20.0 20 20.0 13.5 15.0 13.5 15 20.0 25 [3,] 22 23.5 25 25.0 23.5 16.5 20.0 22 22.0 25 [4,] 25 27.5 40 27.5 27.0 23.5 20.0 25 27.5 35 [5,] 30 30.0 50 30.0 40.0 25.0 25.0 30 30.0 50 $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] 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 20509.49 24945.2 17194.02 19853.59 18114.93 17373.85 20.43893 16.5792 [2,] 35815.51 32212.8 31336.98 26386.41 27346.07 24816.15 29.56107 23.4208 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 18.71947 22.71947 22.63168 19.71947 21.35115 21.21947 13.47424 17.93893 [2,] 23.28053 27.28053 25.36832 24.28053 28.64885 25.78053 28.52576 27.06107 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 19.71947 23.17557 23.63168 21.49141 23.17557 23.17557 20.43893 17.71947 [2,] 24.28053 26.82443 26.36832 26.50859 26.82443 26.82443 29.56107 22.28053 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 22.71947 18.40363 21.5792 22.21947 22.71947 20.43893 18.71947 18.93893 [2,] 27.28053 21.59637 28.4208 26.78053 27.28053 29.56107 23.28053 28.06107 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 22.71947 19.54389 22.03531 18.5792 11.35115 20.0792 17.71947 18.5792 [2,] 27.28053 20.45611 27.96469 25.4208 18.64885 26.9208 22.28053 25.4208 [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 16.80725 20.43893 20.21947 16.80725 16.12309 18.71947 17.03531 18.12309 [2,] 23.19275 29.56107 24.78053 23.19275 23.87691 23.28053 22.96469 25.87691 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [1,] 21.5792 23.17557 17.5792 15.43893 21.5792 23.85973 22.13168 23.63168 [2,] 28.4208 26.82443 24.4208 24.56107 28.4208 26.14027 24.86832 26.36832 [,57] [,58] [,59] [,60] [,61] [,62] [,63] [,64] [1,] 22.49141 15.43893 14.43893 24.31584 23.17557 22.49141 20.17557 22.49141 [2,] 27.50859 24.56107 23.56107 25.68416 26.82443 27.50859 23.82443 27.50859 [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [1,] 23.17557 18.71947 17.71947 18.5792 22.26336 18.71947 21.21947 19.71947 [2,] 26.82443 23.28053 22.28053 25.4208 27.73664 23.28053 25.78053 24.28053 [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [1,] 20.0792 15.87787 21.5792 17.34256 12.62309 17.03531 17.43893 18.5792 [2,] 26.9208 34.12213 28.4208 29.65744 20.37691 22.96469 26.56107 25.4208 [,81] [1,] 20.23612 [2,] 29.76388 $out [1] 2541086 10698 5 35 35 10 30 12 50 [10] 40 40 35 40 30 35 35 10 10 [19] 50 30 12 40 50 50 50 25 25 [28] 10 50 50 40 10 38 40 40 50 [37] 50 35 50 12 65 40 35 35 40 [46] 20 30 35 35 10 35 12 12 10 [55] 40 50 50 50 $group [1] 1 7 9 9 10 11 11 12 13 17 17 18 18 19 19 20 21 21 23 26 28 28 30 30 32 [26] 34 34 34 35 35 36 39 39 39 41 42 43 44 49 50 51 55 56 56 57 60 60 64 66 67 [51] 69 71 71 72 73 79 80 80 $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] "76" "77" "78" "79" "80" NA Warning message: In bxp(list(stats = c(20, 21997.5, 28162.5, 38776.5, 43129, 13253, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6tgy81231599468.ps",horizontal=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,] 3482.482 15543.0 5 [2,] 5320.455 22310.0 5 [3,] 6731.528 28962.0 5 [4,] 8320.906 38765.5 5 [5,] 9124.545 43124.0 5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 5363.002 21456.54 5 [2,] 8100.054 36467.46 5 $out [1] 282225.9 2541076.0 8.0 10.0 $group [1] 1 2 3 3 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(3482.48166094628, 5320.45515881554, 6731.52800868184, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1t6i51231599468.ps tmp/1t6i51231599468.png") > system("convert tmp/2filj1231599468.ps tmp/2filj1231599468.png") > system("convert tmp/3twes1231599468.ps tmp/3twes1231599468.png") > system("convert tmp/49h8g1231599468.ps tmp/49h8g1231599468.png") > system("convert tmp/5bb8o1231599468.ps tmp/5bb8o1231599468.png") > system("convert tmp/6tgy81231599468.ps tmp/6tgy81231599468.png") > > > proc.time() user system elapsed 2.829 1.404 3.088