R version 2.12.1 (2010-12-16)
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(20.98,20.1,20.61,20.27,20.08,23.58,22.31,22.89,21.78,22.19,22.58,22.78,25.06,25.16,25.47,25.34,24.2,25.32,25.57,25.76,24.79,23.14,22.66,22.06,24.26,23.15,22.92,21.43,21.56,23.48,24.35,24.83,24.19,23.58,23.58,24.35,27.18,25.69,24.81,23.26,23.49,26.86,27.12,27.66,26.26,25.51,24.63,23.57,27.63,25.85,26.09,24.47,24.19,25.09,25.26,25.58,24.76,25.02,24.24,24.14,28.69,26.74,26.48,24.45,23.88,26.58,26.23,28.63,26.81,26.56,26.64,26.8,28.37,27.13,28.44,28.62,27.28,31.32,31.26,31.41,31.76,32.72,32.15,33.62,35.97,33.78,33.77,32.75,32.55,33.22,32.88,31.56,30.27,28.65,27.89,27.07,30.8,28.38,27.5,28,28.02,29.2,27.59,27.22,27.16,26.31,25.67,26.41,28.34,25.43,23.72,23.33,23.8,27.7,26.28,27.51,27.93,28.76,28.65,29.52,31.23,27.9,27.87,27.52,27.59,31.2,30.22,30.62,31.52,30.59,31.42,31.95)
> 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] 132
> (np <- floor(n / par1))
[1] 11
> 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] 11 11 11 11 11 11 11 11 11 11 11 11
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 20.98 25.06 24.26 27.18 27.63 28.69 28.37 35.97 30.80 28.34 31.23    NA
 [2,] 20.10 25.16 23.15 25.69 25.85 26.74 27.13 33.78 28.38 25.43 27.90    NA
 [3,] 20.61 25.47 22.92 24.81 26.09 26.48 28.44 33.77 27.50 23.72 27.87    NA
 [4,] 20.27 25.34 21.43 23.26 24.47 24.45 28.62 32.75 28.00 23.33 27.52    NA
 [5,] 20.08 24.20 21.56 23.49 24.19 23.88 27.28 32.55 28.02 23.80 27.59    NA
 [6,] 23.58 25.32 23.48 26.86 25.09 26.58 31.32 33.22 29.20 27.70 31.20    NA
 [7,] 22.31 25.57 24.35 27.12 25.26 26.23 31.26 32.88 27.59 26.28 30.22    NA
 [8,] 22.89 25.76 24.83 27.66 25.58 28.63 31.41 31.56 27.22 27.51 30.62    NA
 [9,] 21.78 24.79 24.19 26.26 24.76 26.81 31.76 30.27 27.16 27.93 31.52    NA
[10,] 22.19 23.14 23.58 25.51 25.02 26.56 32.72 28.65 26.31 28.76 30.59    NA
[11,] 22.58 22.66 23.58 24.63 24.24 26.64 32.15 27.89 25.67 28.65 31.42    NA
[12,] 22.78 22.06 24.35 23.57 24.14 26.80 33.62 27.07 26.41 29.52 31.95    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/153d61322591947.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/2nde91322591947.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/39wmh1322591947.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/4dss31322591947.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]  [,10]
[1,] 20.980 23.150 20.610 20.270 20.080 23.480 22.310 22.890 21.780 22.190
[2,] 26.120 25.295 24.265 23.295 23.645 25.205 25.415 25.670 24.775 24.300
[3,] 28.340 25.850 26.090 24.470 24.190 26.860 26.280 27.510 26.810 26.310
[4,] 29.745 27.515 27.685 27.760 27.435 30.200 28.905 29.625 29.100 28.705
[5,] 31.230 28.380 28.440 32.750 32.550 33.220 32.880 31.560 31.760 32.720
     [,11]  [,12]
[1,] 22.58 22.060
[2,] 23.91 23.855
[3,] 25.67 26.410
[4,] 28.27 28.295
[5,] 32.15 33.620

$n
 [1] 11 11 11 11 11 11 11 11 11 11 11 11

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 26.61309 24.79242 24.46075 22.34293 22.38449 24.48044 24.61741 25.62589
[2,] 30.06691 26.90758 27.71925 26.59707 25.99551 29.23956 27.94259 29.39411
         [,9]    [,10]    [,11]    [,12]
[1,] 24.74962 24.21151 23.59295 24.29484
[2,] 28.87038 28.40849 27.74705 28.52516

$out
[1] 35.97 20.10 33.78 33.77

$group
[1] 1 2 2 3

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

Warning message:
In bxp(list(stats = c(20.98, 26.12, 28.34, 29.745, 31.23, 23.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/53b441322591947.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]  [,10]
[1,] 20.080 22.060 21.430 23.26 24.140 26.230 27.130 27.070 25.670 23.330
[2,] 20.440 23.670 23.035 24.10 24.355 26.355 28.405 29.460 26.785 24.615
[3,] 21.985 25.110 23.580 25.60 25.055 26.610 31.290 32.650 27.545 27.605
[4,] 22.680 25.405 24.305 26.99 25.715 26.805 31.955 33.495 28.200 28.495
[5,] 23.580 25.760 24.830 27.66 27.630 26.810 33.620 35.970 29.200 29.520
      [,11] [,12]
[1,] 27.520    NA
[2,] 27.885    NA
[3,] 30.605    NA
[4,] 31.325    NA
[5,] 31.950    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 20.96332 24.31865 23.00074 24.28185 24.43469 26.40475 29.67082 30.80961
[2,] 23.00668 25.90135 24.15926 26.91815 25.67531 26.81525 32.90918 34.49039
         [,9]    [,10]    [,11] [,12]
[1,] 26.89961 25.83531 29.03599    NA
[2,] 28.19039 29.37469 32.17401    NA

$out
[1] 28.69 24.45 23.88 28.63 30.80

$group
[1] 6 6 6 6 9

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

Warning message:
In bxp(list(stats = c(20.08, 20.44, 21.985, 22.68, 23.58, 22.06,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6l5i91322591947.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,] 2.790298  8.670 3.4200
[2,] 3.215474  9.860 3.5575
[3,] 3.355235 11.065 4.1400
[4,] 3.524097 12.820 4.4225
[5,] 3.939793 14.990 4.9950

$n
[1] 12 12 12

$conf
        [,1]      [,2]     [,3]
[1,] 3.21447  9.714924 3.745468
[2,] 3.49600 12.415076 4.534532

$out
[1] 2.22

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(2.79029845259346, 3.21547440503535, 3.35523512759515,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/153d61322591947.ps tmp/153d61322591947.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nde91322591947.ps tmp/2nde91322591947.png",intern=TRUE))
character(0)
> try(system("convert tmp/39wmh1322591947.ps tmp/39wmh1322591947.png",intern=TRUE))
character(0)
> try(system("convert tmp/4dss31322591947.ps tmp/4dss31322591947.png",intern=TRUE))
character(0)
> try(system("convert tmp/53b441322591947.ps tmp/53b441322591947.png",intern=TRUE))
character(0)
> try(system("convert tmp/6l5i91322591947.ps tmp/6l5i91322591947.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.860   0.368   2.206