R version 2.6.2 (2008-02-08)
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.

  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(2815,2672,2755,2721,2946,3036,2282,2212,2922,4301,5764,7312,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,4739,5428,8314,10651,3633,4292,4154,4121,4647,4753,3965,1723,5048,6922,9858,11331,4016,4276,4968,4677,3523,1821,5222,6872,10803,13916,2639,2899,3370,3740,2927,3986,4217,1738,5221,6424,9842,13076)
> 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] 82
> (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] 7 7 7 7 7 7 7 7 7 7 6 6
> arr
      [,1] [,2] [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 2815 2541 3113  5375  3633  4016  3370
 [2,] 2672 2475 3006  3088  4292  4276  3740
 [3,] 2755 3031 4047  3718  4154  4968  2927
 [4,] 2721 3266 3523  4514  4121  4677  3986
 [5,] 2946 3776 3937  4520  4647  3523  4217
 [6,] 3036 3230 3986  4539  4753  1821  1738
 [7,] 2282 3028 3260  3663  3965  5222  5221
 [8,] 2212 1759 1573  1643  1723  6872  6424
 [9,] 2922 3595 3528  4739  5048 10803  9842
[10,] 4301 4474 5211  5428  6922 13916 13076
[11,] 5764 6838 7614  8314  9858  2639    NA
[12,] 7312 8357 9254 10651 11331  2899    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/1ml2w1211126337.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/2ktu31211126337.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/3cn641211126337.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/4ta1j1211126337.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]
[1,] 2541.0 2475 2755.0 2721.0 2946.0 1738.0 2282 1573  2922.0  4301.0  2639
[2,] 2964.0 2839 2979.0 3394.5 3649.5 2428.5 3144 1683  3561.5  4842.5  5764
[3,] 3370.0 3088 3718.0 3986.0 3937.0 3230.0 3663 1759  4739.0  5428.0  7226
[4,] 3824.5 4008 4100.5 4317.5 4368.5 4262.5 4593 4318  7445.0  9999.0  8314
[5,] 4016.0 4292 4968.0 4677.0 4647.0 4753.0 5222 6872 10803.0 13916.0  9858
       [,12]
[1,]  2899.0
[2,]  7312.0
[3,]  8805.5
[4,] 10651.0
[5,] 11331.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]      [,8]
[1,] 2856.123 2389.892 3048.258 3434.799 3507.625 2134.765 2797.681  185.4205
[2,] 3883.877 3786.108 4387.742 4537.201 4366.375 4325.235 4528.319 3332.5795
         [,9]    [,10]    [,11]     [,12]
[1,] 2419.836 2348.621 5581.168  6651.737
[2,] 7058.164 8507.379 8870.832 10959.263

$out
[1] 5375

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(2541, 2964, 3370, 3824.5, 4016, 2475, 2839,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/50vzs1211126337.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,] 2212.0 1759.0 1573.0 1643.0 1723.0 1821.0 1738.0
[2,] 2696.5 2784.5 3186.5 3690.5 4043.0 3211.0 3370.0
[3,] 2868.5 3248.0 3732.5 4529.5 4469.5 4476.5 4101.5
[4,] 3668.5 4125.0 4629.0 5401.5 5985.0 6047.0 6424.0
[5,] 4301.0 4474.0 5211.0 5428.0 6922.0 6872.0 9842.0

$n
[1] 12 12 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]   [,7]
[1,] 2425.164 2636.589 3074.566 3749.101 3583.741 3182.981 2575.6
[2,] 3311.836 3859.411 4390.434 5309.899 5355.259 5770.019 5627.4

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

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

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

Warning message:
In bxp(list(stats = c(2212, 2696.5, 2868.5, 3668.5, 4301, 1759,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6m7s01211126337.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,]  591.8587 1701.0  719.00
[2,]  772.2862 2084.5 1022.25
[3,] 1156.8955 2977.5 1641.50
[4,] 2751.5260 7550.0 2681.75
[5,] 4111.2844 9615.0 5156.50

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,]  254.1509  484.6488  884.5909
[2,] 2059.6400 5470.3512 2398.4091

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(591.858654297347, 772.286236976461, 1156.89546785552,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ml2w1211126337.ps tmp/1ml2w1211126337.png")
> system("convert tmp/2ktu31211126337.ps tmp/2ktu31211126337.png")
> system("convert tmp/3cn641211126337.ps tmp/3cn641211126337.png")
> system("convert tmp/4ta1j1211126337.ps tmp/4ta1j1211126337.png")
> system("convert tmp/50vzs1211126337.ps tmp/50vzs1211126337.png")
> system("convert tmp/6m7s01211126337.ps tmp/6m7s01211126337.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.388   1.322   2.600