R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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.

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(192528,190773,188996,185320,221698,219771,192528,174414,176163,176163,178112,181616,192528,188996,194449,203412,254400,254400,243516,232604,241567,252473,254400,259853,276217,265306,265306,281670,327033,330709,321580,299762,316121,316121,317876,327033,334241,337917,337917,348823,390681,401565,403314,376071,390681,385228,374317,397889,403314,394185,396112,408773,456085,479624,479624,468740,485082,468740,459588,494239,499664,486832,519533,532366,570521,595842,592339,590389,604999,603222,581432,614128,625040,614128,659491,681309,732102,752143,746712,735800,744935,755841,719441,748461,766753,759373,806657,822993,892101,904761,888424,897554,903007,908460,873809,906511,924624,906511,959275,975618,1046468,1057380,1060884,1079170,1079170,1086378,1053676,1070041,1080925,1060884,1119079,1129986,1202597,1215430,1233543,1249908,1251657,1253584,1220883,1253584)
> 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] 120
> (np <- floor(n / par1))
[1] 10
> 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] 10 10 10 10 10 10 10 10 10 10 10 10
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]    [,9]   [,10]
 [1,] 192528 192528 276217 334241 403314 499664 625040 766753  924624 1080925
 [2,] 190773 188996 265306 337917 394185 486832 614128 759373  906511 1060884
 [3,] 188996 194449 265306 337917 396112 519533 659491 806657  959275 1119079
 [4,] 185320 203412 281670 348823 408773 532366 681309 822993  975618 1129986
 [5,] 221698 254400 327033 390681 456085 570521 732102 892101 1046468 1202597
 [6,] 219771 254400 330709 401565 479624 595842 752143 904761 1057380 1215430
 [7,] 192528 243516 321580 403314 479624 592339 746712 888424 1060884 1233543
 [8,] 174414 232604 299762 376071 468740 590389 735800 897554 1079170 1249908
 [9,] 176163 241567 316121 390681 485082 604999 744935 903007 1079170 1251657
[10,] 176163 252473 316121 385228 468740 603222 755841 908460 1086378 1253584
[11,] 178112 254400 317876 374317 459588 581432 719441 873809 1053676 1220883
[12,] 181616 259853 327033 397889 494239 614128 748461 906511 1070041 1253584
      [,11]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    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/wessaorg/rcomp/tmp/1yiej1312382321.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/wessaorg/rcomp/tmp/23b1n1312382321.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/wessaorg/rcomp/tmp/3o4vz1312382321.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/wessaorg/rcomp/tmp/4s70e1312382321.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]
[1,]  192528  188996.0  188996.0  185320.0  221698  219771  192528.0  174414.0
[2,]  276217  265306.0  265306.0  281670.0  327033  330709  321580.0  299762.0
[3,]  451489  440508.5  457822.5  470569.5  513303  537733  535981.5  529564.5
[4,]  766753  759373.0  806657.0  822993.0  892101  904761  888424.0  897554.0
[5,] 1080925 1060884.0 1119079.0 1129986.0 1202597 1215430 1233543.0 1249908.0
          [,9]   [,10]   [,11]     [,12]
[1,]  176163.0  176163  178112  181616.0
[2,]  316121.0  316121  317876  327033.0
[3,]  545040.5  535981  520510  554183.5
[4,]  903007.0  908460  873809  906511.0
[5,] 1251657.0 1253584 1220883 1253584.0

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

$conf
         [,1]     [,2]   [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 206397.7 193652.9 187342 200102.9 230972.5 250913.7 252763.6 230883.8
[2,] 696580.3 687364.1 728303 741036.1 795633.5 824552.3 819199.4 828245.2
         [,9]    [,10]    [,11]    [,12]
[1,] 251808.9 240024.8 242743.7 264653.2
[2,] 838272.1 831937.2 798276.3 843713.8

$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(192528, 276217, 451489, 766753, 1080925, 188996,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5ywjj1312382321.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]
[1,] 174414.0 188996.0 265306.0 334241.0 394185.0 486832.0 614128.0 759373.0
[2,] 177137.5 198930.5 278943.5 343370.0 406043.5 525949.5 670400.0 814825.0
[3,] 187158.0 242541.5 316121.0 380649.5 464164.0 585910.5 733951.0 890262.5
[4,] 192528.0 254400.0 324306.5 394285.0 479624.0 599532.0 747586.5 903884.0
[5,] 192528.0 259853.0 330709.0 403314.0 494239.0 614128.0 755841.0 908460.0
          [,9]   [,10] [,11]
[1,]  906511.0 1060884    NA
[2,]  967446.5 1124532    NA
[3,] 1055528.0 1218156    NA
[4,] 1074605.5 1250782    NA
[5,] 1086378.0 1253584    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]   [,6]     [,7]     [,8]
[1,] 180138.3 217241.5 295430.6 357426.8 430603.4 552349 698745.7 849642.1
[2,] 194177.7 267841.5 336811.4 403872.2 497724.6 619472 769156.3 930882.9
        [,9]   [,10] [,11]
[1,] 1006652 1160573    NA
[2,] 1104404 1275740    NA

$out
[1] 221698 219771

$group
[1] 1 1

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

Warning message:
In bxp(list(stats = c(174414, 177137.5, 187158, 192528, 192528,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/656xj1312382321.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,] 305688.7  871888.0 439603.0
[2,] 330901.5  937374.5 487760.2
[3,] 351046.5 1018337.0 510069.4
[4,] 363939.6 1073731.0 525489.8
[5,] 369922.8 1077421.0 538276.2

$n
[1] 12 12 12

$conf
         [,1]      [,2]     [,3]
[1,] 335977.6  956143.9 492860.7
[2,] 366115.3 1080530.1 527278.1

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(305688.696172084, 330901.493576298, 351046.464219944,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1yiej1312382321.ps tmp/1yiej1312382321.png",intern=TRUE))
character(0)
> try(system("convert tmp/23b1n1312382321.ps tmp/23b1n1312382321.png",intern=TRUE))
character(0)
> try(system("convert tmp/3o4vz1312382321.ps tmp/3o4vz1312382321.png",intern=TRUE))
character(0)
> try(system("convert tmp/4s70e1312382321.ps tmp/4s70e1312382321.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ywjj1312382321.ps tmp/5ywjj1312382321.png",intern=TRUE))
character(0)
> try(system("convert tmp/656xj1312382321.ps tmp/656xj1312382321.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.368   0.299   1.670