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(274412,272433,268361,268586,264768,269974,304744,309365,308347,298427,289231,291975,294912,293488,290555,284736,281818,287854,316263,325412,326011,328282,317480,317539,313737,312276,309391,302950,300316,304035,333476,337698,335932,323931,313927,314485,313218,309664,302963,298989,298423,301631,329765,335083,327616,309119,295916,291413,291542,284678,276475,272566,264981,263290,296806,303598,286994,276427,266424,267153,268381,262522,255542,253158,243803,250741,280445,285257,270976,261076,255603,260376,263903,264291,263276,262572,256167,264221,293860)
> 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] 79
> (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] 7 7 7 7 7 7 7 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 274412 294912 313737 313218 291542 268381 263903
 [2,] 272433 293488 312276 309664 284678 262522 264291
 [3,] 268361 290555 309391 302963 276475 255542 263276
 [4,] 268586 284736 302950 298989 272566 253158 262572
 [5,] 264768 281818 300316 298423 264981 243803 256167
 [6,] 269974 287854 304035 301631 263290 250741 264221
 [7,] 304744 316263 333476 329765 296806 280445 293860
 [8,] 309365 325412 337698 335083 303598 285257     NA
 [9,] 308347 326011 335932 327616 286994 270976     NA
[10,] 298427 328282 323931 309119 276427 261076     NA
[11,] 289231 317480 313927 295916 266424 255603     NA
[12,] 291975 317539 314485 291413 267153 260376     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/1c6171257969391.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/2unf41257969391.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/3bw6f1257969391.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/4lhk01257969391.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,] 263903.0 262522 255542.0 253158.0 243803.0 250741.0 280445 285257.0 270976
[2,] 271396.5 268362 265818.5 265579.0 260467.5 263755.5 295333 303598.0 286994
[3,] 291542.0 284678 276475.0 272566.0 264981.0 269974.0 304744 317388.5 317179
[4,] 304065.0 301576 296759.0 291862.5 290120.5 294742.5 323014 335083.0 327616
[5,] 313737.0 312276 309391.0 302950.0 300316.0 304035.0 333476 337698.0 335932
      [,10]    [,11]  [,12]
[1,] 261076 255603.0 260376
[2,] 276427 266424.0 267153
[3,] 303773 292573.5 291694
[4,] 323931 313927.0 314485
[5,] 328282 317480.0 317539

$n
 [1] 7 7 7 7 7 7 7 6 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 272032.9 264843.1 257997.8 256869.9 247272.7 251469.1 288213.4 297079.7
[2,] 311051.1 304512.9 294952.2 288262.1 282689.3 288478.9 321274.6 337697.3
         [,9]    [,10]    [,11]    [,12]
[1,] 290976.5 273131.4 261932.5 261163.3
[2,] 343381.5 334414.6 323214.5 322224.7

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(263903, 271396.5, 291542, 304065, 313737,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5r4nl1257969391.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,] 264768.0 281818.0 300316.0 291413 263290.0 243803.0 262572
[2,] 269280.0 289204.5 306713.0 298706 266788.5 254350.0 262924
[3,] 281821.5 305587.5 313832.0 306041 276451.0 260726.0 263903
[4,] 301585.5 321475.5 328703.5 320417 289268.0 269678.5 264256
[5,] 309365.0 328282.0 337698.0 335083 303598.0 285257.0 264291

$n
[1] 12 12 12 12 12 12  7

$conf
         [,1]     [,2]   [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 267086.7 290868.5 303802 296138.5 266197.9 253734.6 263107.6
[2,] 296556.3 320306.5 323862 315943.5 286704.1 267717.4 264698.4

$out
[1] 256167 293860

$group
[1] 7 7

$names
[1] "1" "2" "3" "4" "5" "6" NA 

Warning message:
In bxp(list(stats = c(264768, 269280, 281821.5, 301585.5, 309365,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6pdqt1257969391.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,] 18622.84 49754.0 26283.50
[2,] 20384.24 51137.5 28667.00
[3,] 20547.94 53571.5 31827.75
[4,] 24205.15 59520.0 35261.00
[5,] 26562.04 67206.0 38301.00

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 18805.20 49748.19 28820.18
[2,] 22290.68 57394.81 34835.32

$out
numeric(0)

$group
numeric(0)

$names
[1] "S.D."  "Range" "IQR"  

Warning message:
In bxp(list(stats = c(18622.8426491052, 20384.2404238633, 20547.9408868397,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/1c6171257969391.ps tmp/1c6171257969391.png")
> system("convert tmp/2unf41257969391.ps tmp/2unf41257969391.png")
> system("convert tmp/3bw6f1257969391.ps tmp/3bw6f1257969391.png")
> system("convert tmp/4lhk01257969391.ps tmp/4lhk01257969391.png")
> system("convert tmp/5r4nl1257969391.ps tmp/5r4nl1257969391.png")
> system("convert tmp/6pdqt1257969391.ps tmp/6pdqt1257969391.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.087   0.835   1.528