R version 2.12.2 (2011-02-25)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-bit)

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.

  Natural language support but running in an English locale

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(2851,2672,2755,2721,2946,3036,2282,2212,2922,4301,5764,7132,2541,2475,3031,3266,3776,3230,3028,1759,3595,4474,6838,8357,3113,3006,4047,3523,3937,3986,3260,1573,3528,5211,7614,9254,5375,3088,3718,4514,4520,4539,3663,1643,4734,5428,8314,10651,3633,4292,4154,4121,4647,4753,3965,1723,5048,6923,9858,11331,4016,3957,4510,4276,4968,4677,3523,1821,5222,6872,10803,13916,2639,2899,3370,3740,2927,3986,4217,1738,5221,6424,9842,13076,3934,3162,4286,4676,5010,4874,4633,1659,5951,6981,9851,12670)
> 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] 96
> (np <- floor(n / par1))
[1] 8
> 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] 8 8 8 8 8 8 8 8 8 8 8 8
> arr
      [,1] [,2] [,3]  [,4]  [,5]  [,6]  [,7]  [,8] [,9]
 [1,] 2851 2541 3113  5375  3633  4016  2639  3934   NA
 [2,] 2672 2475 3006  3088  4292  3957  2899  3162   NA
 [3,] 2755 3031 4047  3718  4154  4510  3370  4286   NA
 [4,] 2721 3266 3523  4514  4121  4276  3740  4676   NA
 [5,] 2946 3776 3937  4520  4647  4968  2927  5010   NA
 [6,] 3036 3230 3986  4539  4753  4677  3986  4874   NA
 [7,] 2282 3028 3260  3663  3965  3523  4217  4633   NA
 [8,] 2212 1759 1573  1643  1723  1821  1738  1659   NA
 [9,] 2922 3595 3528  4734  5048  5222  5221  5951   NA
[10,] 4301 4474 5211  5428  6923  6872  6424  6981   NA
[11,] 5764 6838 7614  8314  9858 10803  9842  9851   NA
[12,] 7132 8357 9254 10651 11331 13916 13076 12670   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/wessaorg/rcomp/tmp/1navz1304972121.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2083j1304972121.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3llzt1304972121.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4pkcv1304972121.ps",horizontal=F,onefile=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]
[1,] 2541 2475.0 2755.0 2721.0 2927.0 3036.0 2282 1573.0 2922.0 4301.0  5764.0
[2,] 2745 2785.5 3200.5 3394.5 3361.0 3608.0 3144 1651.0 3561.5 4842.5  7226.0
[3,] 3373 3047.0 3882.5 3930.5 4228.5 4262.5 3593 1730.5 4891.0 5926.0  9078.0
[4,] 3975 3559.5 4220.0 4395.0 4807.5 4715.0 4091 1790.0 5221.5 6897.5  9854.5
[5,] 5375 4292.0 4510.0 4676.0 5010.0 4874.0 4633 1821.0 5951.0 6981.0 10803.0
       [,12]
[1,]  7132.0
[2,]  8805.5
[3,] 10991.0
[4,] 12873.0
[5,] 13916.0

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2685.904 2614.632 3312.993 3371.606 3420.464 3644.114 3063.992 1652.853
[2,] 4060.096 3479.368 4452.007 4489.394 5036.536 4880.886 4122.008 1808.147
       [,9]    [,10]     [,11]     [,12]
[1,] 3963.7 4778.047  7609.682  8718.836
[2,] 5818.3 7073.953 10546.318 13263.164

$out
[1] 2212

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(2541, 2745, 3373, 3975, 5375, 2475, 2785.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/5l02e1304972121.ps",horizontal=F,onefile=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] [,8] [,9]
[1,] 2212.0 1759.0 1573.0 1643.0 1723.0 1821.0 1738.0 1659   NA
[2,] 2696.5 2784.5 3186.5 3690.5 4043.0 3986.5 2913.0 4110   NA
[3,] 2886.5 3248.0 3732.5 4529.5 4469.5 4593.5 3863.0 4775   NA
[4,] 3668.5 4125.0 4629.0 5401.5 5985.5 6047.0 5822.5 6466   NA
[5,] 4301.0 4474.0 5211.0 5428.0 6923.0 6872.0 9842.0 9851   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2443.164 2636.589 3074.566 3749.101 3583.513 3653.692 2535.958 3700.413
[2,] 3329.836 3859.411 4390.434 5309.899 5355.487 5533.308 5190.042 5849.587
     [,9]
[1,]   NA
[2,]   NA

$out
 [1]  5764  7132  6838  8357  7614  9254  8314 10651  9858 11331 10803 13916
[13] 13076 12670

$group
 [1] 1 1 2 2 3 3 4 4 5 5 6 6 7 8

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

Warning message:
In bxp(list(stats = c(2212, 2696.5, 2886.5, 3668.5, 4301, 1759,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/6gtv31304972121.ps",horizontal=F,onefile=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,]  195.8316  639.0  119.500
[2,]  647.7544 1827.5  851.375
[3,]  784.8050 2217.0 1029.125
[4,] 1082.8373 2931.5 1750.500
[5,] 1113.5085 3029.0 2432.750

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 586.3608 1713.458  619.028
[2,] 983.2492 2720.542 1439.222

$out
[1] 1762.465 2405.448 5039.000 6784.000 3741.750

$group
[1] 1 1 2 2 3

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

Warning message:
In bxp(list(stats = c(195.831560275661, 647.754360654914, 784.804995998361,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1navz1304972121.ps tmp/1navz1304972121.png",intern=TRUE))
character(0)
> try(system("convert tmp/2083j1304972121.ps tmp/2083j1304972121.png",intern=TRUE))
character(0)
> try(system("convert tmp/3llzt1304972121.ps tmp/3llzt1304972121.png",intern=TRUE))
character(0)
> try(system("convert tmp/4pkcv1304972121.ps tmp/4pkcv1304972121.png",intern=TRUE))
character(0)
> try(system("convert tmp/5l02e1304972121.ps tmp/5l02e1304972121.png",intern=TRUE))
character(0)
> try(system("convert tmp/6gtv31304972121.ps tmp/6gtv31304972121.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.940   0.490   4.588