R version 2.12.1 (2010-12-16)
Copyright (C) 2010 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(5939520.00,89948768.00,80953652.00,85942882.00,8944937.00,82975432.00,24940816.00,21973899.00,37950221.00,45949881.00,85950373.00,48960313.00,81954506.00,24960419.00,65973338.00,22950513.00,54963528.00,90995659.00,91967517.00,28999053.00,96990529.00,38979852.00,81496957.00,74982424.00,70976192.00,90990000.00,12998850.00,92986156.00,67994976.00,91022206.00,87992489.00,421022698.00,11018942.00,79100042.00,65996442.00,51000620.00,12996871.00,44994249.00,99996135.00,91977037.00,63974211.00,15998036.00,65974265.00,33984410.00,45939098.00,67935827.00,66921032.00,89911836.00,71890975.00,72880342.00,28871286.00,41844334.00,82847667.00,24871401.00,3867451.00,99896846.00,41890361.00,45884264.00,69884586.00,95896400.00,39904491.00,81900399.00,27909863.00,88900470.00,89917101.00,2945005.00,4934411.00,61957264.00,31946515.00,3938309.00,52933321.00,21947613.00)
> 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] 72
> (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] 6 6 6 6 6 6 6 6 6 6 6 6
> arr
          [,1]     [,2]      [,3]     [,4]     [,5]     [,6] [,7]
 [1,]  5939520 81954506  70976192 12996871 71890975 39904491   NA
 [2,] 89948768 24960419  90990000 44994249 72880342 81900399   NA
 [3,] 80953652 65973338  12998850 99996135 28871286 27909863   NA
 [4,] 85942882 22950513  92986156 91977037 41844334 88900470   NA
 [5,]  8944937 54963528  67994976 63974211 82847667 89917101   NA
 [6,] 82975432 90995659  91022206 15998036 24871401  2945005   NA
 [7,] 24940816 91967517  87992489 65974265  3867451  4934411   NA
 [8,] 21973899 28999053 421022698 33984410 99896846 61957264   NA
 [9,] 37950221 96990529  11018942 45939098 41890361 31946515   NA
[10,] 45949881 38979852  79100042 67935827 45884264  3938309   NA
[11,] 85950373 81496957  65996442 66921032 69884586 52933321   NA
[12,] 48960313 74982424  51000620 89911836 95896400 21947613   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/rcomp/tmp/196a71305049469.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/rcomp/tmp/21b7l1305049469.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/rcomp/tmp/3p2dv1305049469.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/rcomp/tmp/4up2n1305049469.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,]  5939520 24960419 12998850 22950513 54963528  2945005  3867451 21973899
[2,] 12996871 44994249 27909863 41844334 54963528 15998036  4934411 28999053
[3,] 55440342 77390370 47422312 87421676 65984594 53923416 45457540 47970837
[4,] 71890975 89948768 80953652 91977037 82847667 90995659 87992489 99896846
[5,] 81954506 90990000 99996135 92986156 89917101 91022206 91967517 99896846
         [,9]    [,10]    [,11]    [,12]
[1,] 11018942  3938309 52933321 21947613
[2,] 31946515 38979852 65996442 48960313
[3,] 39920291 45917072 68402809 62991522
[4,] 45939098 67935827 81496957 89911836
[5,] 45939098 79100042 85950373 95896400

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

$conf
         [,1]      [,2]     [,3]      [,4]     [,5]      [,6]     [,7]     [,8]
[1,] 17451742  48393254 13207355  55084463 47998423   5547527 -8117602  2239471
[2,] 93428941 106387487 81637269 119758889 83970764 102299306 99032683 93702203
         [,9]    [,10]    [,11]    [,12]
[1,] 30894623 27239533 58404476 36576467
[2,] 48945959 64594612 78401142 89406577

$out
[1]   8944937 421022698  96990529

$group
[1] 5 8 9

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

Warning message:
In bxp(list(stats = c(5939520, 12996871, 55440341.5, 71890975, 81954506,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5nsk81305049469.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]
[1,]  5939520 22950513 11018942 12996871  3867451  2945005   NA
[2,] 23457358 33989452 58498531 39489330 35357810 13441012   NA
[3,] 47455097 70477881 75038117 64974238 57884425 35925503   NA
[4,] 84459157 86475082 91006103 78923832 77864004 71928832   NA
[5,] 89948768 96990529 92986156 99996135 99896846 89917101   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 19631767 46538833 60211195 46987897 38497064  9248816   NA
[2,] 75278427 94416929 89865039 82960579 77271786 62602190   NA

$out
[1] 421022698

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(5939520, 23457357.5, 47455097, 84459157, 89948768,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/64k911305049469.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,] 26037257 66029581 11479472
[2,] 27525304 71992215 21825840
[3,] 29554185 78493575 37524009
[4,] 37498396 87537243 56052531
[5,] 41084469 88100066 72551921

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 25005391 71403383 21912986
[2,] 34102979 85583767 53135032

$out
[1] 154417088  11846514 399048799  33017052

$group
[1] 1 1 2 2

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

Warning message:
In bxp(list(stats = c(26037256.7400171, 27525303.9002621, 29554185.1288594,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/196a71305049469.ps tmp/196a71305049469.png",intern=TRUE))
character(0)
> try(system("convert tmp/21b7l1305049469.ps tmp/21b7l1305049469.png",intern=TRUE))
character(0)
> try(system("convert tmp/3p2dv1305049469.ps tmp/3p2dv1305049469.png",intern=TRUE))
character(0)
> try(system("convert tmp/4up2n1305049469.ps tmp/4up2n1305049469.png",intern=TRUE))
character(0)
> try(system("convert tmp/5nsk81305049469.ps tmp/5nsk81305049469.png",intern=TRUE))
character(0)
> try(system("convert tmp/64k911305049469.ps tmp/64k911305049469.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.016   0.616   3.592