R version 2.6.2 (2008-02-08)
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.

  Natural language support but running in an English locale

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

$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 487.2019 478.4074 496.6808 504.9052 518.8562 504.2271 442.2135
[2,] 555.7463 581.8441 578.7826 595.8532 614.7438 625.4488 620.4579 563.4195
         [,9]    [,10]    [,11]    [,12]
[1,] 448.0144 461.9856 464.1715 478.5596
[2,] 563.8806 573.0604 577.4345 587.3034

$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/rcomp/tmp/5seo41210412371.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,] 498.6620 557.2740 567.456 536.6620 500.8570 454.7240 431.3280 439.1480
[2,] 512.6480 569.2585 578.764 549.4655 509.5300 465.8450 445.4275 458.3890
[3,] 526.4945 577.0020 594.406 592.1975 551.0660 510.7385 471.0255 465.6165
[4,] 545.2555 604.8195 616.426 603.1085 567.3935 518.8745 493.4230 486.6025
[5,] 564.5910 625.8090 628.884 612.6130 579.7140 525.1920 512.9270 520.5640
        [,9]
[1,] 492.865
[2,] 492.865
[3,] 492.865
[4,] 492.865
[5,] 492.865

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 511.622 560.7824 577.2281 567.7306 524.6741 486.5514 449.1344 452.7481
[2,] 541.367 593.2216 611.5839 616.6644 577.4579 534.9256 492.9166 478.4849
        [,9]
[1,] 492.865
[2,] 492.865

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(498.662, 512.648, 526.4945, 545.2555, 564.591,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6f2021210412371.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/1zpo81210412371.ps tmp/1zpo81210412371.png")
> system("convert tmp/2o8671210412371.ps tmp/2o8671210412371.png")
> system("convert tmp/356651210412371.ps tmp/356651210412371.png")
> system("convert tmp/4uq0v1210412371.ps tmp/4uq0v1210412371.png")
> system("convert tmp/5seo41210412371.ps tmp/5seo41210412371.png")
> system("convert tmp/6f2021210412371.ps tmp/6f2021210412371.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.398   1.370   2.602