R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(562674,599000,668516,597798,579889,668233,499232,215187,555813,586935,546136,571111,634712,639283,712182,621557,621000,675989,501322,220286,560727,602530,626379,605508,646783,658442,712906,687714,723916,707183,629000,237530,613296,730444,734925,651812,676155,748183,810681,729363,701108,790079,594621,230716,617189,691389,701067,705777,747636,773392,813788,766713,728875,749197,680954,241424,680234,708326,694238,772071,795337,788421,889968,797393,751000,821255,691605,290655,727147,868355,812390,799556,843038,847000,941952,804309,840307,871528,656330,370508,742000,847152,731675,898527,778139,856075,938833,813023,783417,828110,657311,310032,780000,860000,780000,807993,895217,856075,893268,875000,835088,934595,832500,300000,791443,900000,781729,880000,875024,992968,976804,968697,871675,1006852,832037,345587,849528,913871,868746,993733)
> 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,] 562674 634712 646783 676155 747636 795337 843038 778139 895217  875024
 [2,] 599000 639283 658442 748183 773392 788421 847000 856075 856075  992968
 [3,] 668516 712182 712906 810681 813788 889968 941952 938833 893268  976804
 [4,] 597798 621557 687714 729363 766713 797393 804309 813023 875000  968697
 [5,] 579889 621000 723916 701108 728875 751000 840307 783417 835088  871675
 [6,] 668233 675989 707183 790079 749197 821255 871528 828110 934595 1006852
 [7,] 499232 501322 629000 594621 680954 691605 656330 657311 832500  832037
 [8,] 215187 220286 237530 230716 241424 290655 370508 310032 300000  345587
 [9,] 555813 560727 613296 617189 680234 727147 742000 780000 791443  849528
[10,] 586935 602530 730444 691389 708326 868355 847152 860000 900000  913871
[11,] 546136 626379 734925 701067 694238 812390 731675 780000 781729  868746
[12,] 571111 605508 651812 705777 772071 799556 898527 807993 880000  993733
      [,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/html/freestat/rcomp/tmp/1f3s01274461627.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/freestat/rcomp/tmp/28v931274461627.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/freestat/rcomp/tmp/38v931274461627.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/freestat/rcomp/tmp/4cdqr1274461627.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]
[1,] 562674.0 599000.0 668516 597798 579889.0  668233 499232.0 215187.0
[2,] 646783.0 658442.0 712906 687714 701108.0  707183 594621.0 230716.0
[3,] 762887.5 780906.5 851878 782053 739937.5  805667 656820.5 266039.5
[4,] 843038.0 856075.0 938833 813023 835088.0  871528 691605.0 310032.0
[5,] 895217.0 992968.0 976804 968697 871675.0 1006852 832500.0 370508.0
         [,9]  [,10]  [,11]    [,12]
[1,] 555813.0 586935 626379 571111.0
[2,] 613296.0 691389 694238 651812.0
[3,] 703690.5 788798 733300 785813.5
[4,] 780000.0 868355 781729 880000.0
[5,] 849528.0 913871 868746 993733.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,] 664830.7 682161.2 738995.9 719443.6 672995.8 723553.7 608363.4 226410.1
[2,] 860944.3 879651.8 964760.1 844662.4 806879.2 887780.3 705277.6 305668.9
         [,9]    [,10]  [,11]    [,12]
[1,] 620398.5 700378.7 689586 671801.7
[2,] 786982.5 877217.3 777014 899825.3

$out
[1] 546136

$group
[1] 11

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

Warning message:
In bxp(list(stats = c(562674, 646783, 762887.5, 843038, 895217,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5cdqr1274461627.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,] 499232.0 501322.0 613296.0 594621.0 680234.0 691605.0 656330.0 778139.0
[2,] 550974.5 581628.5 637891.5 646672.0 687596.0 739073.5 736837.5 779069.5
[3,] 575500.0 621278.5 673078.0 701087.5 738255.5 796365.0 841672.5 795705.0
[4,] 598399.0 636997.5 718411.0 738773.0 769392.0 816822.5 859340.0 842092.5
[5,] 668516.0 712182.0 734925.0 810681.0 813788.0 889968.0 941952.0 860000.0
         [,9]     [,10] [,11]
[1,] 781729.0  832037.0    NA
[2,] 811971.5  859137.0    NA
[3,] 865537.5  894447.5    NA
[4,] 894242.5  984886.0    NA
[5,] 934595.0 1006852.0    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 553869.4 596024.3 636352.5 659079.6 700947.8 760903.2 785798.3 766959.8
[2,] 597130.6 646532.7 709803.5 743095.4 775563.2 831826.8 897546.7 824450.2
         [,9]    [,10] [,11]
[1,] 828013.1 837092.5    NA
[2,] 903061.9 951802.5    NA

$out
 [1] 215187 220286 237530 230716 241424 290655 370508 938833 657311 310032
[11] 300000 345587

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

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

Warning message:
In bxp(list(stats = c(499232, 550974.5, 575500, 598399, 668516,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6cdqr1274461627.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,]  93004.92 291786.0  75104.50
[2,]  98212.22 301001.5  99222.38
[3,] 111436.52 329739.5 149608.88
[4,] 117037.68 354759.0 174957.88
[5,] 135681.43 422622.0 196695.00

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 102850.1 305220.3 115065.4
[2,] 120022.9 354258.7 184152.3

$out
[1]  55026.3 155321.0

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(93004.9218855886, 98212.2247375678, 111436.522468122,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1f3s01274461627.ps tmp/1f3s01274461627.png",intern=TRUE))
character(0)
> try(system("convert tmp/28v931274461627.ps tmp/28v931274461627.png",intern=TRUE))
character(0)
> try(system("convert tmp/38v931274461627.ps tmp/38v931274461627.png",intern=TRUE))
character(0)
> try(system("convert tmp/4cdqr1274461627.ps tmp/4cdqr1274461627.png",intern=TRUE))
character(0)
> try(system("convert tmp/5cdqr1274461627.ps tmp/5cdqr1274461627.png",intern=TRUE))
character(0)
> try(system("convert tmp/6cdqr1274461627.ps tmp/6cdqr1274461627.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.631   1.234   1.838