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(3592.21,5955.74,4652.25,4211.65,4787.85,3599.73,4174.27,5106.33,5325.75,6604.61,5711.90,6919.30,7048.76,8655.98,6658.53,7247.03,8779.57,6602.49,9832.48,9369.49,8582.76,8206.94,6515.83,8618.10,8505.39,9881.64,9375.29,15642.50,12232.73,6288.93,12473.94,11142.82,10236.32,10581.51,8763.71,10819.04,11636.25,14650.13,10671.38,17468.63,13873.19,13077.58,16866.81,14186.64,19919.87,17681.78,9984.28,17423.09,13514.45,12334.57,12274.56,11752.23,13054.00,12460.98,8626.68,13722.62,12066.22,6798.83,6593.82,6606.19,6315.28,7232.95,6747.44,7803.61,6700.31,5369.53,8081.19,10718.39,9447.21,6815.10,5497.80,6805.31)
> 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,] 3592.21 7048.76  8505.39 11636.25 13514.45  6315.28   NA
 [2,] 5955.74 8655.98  9881.64 14650.13 12334.57  7232.95   NA
 [3,] 4652.25 6658.53  9375.29 10671.38 12274.56  6747.44   NA
 [4,] 4211.65 7247.03 15642.50 17468.63 11752.23  7803.61   NA
 [5,] 4787.85 8779.57 12232.73 13873.19 13054.00  6700.31   NA
 [6,] 3599.73 6602.49  6288.93 13077.58 12460.98  5369.53   NA
 [7,] 4174.27 9832.48 12473.94 16866.81  8626.68  8081.19   NA
 [8,] 5106.33 9369.49 11142.82 14186.64 13722.62 10718.39   NA
 [9,] 5325.75 8582.76 10236.32 19919.87 12066.22  9447.21   NA
[10,] 6604.61 8206.94 10581.51 17681.78  6798.83  6815.10   NA
[11,] 5711.90 6515.83  8763.71  9984.28  6593.82  5497.80   NA
[12,] 6919.30 8618.10 10819.04 17423.09  6606.19  6805.31   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/18jll1273933239.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/2is361273933239.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/3is361273933239.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/4is361273933239.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]
[1,]  3592.210  5955.74  4652.250  4211.65  4787.85  3599.73  4174.27  5106.33
[2,]  6315.280  7232.95  6658.530  7247.03  6700.31  5369.53  8081.19  9369.49
[3,]  7777.075  9268.81  8061.365  9777.92 10506.15  6445.71  9229.58 10930.60
[4,] 11636.250 12334.57 10671.380 15642.50 13054.00 12460.98 12473.94 13722.62
[5,] 13514.450 14650.13 12274.560 17468.63 13873.19 13077.58 16866.81 14186.64
          [,9]    [,10]    [,11]    [,12]
[1,]  5325.750  6604.61 5497.800  6606.19
[2,]  8582.760  6798.83 5711.900  6805.31
[3,]  9841.765  7511.02 6554.825  7768.70
[4,] 12066.220 10581.51 8763.710 10819.04
[5,] 12066.220 10581.51 9984.280 10819.04

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

$conf
          [,1]     [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  4344.877  5978.10  5472.947  4362.571  6407.815  1871.496  6396.114
[2,] 11209.273 12559.52 10649.783 15193.269 14604.485 11019.924 12063.046
          [,8]      [,9]    [,10]    [,11]     [,12]
[1,]  8122.696  7594.821 5071.069 4586.309  5179.715
[2,] 13738.514 12088.709 9950.971 8523.341 10357.685

$out
[1] 19919.87 17681.78 17423.09

$group
[1]  9 10 12

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

Warning message:
In bxp(list(stats = c(3592.21, 6315.28, 7777.075, 11636.25, 13514.45,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5bj2r1273933239.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,] 3592.21 6515.830  6288.93  9984.28  6593.820 5369.530   NA
[2,] 4192.96 6853.645  9069.50 12356.92  7712.755 6507.795   NA
[3,] 4947.09 8394.850 10408.92 14418.38 12170.390 6810.205   NA
[4,] 5833.82 8717.775 11687.77 17445.86 12757.490 7942.400   NA
[5,] 6919.30 9832.480 12473.94 19919.87 13722.620 9447.210   NA

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

$conf
         [,1]     [,2]      [,3]     [,4]      [,5]     [,6] [,7]
[1,] 4198.683 7544.608  9214.702 12097.28  9869.453 6155.872   NA
[2,] 5695.497 9245.092 11603.128 16739.49 14471.327 7464.538   NA

$out
[1] 15642.50 10718.39

$group
[1] 3 6

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

Warning message:
In bxp(list(stats = c(3592.21, 4192.96, 4947.09, 5833.82, 6919.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6bj2r1273933239.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,] 1796.772  4486.480 2308.355
[2,] 3273.680  8887.350 3277.963
[3,] 3814.158  9700.045 3631.306
[4,] 4306.545 11884.855 4875.931
[5,] 5172.228 14594.120 5628.557

$n
[1] 12 12 12

$conf
         [,1]      [,2]     [,3]
[1,] 3343.061  8332.863 2902.462
[2,] 4285.254 11067.227 4360.151

$out
[1] 7283.757

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(1796.77156039381, 3273.67974588052, 3814.15767992764,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/18jll1273933239.ps tmp/18jll1273933239.png",intern=TRUE))
character(0)
> try(system("convert tmp/2is361273933239.ps tmp/2is361273933239.png",intern=TRUE))
character(0)
> try(system("convert tmp/3is361273933239.ps tmp/3is361273933239.png",intern=TRUE))
character(0)
> try(system("convert tmp/4is361273933239.ps tmp/4is361273933239.png",intern=TRUE))
character(0)
> try(system("convert tmp/5bj2r1273933239.ps tmp/5bj2r1273933239.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bj2r1273933239.ps tmp/6bj2r1273933239.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.090   0.845  33.762