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(2849.27,2921.44,2981.85,3080.58,3106.22,3119.31,3061.26,3097.31,3161.69,3257.16,3277.01,3295.32,3363.99,3494.17,3667.03,3813.06,3917.96,3895.51,3801.06,3570.12,3701.61,3862.27,3970.1,4138.52,4199.75,4290.89,4443.91,4502.64,4356.98,4591.27,4696.96,4621.4,4562.84,4202.52,4296.49,4435.23,4105.18,4116.68,3844.49,3720.98,3674.4,3857.62,3801.06,3504.37,3032.6,3047.03,2962.34,2197.82,2014.45,1862.83,1905.41,1810.99,1670.07,1864.44,2052.02,2029.6,2070.83,2293.41,2443.27,2513.17)
> 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,] 2849.27 3363.99 4199.75 4105.18 2014.45   NA
 [2,] 2921.44 3494.17 4290.89 4116.68 1862.83   NA
 [3,] 2981.85 3667.03 4443.91 3844.49 1905.41   NA
 [4,] 3080.58 3813.06 4502.64 3720.98 1810.99   NA
 [5,] 3106.22 3917.96 4356.98 3674.40 1670.07   NA
 [6,] 3119.31 3895.51 4591.27 3857.62 1864.44   NA
 [7,] 3061.26 3801.06 4696.96 3801.06 2052.02   NA
 [8,] 3097.31 3570.12 4621.40 3504.37 2029.60   NA
 [9,] 3161.69 3701.61 4562.84 3032.60 2070.83   NA
[10,] 3257.16 3862.27 4202.52 3047.03 2293.41   NA
[11,] 3277.01 3970.10 4296.49 2962.34 2443.27   NA
[12,] 3295.32 4138.52 4435.23 2197.82 2513.17   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/1c9l31259511233.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/271fk1259511233.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/366r31259511233.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/4itb71259511233.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,] 2014.45 1862.83 1905.41 3080.58 3106.22 3119.31 2052.02 3097.31 2070.83
[2,] 2849.27 2921.44 2981.85 3080.58 3106.22 3119.31 3061.26 3097.31 3032.60
[3,] 3363.99 3494.17 3667.03 3720.98 3674.40 3857.62 3801.06 3504.37 3161.69
[4,] 4105.18 4116.68 3844.49 3813.06 3917.96 3895.51 3801.06 3570.12 3701.61
[5,] 4199.75 4290.89 4443.91 4502.64 4356.98 4591.27 4696.96 3570.12 4562.84
       [,10]   [,11]   [,12]
[1,] 2293.41 2443.27 2197.82
[2,] 3047.03 2962.34 2513.17
[3,] 3257.16 3277.01 3295.32
[4,] 3862.27 3970.10 4138.52
[5,] 4202.52 4296.49 4435.23

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2476.567 2649.616 3057.491 3203.411 3100.827 3309.159 3278.319 3170.284
[2,] 4251.413 4338.724 4276.569 4238.549 4247.973 4406.081 4323.801 3838.456
         [,9]    [,10]    [,11]    [,12]
[1,] 2688.969 2681.113 2564.929 2146.852
[2,] 3634.411 3833.207 3989.091 4443.788

$out
[1] 1810.99 1670.07 1864.44 4621.40 2029.60

$group
[1] 4 5 6 8 8

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

Warning message:
In bxp(list(stats = c(2014.45, 2849.27, 3363.99, 4105.18, 4199.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5mjr51259511233.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,] 2849.270 3363.990 4199.750 2197.820 1670.070   NA
[2,] 3021.555 3618.575 4293.690 3039.815 1863.635   NA
[3,] 3101.765 3807.060 4439.570 3697.690 2022.025   NA
[4,] 3209.425 3906.735 4577.055 3851.055 2182.120   NA
[5,] 3295.320 4138.520 4696.960 4116.680 2513.170   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 3016.076 3675.628 4310.325 3327.678 1876.762   NA
[2,] 3187.454 3938.492 4568.815 4067.702 2167.288   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2849.27, 3021.555, 3101.765, 3209.425, 3295.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/63knm1259511233.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,]  910.7062 1853.220  472.81
[2,]  913.4314 2211.355  736.14
[3,]  973.8263 2515.255  813.49
[4,] 1000.9732 2665.925 1101.50
[5,] 1039.6270 2726.830 1625.35

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,]  933.898 2307.923 646.8469
[2,] 1013.755 2722.587 980.1331

$out
[1] 742.3726 750.0329

$group
[1] 1 1

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

Warning message:
In bxp(list(stats = c(910.706165082899, 913.431397742573, 973.826318449808,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1c9l31259511233.ps tmp/1c9l31259511233.png")
> system("convert tmp/271fk1259511233.ps tmp/271fk1259511233.png")
> system("convert tmp/366r31259511233.ps tmp/366r31259511233.png")
> system("convert tmp/4itb71259511233.ps tmp/4itb71259511233.png")
> system("convert tmp/5mjr51259511233.ps tmp/5mjr51259511233.png")
> system("convert tmp/63knm1259511233.ps tmp/63knm1259511233.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.025   0.799   1.281