R version 2.8.0 (2008-10-20)
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(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,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] 60
> (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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5] [,6]
 [1,] 328282 323931 309119 276427 261076   NA
 [2,] 317480 313927 295916 266424 255603   NA
 [3,] 317539 314485 291413 267153 260376   NA
 [4,] 313737 313218 291542 268381 263903   NA
 [5,] 312276 309664 284678 262522 264291   NA
 [6,] 309391 302963 276475 255542 263276   NA
 [7,] 302950 298989 272566 253158 262572   NA
 [8,] 300316 298423 264981 243803 256167   NA
 [9,] 304035 301631 263290 250741 264221   NA
[10,] 333476 329765 296806 280445 293860   NA
[11,] 337698 335083 303598 285257 300713   NA
[12,] 335932 327616 286994 270976 287224   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/freestat/rcomp/tmp/1apn61257933301.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/freestat/rcomp/tmp/2psu51257933301.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/freestat/rcomp/tmp/3f0r21257933301.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/freestat/rcomp/tmp/48tcs1257933301.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,] 261076 255603 260376 263903 262522 255542 253158 243803 250741 280445
[2,] 276427 266424 267153 268381 264291 263276 262572 256167 263290 293860
[3,] 309119 295916 291413 291542 284678 276475 272566 264981 264221 296806
[4,] 323931 313927 314485 313218 309664 302963 298989 298423 301631 329765
[5,] 328282 317480 317539 313737 312276 309391 302950 300316 304035 333476
      [,11]  [,12]
[1,] 285257 270976
[2,] 300713 286994
[3,] 303598 287224
[4,] 335083 327616
[5,] 337698 335932

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]   [,8]
[1,] 275552.8 262350.5 257968.3 259860.3 252617.6 248432.3 246833.8 235123
[2,] 342685.2 329481.5 324857.7 323223.7 316738.4 304517.7 298298.2 294839
         [,9]    [,10]    [,11]    [,12]
[1,] 237129.3 271435.6 279312.2 258520.6
[2,] 291312.7 322176.4 327883.8 315927.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(261076, 276427, 309119, 323931, 328282, 255603,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/52wso1257933301.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   NA
[2,] 306713.0 302297.0 274520.5 254350.0 260726.0   NA
[3,] 315608.5 313572.5 289203.5 266788.5 263589.5   NA
[4,] 330879.0 325773.5 296361.0 273701.5 275757.5   NA
[5,] 337698.0 335083.0 309119.0 285257.0 293860.0   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 304586.2 302864.7 279241.9 257962.2 256733.5   NA
[2,] 326630.8 324280.3 299165.1 275614.8 270445.5   NA

$out
[1] 300713

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(300316, 306713, 315608.5, 330879, 337698,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6jxtu1257933301.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,] 22066.41 49754.0 34370.0
[2,] 23595.40 51137.5 37379.0
[3,] 24108.73 53571.5 41439.0
[4,] 27055.12 59520.0 46352.5
[5,] 29690.96 67206.0 47504.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 22530.73 49748.19 37346.13
[2,] 25686.73 57394.81 45531.87

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(22066.4115795931, 23595.397940133, 24108.7290444142,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1apn61257933301.ps tmp/1apn61257933301.png")
> system("convert tmp/2psu51257933301.ps tmp/2psu51257933301.png")
> system("convert tmp/3f0r21257933301.ps tmp/3f0r21257933301.png")
> system("convert tmp/48tcs1257933301.ps tmp/48tcs1257933301.png")
> system("convert tmp/52wso1257933301.ps tmp/52wso1257933301.png")
> system("convert tmp/6jxtu1257933301.ps tmp/6jxtu1257933301.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.639   1.265   2.042