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(15136
+ ,16733
+ ,20016
+ ,17708
+ ,18019
+ ,19227
+ ,22893
+ ,23739
+ ,21133
+ ,22591
+ ,26786
+ ,29740
+ ,15028
+ ,17977
+ ,20008
+ ,21354
+ ,19498
+ ,22125
+ ,25817
+ ,28779
+ ,20960
+ ,22254
+ ,27392
+ ,29945
+ ,16933
+ ,17892
+ ,20533
+ ,23569
+ ,22417
+ ,22084
+ ,26580
+ ,27454
+ ,24081
+ ,23451
+ ,28991
+ ,31386
+ ,16896
+ ,20045
+ ,23471
+ ,21747
+ ,25621
+ ,23859
+ ,25500
+ ,30998
+ ,24475
+ ,23145
+ ,29701
+ ,34365
+ ,17556
+ ,22077
+ ,25702
+ ,22214
+ ,26886
+ ,23191
+ ,27831
+ ,35406
+ ,23195
+ ,25110
+ ,30009
+ ,36242
+ ,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
+ ,13652
+ ,22784
+ ,23565
+ ,26323
+ ,23779
+ ,27549
+ ,29660
+ ,23356)
> 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] 176
> (np <- floor(n / par1))
[1] 14
> 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] 15 15 15 15 15 15 15 15 14 14 14 14
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 15136 15028 16933 16896 17556 18450 19969 21337 20504 16991 14672 16052
 [2,] 16733 17977 17892 20045 22077 21845 21701 19419 23060 21109 17543 22146
 [3,] 20016 20008 20533 23471 25702 26488 26249 23166 23562 23740 25453 21198
 [4,] 17708 21354 23569 21747 22214 22394 24493 28286 27562 25552 32683 19543
 [5,] 18019 19498 22417 25621 26886 28057 24603 24570 23940 21752 22449 22084
 [6,] 19227 22125 22084 23859 23191 25451 26485 24001 24584 20294 22316 23816
 [7,] 22893 25817 26580 25500 27831 24872 30723 33151 34303 29009 27595 29961
 [8,] 23739 28779 27454 30998 35406 33424 34569 24878 25517 25500 25451 26773
 [9,] 21133 20960 24081 24475 23195 24052 26689 26804 23494 24166 25421 26635
[10,] 22591 22254 23451 23145 25110 28449 26157 28967 29095 26960 25288 26972
[11,] 26786 27392 28991 29701 30009 33533 32064 33311 32903 31222 32568 30207
[12,] 29740 29945 31386 34365 36242 37351 38870 40226 34379 38641 35110 38687
      [,13] [,14] [,15]
 [1,] 16974 17466 13652
 [2,] 21697 19463 22784
 [3,] 24179 24352 23565
 [4,] 23757 26805 26323
 [5,] 25013 25236 23779
 [6,] 24019 24735 27549
 [7,] 30345 29356 29660
 [8,] 24488 31234 23356
 [9,] 25156 22724    NA
[10,] 25650 28496    NA
[11,] 30923 32857    NA
[12,] 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/1rhnh1273507754.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/2rhnh1273507754.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/3k84k1273507754.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/4k84k1273507754.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,] 13652 16733 20008.0 17708.0 18019.0 19227.0 22893.0 23356.0 20960.0
[2,] 15594 18698 22182.0 21980.5 22250.5 22220.5 26198.5 25164.5 23195.0
[3,] 16974 21109 23565.0 23757.0 23940.0 23859.0 29009.0 26773.0 24123.5
[4,] 18003 21961 24902.5 26564.0 25124.5 24659.5 30153.0 31116.0 25421.0
[5,] 21337 23060 26488.0 32683.0 28057.0 27549.0 34303.0 35406.0 26804.0
       [,10]   [,11] [,12]
[1,] 22254.0 26786.0 29740
[2,] 23451.0 29701.0 34365
[3,] 25903.5 31072.5 36720
[4,] 28449.0 32857.0 38641
[5,] 29095.0 33533.0 40226

$n
 [1] 15 15 15 15 15 15 15 15 14 14 14 14

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]  [,6]     [,7]     [,8]
[1,] 15991.24 19777.85 22455.16 21887.14 22767.54 22864 27395.74 24345.06
[2,] 17956.76 22440.15 24674.84 25626.86 25112.46 24854 30622.26 29200.94
         [,9]    [,10]    [,11]    [,12]
[1,] 23183.52 23792.98 29739.81 34914.36
[2,] 25063.48 28014.02 32405.19 38525.64

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(13652, 15594, 16974, 18003, 21337, 16733,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5k84k1273507754.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,] 15136.0 15028.0 16933.0 16896 17556.0 18450.0 19969.0 19419.0 20504.0
[2,] 17863.5 19753.0 21308.5 22446 22702.5 23223.0 24548.0 23583.5 23528.0
[3,] 20574.5 21739.5 23510.0 24167 25406.0 25969.5 26367.0 25841.0 25050.5
[4,] 23316.0 26604.5 27017.0 27661 28920.0 30936.5 31393.5 31059.0 30999.0
[5,] 29740.0 29945.0 31386.0 34365 36242.0 37351.0 38870.0 40226.0 34379.0
       [,10]   [,11]   [,12]   [,13]   [,14] [,15]
[1,] 16991.0 14672.0 16052.0 21697.0 17466.0 22784
[2,] 21430.5 22382.5 21641.0 23888.0 23538.0 23070
[3,] 24833.0 25436.0 25225.5 24750.5 26020.5 23672
[4,] 27984.5 30081.5 28466.5 27997.5 30295.0 26936
[5,] 31222.0 35110.0 38687.0 30923.0 37198.0 29660

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 18087.58 18614.48 20906.31 21788.40 22570.16 22451.32 23244.72 22431.37
[2,] 23061.42 24864.52 26113.69 26545.60 28241.84 29487.68 29489.28 29250.63
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]
[1,] 21642.93 21843.68 21924.43 22112.34 22876.13 22938.59 21512.40
[2,] 28458.07 27822.32 28947.57 28338.66 26624.87 29102.41 25831.60

$out
[1] 38641 16974 37240 13652

$group
[1] 10 13 13 15

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

Warning message:
In bxp(list(stats = c(15136, 17863.5, 20574.5, 23316, 29740, 15028,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6v0mn1273507754.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,] 1860.758  5844.0 2085.000
[2,] 2154.439  6613.5 2579.750
[3,] 2287.250  8003.5 3134.875
[4,] 3226.689 10948.0 4084.250
[5,] 4051.725 14975.0 5951.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1798.190 6026.505 2448.662
[2,] 2776.311 9980.495 3821.088

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1860.75804150109, 2154.43904226937, 2287.2502746719,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1rhnh1273507754.ps tmp/1rhnh1273507754.png",intern=TRUE))
character(0)
> try(system("convert tmp/2rhnh1273507754.ps tmp/2rhnh1273507754.png",intern=TRUE))
character(0)
> try(system("convert tmp/3k84k1273507754.ps tmp/3k84k1273507754.png",intern=TRUE))
character(0)
> try(system("convert tmp/4k84k1273507754.ps tmp/4k84k1273507754.png",intern=TRUE))
character(0)
> try(system("convert tmp/5k84k1273507754.ps tmp/5k84k1273507754.png",intern=TRUE))
character(0)
> try(system("convert tmp/6v0mn1273507754.ps tmp/6v0mn1273507754.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.077   0.844   1.381