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(18450,21845,26488,22394,28057,25451,24872,33424,24052,28449,33533,37351,19969,21701,26249,24493,24603,26485,30723,34569,26689,26157,32064,38870,21337,19419,23166,28286,24570,24001,33151,24878,26804,28967,33311,40226,20504,23060,23562,27562,23940,24584,34303,25517,23494,29095,32903,34379,16991,21109,23740,25552,21752,20294,29009,25500,24166,26960,31222,38641,14672,17543,25453,32683,22449,22316,27595,25451,25421,25288,32568,35110,16052,22146,21198,19543,22084,23816,29961,26773,26635,26972,30207,38687,16974,21697,24179,23757,25013,24019,30345,24488,25156,25650,30923,37240,17466,19463,24352,26805,25236,24735,29356,31234,22724,28496,32857,37198)
> 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] 108
> (np <- floor(n / par1))
[1] 9
> 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] 9 9 9 9 9 9 9 9 9 9 9 9
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10]
 [1,] 18450 19969 21337 20504 16991 14672 16052 16974 17466    NA
 [2,] 21845 21701 19419 23060 21109 17543 22146 21697 19463    NA
 [3,] 26488 26249 23166 23562 23740 25453 21198 24179 24352    NA
 [4,] 22394 24493 28286 27562 25552 32683 19543 23757 26805    NA
 [5,] 28057 24603 24570 23940 21752 22449 22084 25013 25236    NA
 [6,] 25451 26485 24001 24584 20294 22316 23816 24019 24735    NA
 [7,] 24872 30723 33151 34303 29009 27595 29961 30345 29356    NA
 [8,] 33424 34569 24878 25517 25500 25451 26773 24488 31234    NA
 [9,] 24052 26689 26804 23494 24166 25421 26635 25156 22724    NA
[10,] 28449 26157 28967 29095 26960 25288 26972 25650 28496    NA
[11,] 33533 32064 33311 32903 31222 32568 30207 30923 32857    NA
[12,] 37351 38870 40226 34379 38641 35110 38687 37240 37198    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/1b9wy1273854246.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/2b9wy1273854246.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/3b9wy1273854246.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/43iej1273854246.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,] 14672 17543 21198 19543 21752 23816 27595 24488 22724 25288 30207 35110
[2,] 16974 19463 23562 23757 22449 23816 29009 25451 24052 26157 31222 37198
[3,] 17466 21697 24179 25552 24570 24019 29961 25517 25156 26972 32568 37351
[4,] 19969 21845 25453 27562 25013 24735 30723 31234 26635 28496 32903 38687
[5,] 21337 23060 26488 32683 28057 25451 33151 34569 26804 29095 33533 40226

$n
 [1] 9 9 9 9 9 9 9 9 9 9 9 9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 15888.63 20442.48 23183.07 23548.03 23219.63 23534.99 29058.29 22471.29
[2,] 19043.37 22951.52 25174.93 27555.97 25920.37 24503.01 30863.71 28562.71
         [,9]    [,10]    [,11]    [,12]
[1,] 23795.62 25740.13 31682.67 36566.79
[2,] 26516.38 28203.87 33453.33 38135.21

$out
[1] 26485 20294 22316 24872 34303 34379

$group
[1]  6  6  6  7  7 12

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

Warning message:
In bxp(list(stats = c(14672, 16974, 17466, 19969, 21337, 17543,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/53iej1273854246.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]    [,8]    [,9]
[1,] 18450.0 19969.0 19419.0 20504.0 16991.0 14672.0 16052.0 21697.0 17466.0
[2,] 23223.0 24548.0 23583.5 23528.0 21430.5 22382.5 21641.0 23888.0 23538.0
[3,] 25969.5 26367.0 25841.0 25050.5 24833.0 25436.0 25225.5 24750.5 26020.5
[4,] 30936.5 31393.5 31059.0 30999.0 27984.5 30081.5 28466.5 27997.5 30295.0
[5,] 37351.0 38870.0 40226.0 34379.0 31222.0 35110.0 38687.0 30923.0 37198.0
     [,10]
[1,]    NA
[2,]    NA
[3,]    NA
[4,]    NA
[5,]    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 22451.32 23244.72 22431.37 21642.93 21843.68 21924.43 22112.34 22876.13
[2,] 29487.68 29489.28 29250.63 28458.07 27822.32 28947.57 28338.66 26624.87
         [,9] [,10]
[1,] 22938.59    NA
[2,] 29102.41    NA

$out
[1] 38641 16974 37240

$group
[1] 5 8 8

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

Warning message:
In bxp(list(stats = c(18450, 23223, 25969.5, 30936.5, 37351, 19969,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wsvm1273854246.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,] 1153.541  3326  919.0
[2,] 1573.726  4685 1697.5
[3,] 1826.673  6019 2360.5
[4,] 2497.291  8048 2789.0
[5,] 3778.631 10081 3805.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1405.429 4485.113 1862.660
[2,] 2247.918 7552.887 2858.340

$out
[1]  3961.163 13140.000  5783.000

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(1153.54108196361, 1573.72605592677, 1826.67332707107,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1b9wy1273854246.ps tmp/1b9wy1273854246.png",intern=TRUE))
character(0)
> try(system("convert tmp/2b9wy1273854246.ps tmp/2b9wy1273854246.png",intern=TRUE))
character(0)
> try(system("convert tmp/3b9wy1273854246.ps tmp/3b9wy1273854246.png",intern=TRUE))
character(0)
> try(system("convert tmp/43iej1273854246.ps tmp/43iej1273854246.png",intern=TRUE))
character(0)
> try(system("convert tmp/53iej1273854246.ps tmp/53iej1273854246.png",intern=TRUE))
character(0)
> try(system("convert tmp/6wsvm1273854246.ps tmp/6wsvm1273854246.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.046   0.807   1.330