R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(356445,291705,310900,332340,257166,334551,317365,270863,317904,423141,317684,411063,371161,299023,326964,327146,303447,351994,320317,257151,320274,476982,301723,363567,338831,265802,307691,334207,303127,318863,292123,245155,284794,391604,304982,369552,356021,247577,277885,294032,310845,311023,298462,234188,297478,371017,291128,316374,326001,222302,227424,255428,278250,280335,241894,255075,255115,319482,270694,300209,283531,218924,236466,267980,219994,256052,230444,200778,240960,277837,209776,232065) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 356445 371161 338831 356021 326001 283531 NA [2,] 291705 299023 265802 247577 222302 218924 NA [3,] 310900 326964 307691 277885 227424 236466 NA [4,] 332340 327146 334207 294032 255428 267980 NA [5,] 257166 303447 303127 310845 278250 219994 NA [6,] 334551 351994 318863 311023 280335 256052 NA [7,] 317365 320317 292123 298462 241894 230444 NA [8,] 270863 257151 245155 234188 255075 200778 NA [9,] 317904 320274 284794 297478 255115 240960 NA [10,] 423141 476982 391604 371017 319482 277837 NA [11,] 317684 301723 304982 291128 270694 209776 NA [12,] 411063 363567 369552 316374 300209 232065 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/1ddpg1243852716.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/2jh5w1243852716.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/38udl1243852716.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/4q8rp1243852716.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] [1,] 283531 218924.0 227424 255428 219994.0 256052 230444.0 200778 240960 [2,] 326001 222302.0 236466 267980 257166.0 280335 241894.0 234188 255115 [3,] 347426 256689.5 292788 310589 290688.5 314943 295292.5 250115 291136 [4,] 356445 291705.0 310900 332340 303447.0 334551 317365.0 257151 317904 [5,] 371161 299023.0 326964 334207 310845.0 351994 320317.0 270863 320274 [,10] [,11] [,12] [1,] 277837.0 270694.0 232065.0 [2,] 319482.0 270694.0 300209.0 [3,] 381310.5 296425.5 339970.5 [4,] 423141.0 304982.0 369552.0 [5,] 476982.0 317684.0 411063.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 327788.6 211922.3 244775.7 269074.7 260835.8 279971.9 246611.3 235303.1 [2,] 367063.4 301456.7 340800.3 352103.3 320541.2 349914.1 343973.7 264926.9 [,9] [,10] [,11] [,12] [1,] 250635.1 314447.1 274308.6 295242 [2,] 331636.9 448173.9 318542.4 384699 $out [1] 209776 $group [1] 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(283531, 326001, 347426, 356445, 371161, 218924, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5aisq1243852716.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] [1,] 257166.0 257151.0 245155.0 247577.0 222302.0 200778.0 NA [2,] 301302.5 302585.0 288458.5 284506.5 248484.5 219459.0 NA [3,] 317794.0 323640.5 306336.5 297970.0 263061.0 234265.5 NA [4,] 345498.0 357780.5 336519.0 313698.5 290272.0 262016.0 NA [5,] 411063.0 371161.0 391604.0 356021.0 326001.0 283531.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 297636.1 298465.5 284415.8 284655.3 244001.4 214855.0 NA [2,] 337951.9 348815.5 328257.2 311284.7 282120.6 253676.0 NA $out [1] 423141 476982 234188 371017 $group [1] 1 2 4 4 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(257166, 301302.5, 317794, 345498, 411063, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6z5vs1243852716.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,] 24424.87 70085.0 19702.25 [2,] 33303.95 79706.5 34647.38 [3,] 35185.47 90362.0 53405.62 [4,] 40228.17 103724.0 60732.50 [5,] 41429.90 107908.0 82891.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 32027.28 79407.46 41508.02 [2,] 38343.66 101316.54 65303.23 $out [1] 71435.93 63086.16 199145.00 178998.00 $group [1] 1 1 2 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(24424.8718782037, 33303.9478780668, 35185.469871993, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ddpg1243852716.ps tmp/1ddpg1243852716.png") > system("convert tmp/2jh5w1243852716.ps tmp/2jh5w1243852716.png") > system("convert tmp/38udl1243852716.ps tmp/38udl1243852716.png") > system("convert tmp/4q8rp1243852716.ps tmp/4q8rp1243852716.png") > system("convert tmp/5aisq1243852716.ps tmp/5aisq1243852716.png") > system("convert tmp/6z5vs1243852716.ps tmp/6z5vs1243852716.png") > > > proc.time() user system elapsed 1.053 0.826 1.282