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(31.514,27.071,29.462,26.105,22.397,23.843,21.705,18.089,20.764,25.316,17.704,15.548,28.029,29.383,36.438,32.034,22.679,24.319,18.004,17.537,20.366,22.782,19.169,13.807,29.743,25.591,29.096,26.482,22.405,27.044,17.970,18.730,19.684,19.785,18.479,10.698,31.956,29.506,34.506,27.165,26.736,23.691,18.157,17.328,18.205,20.995,17.382,9.367,31.124,26.551,30.651,25.859,25.100,25.778,20.418,18.688,20.424,24.776,19.814,12.738,31.566,30.111,30.019,31.934,25.826,26.835,20.205,17.789,20.520,22.518,15.572,11.509,25.447,24.090,27.786,26.195,20.516,22.759,19.028,16.971,20.036,22.485,18.730,14.538,27.561,25.985,34.670,32.066,27.186,29.586,21.359,21.553,19.573,24.256,22.380,16.167,27.297,28.287)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> ylab = ''
> xlab = ''
> main = ''
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 98
> (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 9 8 8 8 8 8 8 8 8 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,] 31.514 28.029 29.743 31.956 31.124 31.566 25.447 27.561 27.297
 [2,] 27.071 29.383 25.591 29.506 26.551 30.111 24.090 25.985 28.287
 [3,] 29.462 36.438 29.096 34.506 30.651 30.019 27.786 34.670     NA
 [4,] 26.105 32.034 26.482 27.165 25.859 31.934 26.195 32.066     NA
 [5,] 22.397 22.679 22.405 26.736 25.100 25.826 20.516 27.186     NA
 [6,] 23.843 24.319 27.044 23.691 25.778 26.835 22.759 29.586     NA
 [7,] 21.705 18.004 17.970 18.157 20.418 20.205 19.028 21.359     NA
 [8,] 18.089 17.537 18.730 17.328 18.688 17.789 16.971 21.553     NA
 [9,] 20.764 20.366 19.684 18.205 20.424 20.520 20.036 19.573     NA
[10,] 25.316 22.782 19.785 20.995 24.776 22.518 22.485 24.256     NA
[11,] 17.704 19.169 18.479 17.382 19.814 15.572 18.730 22.380     NA
[12,] 15.548 13.807 10.698  9.367 12.738 11.509 14.538 16.167     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/18wrx1305863343.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/28oi01305863343.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/3yf3c1305863343.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/4b6ky1305863343.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]    [,8]    [,9]
[1,] 25.447 24.090 27.786 25.8590 20.5160 22.7590 17.9700 16.9710 19.5730
[2,] 27.561 25.985 29.279 26.1500 22.4010 23.7670 18.0805 17.4325 19.6285
[3,] 29.743 27.071 30.335 26.8235 23.8895 25.0485 19.6165 17.9390 20.2010
[4,] 31.514 29.383 34.588 31.9840 26.2810 26.9395 20.8885 18.7090 20.4720
[5,] 31.956 30.111 36.438 32.0660 27.1860 29.5860 21.7050 18.7300 20.7640
      [,10]   [,11]   [,12]
[1,] 19.785 15.5720  9.3670
[2,] 21.740 17.5430 11.1035
[3,] 22.650 18.6045 13.2725
[4,] 24.516 19.4915 15.0430
[5,] 25.316 22.3800 16.1670

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 27.66109 25.28139 27.36932 23.56454 21.72208 23.27630 18.04791 17.22593
[2,] 31.82491 28.86061 33.30068 30.08246 26.05692 26.82070 21.18509 18.65207
         [,9]    [,10]    [,11]    [,12]
[1,] 19.72981 21.09929 17.51604 11.07184
[2,] 20.67219 24.20071 19.69296 15.47316

$out
[1] 21.553 18.205

$group
[1] 8 9

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

Warning message:
In bxp(list(stats = c(25.447, 27.561, 29.743, 31.514, 31.956, 24.09,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5fvn01305863343.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]    [,8]   [,9]
[1,] 15.5480 13.8070 10.6980  9.3670 12.738 11.509 14.5380 16.1670 27.297
[2,] 19.4265 18.5865 18.6045 17.7695 20.116 18.997 18.8790 21.4560 27.297
[3,] 23.1200 22.7305 21.0950 22.3430 24.938 24.172 21.5005 25.1205 27.792
[4,] 26.5880 28.7060 26.7630 28.3355 26.205 30.065 24.7685 28.5735 28.287
[5,] 31.5140 36.4380 29.7430 34.5060 31.124 31.934 27.7860 34.6700 28.287

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 19.85359 18.11493 17.37385 17.52378 22.16077 19.12381 18.81426 21.87416
[2,] 26.38641 27.34607 24.81615 27.16222 27.71523 29.22019 24.18674 28.36684
         [,9]
[1,] 26.68594
[2,] 28.89806

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15.548, 19.4265, 23.12, 26.588, 31.514, 13.807,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6ajs61305863343.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.8150755 2.5590 0.791750
[2,] 1.7020743 5.0565 1.990125
[3,] 2.1558607 6.3580 3.240125
[4,] 2.4112474 6.8040 3.801750
[5,] 3.1636349 8.6520 5.786500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.832402 5.560954 2.413831
[2,] 2.479319 7.155046 4.066419

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.81507545662963, 1.70207429793347, 2.15586073836775,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> try(system("convert tmp/18wrx1305863343.ps tmp/18wrx1305863343.png",intern=TRUE))
character(0)
> try(system("convert tmp/28oi01305863343.ps tmp/28oi01305863343.png",intern=TRUE))
character(0)
> try(system("convert tmp/3yf3c1305863343.ps tmp/3yf3c1305863343.png",intern=TRUE))
character(0)
> try(system("convert tmp/4b6ky1305863343.ps tmp/4b6ky1305863343.png",intern=TRUE))
character(0)
> try(system("convert tmp/5fvn01305863343.ps tmp/5fvn01305863343.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ajs61305863343.ps tmp/6ajs61305863343.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.390   0.900   2.318