R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(36.80,35.40,33.00,28.73,26.70,26.46,24.60,28.00,31.60,33.50,34.50,35.00,34.76,33.50,32.74,34.40,31.93,29.24,25.75,26.03,26.08,23.80,20.61,19.70,18.18,19.60,20.60,20.03,23.00,23.60,22.56,22.55,23.75,24.92,24.50,30.58,28.07,27.70,27.00,25.23,26.86,25.60,24.55,23.96,23.50,23.64,21.55,21.05,21.89,21.98,21.45,22.15,22.58,23.80,23.30,22.38,23.00,21.96,22.40,20.80,20.40,16.00,12.78,9.75,7.50,11.24,12.24,12.75,12.52,14.49,14.21,14.32,22.15,22.58,23.80,23.30,22.38,23.00,21.96,22.40,20.80,20.40,16.00,12.78,9.75,7.50,11.24,12.24,12.75,12.52,14.49,14.21,14.32)
> 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] 93
> (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 8 8 8 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 36.80 34.76 18.18 28.07 21.89 20.40 22.15  9.75
 [2,] 35.40 33.50 19.60 27.70 21.98 16.00 22.58  7.50
 [3,] 33.00 32.74 20.60 27.00 21.45 12.78 23.80 11.24
 [4,] 28.73 34.40 20.03 25.23 22.15  9.75 23.30 12.24
 [5,] 26.70 31.93 23.00 26.86 22.58  7.50 22.38 12.75
 [6,] 26.46 29.24 23.60 25.60 23.80 11.24 23.00 12.52
 [7,] 24.60 25.75 22.56 24.55 23.30 12.24 21.96 14.49
 [8,] 28.00 26.03 22.55 23.96 22.38 12.75 22.40 14.21
 [9,] 31.60 26.08 23.75 23.50 23.00 12.52 20.80 14.32
[10,] 33.50 23.80 24.92 23.64 21.96 14.49 20.40    NA
[11,] 34.50 20.61 24.50 21.55 22.40 14.21 16.00    NA
[12,] 35.00 19.70 30.58 21.05 20.80 14.32 12.78    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/1t1hc1244400311.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/2a5rl1244400311.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/3lpqw1244400311.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/4sk8j1244400311.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]   [,9] [,10]  [,11]
[1,]  9.750  7.50 11.240  9.750  7.500 11.24 12.240 12.750 12.520 20.40 14.210
[2,] 19.290 17.80 16.690 16.135 17.565 17.76 18.225 18.295 17.560 21.18 18.305
[3,] 22.020 22.28 22.625 22.725 22.790 23.70 22.930 22.475 23.250 23.64 21.550
[4,] 31.415 30.60 29.870 26.980 26.780 26.03 24.575 24.995 24.915 24.36 23.450
[5,] 36.800 35.40 33.000 34.400 31.930 29.24 25.750 28.000 31.600 24.92 24.500
      [,12]
[1,] 12.780
[2,] 17.010
[3,] 20.800
[4,] 25.815
[5,] 35.000

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 15.2468 15.12974 15.26246 16.66683 17.64237 19.08026 19.3828 18.73228
[2,] 28.7932 29.43026 29.98754 28.78317 27.93763 28.31974 26.4772 26.21772
         [,9]    [,10]    [,11]    [,12]
[1,] 19.14139 21.74096 18.47749 15.54180
[2,] 27.35861 25.53904 24.62251 26.05820

$out
[1] 33.50 14.49 34.50

$group
[1] 10 10 11

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

Warning message:
In bxp(list(stats = c(9.75, 19.29, 22.02, 31.415, 36.8, 7.5, 17.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/509db1244400311.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,] 24.60 19.700 18.180 21.05 20.800  9.750 20.400  7.50
[2,] 27.35 24.775 20.315 23.57 21.925 11.740 20.600 11.24
[3,] 32.30 27.660 22.780 24.89 22.265 12.765 22.265 12.52
[4,] 34.75 33.120 24.125 26.93 22.790 14.405 22.790 14.21
[5,] 36.80 34.760 24.920 28.07 23.800 16.000 23.800 14.49

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 28.92481 23.85379 21.04223 23.35748 21.87047 11.54948 21.26613 10.9558
[2,] 35.67519 31.46621 24.51777 26.42252 22.65953 13.98052 23.26387 14.0842

$out
[1] 30.58 20.40  7.50 16.00 12.78

$group
[1] 3 6 6 7 7

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

Warning message:
In bxp(list(stats = c(24.6, 27.35, 32.3, 34.75, 36.8, 19.7, 24.775,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6lomg1244400311.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,] 5.008922 13.510  3.18000
[2,] 5.927858 18.505  4.80750
[3,] 7.269088 21.025  6.10125
[4,] 8.117747 24.540  9.29750
[5,] 9.166074 27.900 10.45000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 6.270265 18.27240 4.053331
[2,] 8.267911 23.77760 8.149169

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(5.00892186146737, 5.92785825464381, 7.26908804178794,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1t1hc1244400311.ps tmp/1t1hc1244400311.png")
> system("convert tmp/2a5rl1244400311.ps tmp/2a5rl1244400311.png")
> system("convert tmp/3lpqw1244400311.ps tmp/3lpqw1244400311.png")
> system("convert tmp/4sk8j1244400311.ps tmp/4sk8j1244400311.png")
> system("convert tmp/509db1244400311.ps tmp/509db1244400311.png")
> system("convert tmp/6lomg1244400311.ps tmp/6lomg1244400311.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.019   0.807   1.340