R version 2.10.1 (2009-12-14)
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(285708,905858,225733,405481,845758,805651,395747,695853,175625,405534,965639,575634,576023,566089,336141,26271,586226,376484,176583,287042,997142,207694,418003,838258,848182,658215,208304,398599,438399,578393,988390,958304,318251,78307,408520,748640,258520,518618,388588,238842,328957,499266,109011,168896,798921,878732,897576,518317,228370,758167,658491,518170,398212,498286,78136,647990,357927,698061,407932,637934,397784,217980,47737,467672,67651,167524,687406,367345,157553,887453,227566,817279,697059,997185,847075,547122,996977,346998,967154,547097,586853,46728,236883,36784,277085,446998,586725,496845,86765,146966,197113,657096,337200,17273,457284,507696,547628,157435,67793,267631,518397,918560,918895,429509,289569,9010172,1810617,7111400,1611919,9712714,2913310,1013816,6714518,2414721,9114534,8214993,6215159,515612,9415340,3715267)
> 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,] 285708 576023 848182 258520 228370 397784 697059 277085  547628 1611919
 [2,] 905858 566089 658215 518618 758167 217980 997185 446998  157435 9712714
 [3,] 225733 336141 208304 388588 658491  47737 847075 586725   67793 2913310
 [4,] 405481  26271 398599 238842 518170 467672 547122 496845  267631 1013816
 [5,] 845758 586226 438399 328957 398212  67651 996977  86765  518397 6714518
 [6,] 805651 376484 578393 499266 498286 167524 346998 146966  918560 2414721
 [7,] 395747 176583 988390 109011  78136 687406 967154 197113  918895 9114534
 [8,] 695853 287042 958304 168896 647990 367345 547097 657096  429509 8214993
 [9,] 175625 997142 318251 798921 357927 157553 586853 337200  289569 6215159
[10,] 405534 207694  78307 878732 698061 887453  46728  17273 9010172  515612
[11,] 965639 418003 408520 897576 407932 227566 236883 457284 1810617 9415340
[12,] 575634 838258 748640 518317 637934 817279  36784 507696 7111400 3715267
      [,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/www/rcomp/tmp/1s0hv1274454117.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/rcomp/tmp/2s0hv1274454117.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/rcomp/tmp/3s0hv1274454117.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/rcomp/tmp/4lrhg1274454117.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,] 228370 157435  47737.0  26271.0  67651 146966  78136.0 168896.0 157553.0
[2,] 277085 446998 208304.0 267631.0 328957 346998 176583.0 367345.0 289569.0
[3,] 472706 612152 362364.5 436576.5 478398 498776 541576.5 597543.5 347563.5
[4,] 697059 905858 658491.0 518170.0 845758 805651 967154.0 695853.0 798921.0
[5,] 848182 997185 847075.0 547122.0 996977 918560 988390.0 958304.0 997142.0
      [,10]    [,11]  [,12]
[1,]  17273 227566.0 507696
[2,]  78307 407932.0 518317
[3,] 460573 437643.5 693287
[4,] 878732 965639.0 838258
[5,] 887453 965639.0 838258

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 262870.2 382887.2 137433.1 311397.2 220183.6 269614.7 146575.7 433407.8
[2,] 682541.8 841416.8 587295.9 561755.8 736612.4 727937.3 936577.3 761679.2
          [,9]     [,10]    [,11]    [,12]
[1,]  93070.93  60648.76 158990.8 533431.7
[2,] 602056.07 860497.24 716296.2 853142.3

$out
 [1] 1611919 9712714 2913310 1013816 6714518 2414721 9114534 8214993 6215159
[10] 9010172 1810617 9415340   36784 7111400 3715267

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

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

Warning message:
In bxp(list(stats = c(228370, 277085, 472706, 697059, 848182, 157435,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5lrhg1274454117.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]
[1,] 175625.0  26271.0  78307 109011.0  78136.0  47737.0  36784.0  17273.0
[2,] 340727.5 247368.0 358425 248681.0 378069.5 162538.5 291940.5 172039.5
[3,] 490584.0 397243.5 508396 443927.0 508228.0 297455.5 566987.5 392099.0
[4,] 825704.5 581124.5 798411 658769.5 653240.5 577539.0 907114.5 502270.5
[5,] 965639.0 997142.0 988390 897576.0 758167.0 887453.0 997185.0 657096.0
          [,9]   [,10] [,11]
[1,]   67793.0  515612    NA
[2,]  278600.0 2013320    NA
[3,]  533012.5 4965213    NA
[4,] 1364756.0 8664764    NA
[5,] 1810617.0 9712714    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 269382.7 245014.9 307715.4 256882.9 382720.7 108171.0 286402.5 241478.4
[2,] 711785.3 549472.1 709076.6 630971.1 633735.3 486740.0 847572.5 542719.6
           [,9]   [,10] [,11]
[1,]   37609.46 1931445    NA
[2,] 1028415.54 7998981    NA

$out
[1] 9010172 7111400

$group
[1] 9 9

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

Warning message:
In bxp(list(stats = c(175625, 340727.5, 490584, 825704.5, 965639,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6wjg11274454117.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,]  257134.3  987545 212465.8
[2,]  751234.7 2566664 345418.5
[3,] 2097018.6 6860742 415960.2
[4,] 2742967.6 9014648 494854.4
[5,] 2900253.2 9555279 540544.2

$n
[1] 12 12 12

$conf
        [,1]    [,2]     [,3]
[1,] 1188576 3919772 347801.5
[2,] 3005461 9801711 484119.0

$out
[1] 773373.8 722910.5

$group
[1] 3 3

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

Warning message:
In bxp(list(stats = c(257134.259314057, 751234.74157378, 2097018.64381904,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1s0hv1274454117.ps tmp/1s0hv1274454117.png",intern=TRUE))
character(0)
> try(system("convert tmp/2s0hv1274454117.ps tmp/2s0hv1274454117.png",intern=TRUE))
character(0)
> try(system("convert tmp/3s0hv1274454117.ps tmp/3s0hv1274454117.png",intern=TRUE))
character(0)
> try(system("convert tmp/4lrhg1274454117.ps tmp/4lrhg1274454117.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lrhg1274454117.ps tmp/5lrhg1274454117.png",intern=TRUE))
character(0)
> try(system("convert tmp/6wjg11274454117.ps tmp/6wjg11274454117.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.420   1.280   2.427