R version 2.8.1 (2008-12-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(325412,326011,328282,317480,317539,313737,312276,309391,302950,300316,304035,333476,337698,335932,323931,313927,314485,313218,309664,302963,298989,298423,310631,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,300713,287224)
> 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] 62
> (np <- floor(n / par1))
[1] 5
> 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 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]
 [1,] 325412 337698 335083 303598 285257 300713
 [2,] 326011 335932 327616 286994 270976 287224
 [3,] 328282 323931 309119 276427 261076     NA
 [4,] 317480 313927 295916 266424 255603     NA
 [5,] 317539 314485 291413 267153 260376     NA
 [6,] 313737 313218 291542 268381 263903     NA
 [7,] 312276 309664 284678 262522 264291     NA
 [8,] 309391 302963 276475 255542 263276     NA
 [9,] 302950 298989 272566 253158 262572     NA
[10,] 300316 298423 264981 243803 256167     NA
[11,] 304035 310631 263290 250741 264221     NA
[12,] 333476 329765 296806 280445 293860     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/rcomp/tmp/140sn1257350807.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/rcomp/tmp/2kvi81257350807.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/rcomp/tmp/3nm8c1257350807.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/rcomp/tmp/4ov1a1257350807.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]
[1,] 285257 270976.0 261076 255603 260376 263903 262522 255542 253158 243803
[2,] 300713 286994.0 276427 266424 267153 268381 264291 263276 262572 256167
[3,] 314505 306617.5 309119 295916 291413 291542 284678 276475 272566 264981
[4,] 335083 327616.0 323931 313927 314485 313218 309664 302963 298989 298423
[5,] 337698 335932.0 328282 317480 317539 313737 312276 309391 302950 300316
      [,11]  [,12]
[1,] 250741 280445
[2,] 263290 293860
[3,] 264221 296806
[4,] 304035 329765
[5,] 310631 333476

$n
 [1] 6 6 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]   [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 292335.2 280415 275552.8 262350.5 257968.3 259860.3 252617.6 248432.3
[2,] 336674.8 332820 342685.2 329481.5 324857.7 323223.7 316738.4 304517.7
         [,9]  [,10]    [,11]    [,12]
[1,] 246833.8 235123 235430.7 271435.6
[2,] 298298.2 294839 293011.3 322176.4

$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(285257, 300713, 314505, 335083, 337698, 270976,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5jmmd1257350807.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]
[1,] 300316.0 298423.0 263290.0 243803.0 255603.0 287224.0
[2,] 306713.0 306313.5 274520.5 254350.0 260726.0 287224.0
[3,] 315608.5 313572.5 291477.5 266788.5 263589.5 293968.5
[4,] 325711.5 326848.0 302962.5 278436.0 267633.5 300713.0
[5,] 333476.0 337698.0 335083.0 303598.0 270976.0 300713.0

$n
[1] 12 12 12 12 12  2

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 306943.2 304206.6 278504.9 255802.7 260438.9 278898.2
[2,] 324273.8 322938.4 304450.1 277774.3 266740.1 309038.8

$out
[1] 285257 293860

$group
[1] 5 5

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

Warning message:
In bxp(list(stats = c(300316, 306713, 315608.5, 325711.5, 333476,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6t3571257350807.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,] 21189.01 49754.0 31231.0
[2,] 23595.40 51137.5 38052.0
[3,] 24645.75 55181.0 41500.5
[4,] 27048.60 60883.5 46352.5
[5,] 29690.96 67206.0 47504.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 23070.72 50735.78 37714.59
[2,] 26220.78 59626.22 45286.41

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(21189.0081827033, 23595.397940133, 24645.7529311954,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/140sn1257350807.ps tmp/140sn1257350807.png")
> system("convert tmp/2kvi81257350807.ps tmp/2kvi81257350807.png")
> system("convert tmp/3nm8c1257350807.ps tmp/3nm8c1257350807.png")
> system("convert tmp/4ov1a1257350807.ps tmp/4ov1a1257350807.png")
> system("convert tmp/5jmmd1257350807.ps tmp/5jmmd1257350807.png")
> system("convert tmp/6t3571257350807.ps tmp/6t3571257350807.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.390   1.260   2.424