R version 2.12.0 (2010-10-15)
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(00.521505,00.424828,00.425031,00.477194,00.828021,00.615619,00.366627,00.430888,00.281029,00.464625,00.269395,00.577905,00.566115,00.507758,00.750718,00.680840,00.766109,00.456147,00.497750,00.419327,00.609551,00.457337,00.570548,00.347900,00.387499,00.582429,00.239103,00.236745,00.262616,00.424093,00.365275,00.375076,00.409006,00.389168,00.240261,00.158950,00.439337,00.509468,00.374347,00.433983,00.413056,00.328893,00.518665,00.548650,00.546911,00.496349,00.530893,00.595776,00.557058,00.573133,00.500542,00.543127,00.559366,00.691169,00.440349,00.567666,00.596911,00.473554,00.592394,00.597556,00.633413,00.605712,00.704611,00.480526,00.702686,00.700902,00.603085,00.698092,00.597656,00.802342,00.601711,00.599313,00.602563,00.701663,00.499571,00.498092,00.497569,00.600183,00.333954,00.274437,00.320943,00.540667,00.405021,00.288596,00.327594,00.313261,00.257556,00.213839,00.186186,00.159271)
> 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] 90
> (np <- floor(n / par1))
[1] 7
> 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] 8 8 8 8 8 8 7 7 7 7 7 7
> arr
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
 [1,] 0.521505 0.566115 0.387499 0.439337 0.557058 0.633413 0.602563 0.327594
 [2,] 0.424828 0.507758 0.582429 0.509468 0.573133 0.605712 0.701663 0.313261
 [3,] 0.425031 0.750718 0.239103 0.374347 0.500542 0.704611 0.499571 0.257556
 [4,] 0.477194 0.680840 0.236745 0.433983 0.543127 0.480526 0.498092 0.213839
 [5,] 0.828021 0.766109 0.262616 0.413056 0.559366 0.702686 0.497569 0.186186
 [6,] 0.615619 0.456147 0.424093 0.328893 0.691169 0.700902 0.600183 0.159271
 [7,] 0.366627 0.497750 0.365275 0.518665 0.440349 0.603085 0.333954       NA
 [8,] 0.430888 0.419327 0.375076 0.548650 0.567666 0.698092 0.274437       NA
 [9,] 0.281029 0.609551 0.409006 0.546911 0.596911 0.597656 0.320943       NA
[10,] 0.464625 0.457337 0.389168 0.496349 0.473554 0.802342 0.540667       NA
[11,] 0.269395 0.570548 0.240261 0.530893 0.592394 0.601711 0.405021       NA
[12,] 0.577905 0.347900 0.158950 0.595776 0.597556 0.599313 0.288596       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/1pnkd1293625551.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/2iejy1293625551.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/3iejy1293625551.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/4sn001293625551.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]
[1,] 0.3275940 0.3132610 0.2391030 0.2138390 0.1861860 0.159271 0.3339540
[2,] 0.4134180 0.4662930 0.3159515 0.3353640 0.3378360 0.376493 0.3659510
[3,] 0.5392815 0.5413005 0.4623010 0.4788600 0.5284675 0.528165 0.4403490
[4,] 0.5843390 0.5940705 0.6025765 0.5206095 0.7343975 0.653394 0.5082075
[5,] 0.6334130 0.7016630 0.7507180 0.6808400 0.8280210 0.700902 0.6030850
          [,8]      [,9]    [,10]    [,11]    [,12]
[1,] 0.2744370 0.2810290 0.389168 0.240261 0.158950
[2,] 0.3972015 0.3649745 0.460981 0.337208 0.318248
[3,] 0.4308880 0.5469110 0.473554 0.530893 0.577905
[4,] 0.5581580 0.5972835 0.518508 0.581471 0.596666
[5,] 0.6980920 0.6095510 0.540667 0.601711 0.599313

$n
 [1] 8 8 8 8 8 8 7 7 7 7 7 7

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.4438026 0.4699222 0.3021882 0.3753792 0.3069426 0.3734841 0.3553957
[2,] 0.6347604 0.6126788 0.6224138 0.5823408 0.7499924 0.6828459 0.5253023
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.3347674 0.4081798 0.4391998 0.3850231 0.4116383
[2,] 0.5270086 0.6856422 0.5079082 0.6767629 0.7441717

$out
[1] 0.802342

$group
[1] 10

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

Warning message:
In bxp(list(stats = c(0.327594, 0.413418, 0.5392815, 0.584339, 0.633413,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5sn001293625551.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,] 0.2693950 0.3479000 0.1589500 0.3288930 0.4403490 0.4805260 0.2744370
[2,] 0.3957275 0.4567420 0.2396820 0.4235195 0.5218345 0.6005120 0.3274485
[3,] 0.4477565 0.5369365 0.3701755 0.5029085 0.5635160 0.6195625 0.4978305
[4,] 0.5497050 0.6451955 0.3990870 0.5389020 0.5946525 0.7017940 0.5704250
[5,] 0.6156190 0.7661090 0.5824290 0.5957760 0.6911690 0.8023420 0.7016630
          [,8]
[1,] 0.1592710
[2,] 0.1861860
[3,] 0.2356975
[4,] 0.3132610
[5,] 0.3275940

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.3775263 0.4509816 0.2974698 0.4502818 0.5303032 0.5733671 0.3870073
[2,] 0.5179867 0.6228914 0.4428812 0.5555352 0.5967288 0.6657579 0.6086537
          [,8]
[1,] 0.1537300
[2,] 0.3176650

$out
[1] 0.828021

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(0.269395, 0.3957275, 0.4477565, 0.549705,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6sn001293625551.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,] 0.09834873 0.2691310 0.0575270
[2,] 0.12612315 0.3449860 0.1334669
[3,] 0.14747154 0.4184145 0.1836833
[4,] 0.18552112 0.4893080 0.2392383
[5,] 0.23287101 0.6418350 0.3430958

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.1203797 0.3525883 0.1354402
[2,] 0.1745634 0.4842407 0.2319263

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.0983487329332281, 0.126123147809035, 0.147471541457361,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> try(system("convert tmp/1pnkd1293625551.ps tmp/1pnkd1293625551.png",intern=TRUE))
character(0)
> try(system("convert tmp/2iejy1293625551.ps tmp/2iejy1293625551.png",intern=TRUE))
character(0)
> try(system("convert tmp/3iejy1293625551.ps tmp/3iejy1293625551.png",intern=TRUE))
character(0)
> try(system("convert tmp/4sn001293625551.ps tmp/4sn001293625551.png",intern=TRUE))
character(0)
> try(system("convert tmp/5sn001293625551.ps tmp/5sn001293625551.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sn001293625551.ps tmp/6sn001293625551.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.370   0.990   2.359