R version 2.7.0 (2008-04-22)
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(492.865,480.961,461.935,456.608,441.977,439.148,488.180,520.564,501.492,485.025,464.196,460.170,467.037,460.070,447.988,442.867,436.087,431.328,484.015,509.673,512.927,502.831,470.984,471.067,476.049,474.605,470.439,461.251,454.724,455.626,516.847,525.192,522.975,518.585,509.239,512.238,519.164,517.009,509.933,509.127,500.857,506.971,569.323,579.714,577.992,565.464,547.344,554.788,562.325,560.854,555.332,543.599,536.662,542.722,593.530,610.763,612.613,611.324,594.167,595.454,590.865,589.379,584.428,573.100,567.456,569.028,620.735,628.884,628.232,612.117,595.404,597.141,593.408,590.072,579.799,574.205,572.775,572.942,619.567,625.809,619.916,587.625,565.742,557.274,560.576,548.854,531.673,525.919,511.038,498.662,555.362,564.591,541.657,527.070,509.846,514.258,516.922)
> 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] 97
> (np <- floor(n / par1))
[1] 8
> 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] 9 8 8 8 8 8 8 8 8 8 8 8
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
 [1,] 492.865 467.037 476.049 519.164 562.325 590.865 593.408 560.576 516.922
 [2,] 480.961 460.070 474.605 517.009 560.854 589.379 590.072 548.854      NA
 [3,] 461.935 447.988 470.439 509.933 555.332 584.428 579.799 531.673      NA
 [4,] 456.608 442.867 461.251 509.127 543.599 573.100 574.205 525.919      NA
 [5,] 441.977 436.087 454.724 500.857 536.662 567.456 572.775 511.038      NA
 [6,] 439.148 431.328 455.626 506.971 542.722 569.028 572.942 498.662      NA
 [7,] 488.180 484.015 516.847 569.323 593.530 620.735 619.567 555.362      NA
 [8,] 520.564 509.673 525.192 579.714 610.763 628.884 625.809 564.591      NA
 [9,] 501.492 512.927 522.975 577.992 612.613 628.232 619.916 541.657      NA
[10,] 485.025 502.831 518.585 565.464 611.324 612.117 587.625 527.070      NA
[11,] 464.196 470.984 509.239 547.344 594.167 595.404 565.742 509.846      NA
[12,] 460.170 471.067 512.238 554.788 595.454 597.141 557.274 514.258      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/1t9fw1212310625.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/2tg4l1212310625.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/3vjp91212310625.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/4yk7z1212310625.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,] 467.037 460.0700 447.9880 442.8670 436.0870 431.3280 484.0150 509.6730
[2,] 492.865 477.7830 466.1870 458.9295 448.3505 447.3870 502.5135 522.8780
[3,] 519.164 532.9315 520.8030 517.5230 505.9475 502.8165 562.3425 572.1525
[4,] 562.325 575.1165 567.5655 558.3495 552.0590 555.8750 606.5485 618.2860
[5,] 593.408 590.0720 584.4280 574.2050 572.7750 572.9420 620.7350 628.8840
         [,9]    [,10]    [,11]    [,12]
[1,] 501.4920 485.0250 464.1960 460.1700
[2,] 517.9510 510.7080 490.1115 491.6525
[3,] 559.8245 546.2670 528.5950 534.5230
[4,] 616.2645 599.4745 579.9545 576.3640
[5,] 628.2320 612.1170 595.4040 597.1410

$n
 [1] 9 8 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 482.5817 478.5596 464.1715 461.9856 448.0144 442.2135 504.2271 518.8562
[2,] 555.7463 587.3034 577.4345 573.0604 563.8806 563.4195 620.4579 625.4488
         [,9]    [,10]    [,11]    [,12]
[1,] 504.9052 496.6808 478.4074 487.2019
[2,] 614.7438 595.8532 578.7826 581.8441

$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(467.037, 492.865, 519.164, 562.325, 593.408,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5o8kl1212310625.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,] 439.1480 431.3280 454.724 500.8570 536.6620 567.4560 557.2740 498.6620
[2,] 458.3890 445.4275 465.845 509.5300 549.4655 578.7640 572.8585 512.6480
[3,] 472.5785 469.0105 492.644 533.2540 577.9275 593.1345 583.7120 529.3715
[4,] 490.5225 493.4230 517.716 567.3935 603.1085 616.4260 606.4875 552.1080
[5,] 520.5640 512.9270 525.192 579.7140 612.6130 628.8840 625.8090 564.5910
        [,9]
[1,] 516.922
[2,] 516.922
[3,] 516.922
[4,] 516.922
[5,] 516.922

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 457.9222 447.1194 468.9853 506.8621 553.4606 575.9566 568.3736 511.3735
[2,] 487.2348 490.9016 516.3027 559.6459 602.3944 610.3124 599.0504 547.3695
        [,9]
[1,] 516.922
[2,] 516.922

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(439.148, 458.389, 472.5785, 490.5225, 520.564,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6dcf01212310625.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,] 47.74673 119.211 64.87375
[2,] 50.50759 126.916 76.03813
[3,] 52.05310 131.273 90.42425
[4,] 54.07548 136.704 92.97950
[5,] 56.46122 141.614 97.79200

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 50.42576 126.8086 82.69718
[2,] 53.68044 135.7374 98.15132

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(47.746733027274, 50.5075889737488, 52.053101869091,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1t9fw1212310625.ps tmp/1t9fw1212310625.png")
> system("convert tmp/2tg4l1212310625.ps tmp/2tg4l1212310625.png")
> system("convert tmp/3vjp91212310625.ps tmp/3vjp91212310625.png")
> system("convert tmp/4yk7z1212310625.ps tmp/4yk7z1212310625.png")
> system("convert tmp/5o8kl1212310625.ps tmp/5o8kl1212310625.png")
> system("convert tmp/6dcf01212310625.ps tmp/6dcf01212310625.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.845   1.200   1.984