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(6430,5124,4836,4629,4597,4490,4517,4560,4135,4559,4739,4886,5605,4616,4997,4607,4882,4555,4462,4476,4277,4369,4492,5183,6039,4923,4953,4892,4614,4363,4675,4556,4217,4664,4601,5428,5607,4869,5174,5031,4671,4491,4504,4615,4582,4800,4775,5791,5818,4714,4915,4598,4407,4383,4412,4274,4236,4637,4534,5271,5467,5204,5752,4724,4623,4451,4138,4140,4169,4603,4434,5185)
> 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,] 6430 5605 6039 5607 5818 5467   NA
 [2,] 5124 4616 4923 4869 4714 5204   NA
 [3,] 4836 4997 4953 5174 4915 5752   NA
 [4,] 4629 4607 4892 5031 4598 4724   NA
 [5,] 4597 4882 4614 4671 4407 4623   NA
 [6,] 4490 4555 4363 4491 4383 4451   NA
 [7,] 4517 4462 4675 4504 4412 4138   NA
 [8,] 4560 4476 4556 4615 4274 4140   NA
 [9,] 4135 4277 4217 4582 4236 4169   NA
[10,] 4559 4369 4664 4800 4637 4603   NA
[11,] 4739 4492 4601 4775 4534 4434   NA
[12,] 4886 5183 5428 5791 5271 5185   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/1o8c31250538465.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/21x8f1250538465.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/3y5cg1250538465.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/457uh1250538465.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,] 5467.0 4616 4836 4598.0 4597.0 4363.0 4412 4140 4135.0  4559 4434.0  4886
[2,] 5605.0 4714 4915 4607.0 4597.0 4383.0 4412 4274 4169.0  4559 4492.0  5183
[3,] 5712.5 4896 4975 4676.5 4618.5 4470.5 4483 4516 4226.5  4620 4567.5  5228
[4,] 6039.0 5124 5174 4892.0 4671.0 4491.0 4517 4560 4277.0  4664 4739.0  5428
[5,] 6430.0 5204 5174 5031.0 4671.0 4555.0 4517 4615 4277.0  4800 4775.0  5791

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 5432.556 4631.537 4807.937 4492.666 4570.768 4400.837 4415.272 4331.521
[2,] 5992.444 5160.463 5142.063 4860.334 4666.232 4540.163 4550.728 4700.479
         [,9]    [,10]    [,11]    [,12]
[1,] 4156.837 4552.272 4408.177 5069.967
[2,] 4296.163 4687.728 4726.823 5386.033

$out
[1] 5752 4882 4407 4675 4138 4582 4369

$group
[1]  3  5  5  7  7  9 10

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

Warning message:
In bxp(list(stats = c(5467, 5605, 5712.5, 6039, 6430, 4616, 4714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/58u3e1250538465.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,] 4135 4277.0 4217.0 4491.0 4236.0 4138.0   NA
[2,] 4538 4469.0 4578.5 4598.5 4395.0 4301.5   NA
[3,] 4613 4581.0 4669.5 4787.5 4566.0 4613.0   NA
[4,] 4861 4939.5 4938.0 5102.5 4814.5 5194.5   NA
[5,] 5124 5605.0 5428.0 5791.0 5271.0 5752.0   NA

$n
[1] 12 12 12 12 12 12  0

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6] [,7]
[1,] 4465.678 4366.402 4505.53 4557.622 4374.663 4205.697   NA
[2,] 4760.322 4795.598 4833.47 5017.378 4757.337 5020.303   NA

$out
[1] 6430 6039 5818

$group
[1] 1 3 5

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

Warning message:
In bxp(list(stats = c(4135, 4538, 4613, 4861, 5124, 4277, 4469,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6x07i1250538465.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,]  72.36505 192.0  57.750
[2,] 147.08073 432.0  88.250
[3,] 177.45011 475.0 203.625
[4,] 264.79205 746.5 236.000
[5,] 356.62735 963.0 378.250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 123.7612 331.5544 136.2352
[2,] 231.1390 618.4456 271.0148

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(72.3650468112887, 147.080726748142, 177.450105300303,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1o8c31250538465.ps tmp/1o8c31250538465.png")
> system("convert tmp/21x8f1250538465.ps tmp/21x8f1250538465.png")
> system("convert tmp/3y5cg1250538465.ps tmp/3y5cg1250538465.png")
> system("convert tmp/457uh1250538465.ps tmp/457uh1250538465.png")
> system("convert tmp/58u3e1250538465.ps tmp/58u3e1250538465.png")
> system("convert tmp/6x07i1250538465.ps tmp/6x07i1250538465.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.068   0.814   1.743