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(122.6,115.4,109,129.1,102.8,96.2,127.7,128.9,126.5,119.8,113.2,114.1,134.1,130,121.8,132.1,105.3,103,117.1,126.3,138.1,119.5,138,135.5,178.6,162.2,176.9,204.9,132.2,142.5,164.3,174.9,175.4,143,158.7,155.4,176.6,163.3,178.9,182.7,158.9,115.5,169.4,168,159.8,129,157.9,169.5,169.1,183.6,168.9,186.2,227.1,126.4,169.3,175,133.9,110.1,104.3,108.7,112.6)
> 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] 61
> (np <- floor(n / par1))
[1] 5
> 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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]
 [1,] 122.6 134.1 178.6 176.6 169.1 112.6
 [2,] 115.4 130.0 162.2 163.3 183.6    NA
 [3,] 109.0 121.8 176.9 178.9 168.9    NA
 [4,] 129.1 132.1 204.9 182.7 186.2    NA
 [5,] 102.8 105.3 132.2 158.9 227.1    NA
 [6,]  96.2 103.0 142.5 115.5 126.4    NA
 [7,] 127.7 117.1 164.3 169.4 169.3    NA
 [8,] 128.9 126.3 174.9 168.0 175.0    NA
 [9,] 126.5 138.1 175.4 159.8 133.9    NA
[10,] 119.8 119.5 143.0 129.0 110.1    NA
[11,] 113.2 138.0 158.7 157.9 104.3    NA
[12,] 114.1 135.5 155.4 169.5 108.7    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/1q6m61258145479.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/2kvc91258145479.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/34t931258145479.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/4wvom1258145479.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] [,12]
[1,] 112.6 115.4 109.0 129.1 102.8  96.2 117.1 126.3 126.5 110.1 104.3 108.7
[2,] 122.6 130.0 121.8 132.1 105.3 103.0 127.7 128.9 133.9 119.5 113.2 114.1
[3,] 151.6 162.2 168.9 182.7 132.2 115.5 164.3 168.0 138.1 119.8 138.0 135.5
[4,] 176.6 163.3 176.9 186.2 158.9 126.4 169.3 174.9 159.8 129.0 157.9 155.4
[5,] 178.6 183.6 178.9 204.9 227.1 142.5 169.4 175.0 175.4 143.0 158.7 169.5

$n
 [1] 6 5 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]     [,2]     [,3]     [,4]      [,5]      [,6]     [,7]     [,8]
[1,] 116.7683 138.6703 129.9665 144.4731  94.32638  98.96562 134.9055 135.4965
[2,] 186.4317 185.7297 207.8335 220.9269 170.07362 132.03438 193.6945 200.5035
         [,9]    [,10]    [,11]    [,12]
[1,] 119.7991 113.0873 106.4151 106.3175
[2,] 156.4009 126.5127 169.5849 164.6825

$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(112.6, 122.6, 151.6, 176.6, 178.6, 115.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5d4l71258145479.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]
[1,]  96.2 103.00 132.20 157.90 104.30 112.6
[2,] 111.1 118.30 149.20 158.40 118.25 112.6
[3,] 117.6 128.15 163.25 165.65 169.00 112.6
[4,] 127.1 134.80 176.15 173.05 179.30 112.6
[5,] 129.1 138.10 204.90 182.70 227.10 112.6

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]  [,6]
[1,] 110.3023 120.6242 150.9579 158.9680 141.1547 112.6
[2,] 124.8977 135.6758 175.5421 172.3320 196.8453 112.6

$out
[1] 115.5 129.0

$group
[1] 4 4

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

Warning message:
In bxp(list(stats = c(96.2, 111.1, 117.6, 127.1, 129.1, 103, 118.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/63t481258145479.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,] 12.41680 32.90  9.500
[2,] 22.55045 48.80 29.600
[3,] 25.62243 57.60 43.150
[4,] 31.19584 69.05 51.425
[5,] 34.30000 75.80 55.100

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 21.67921 48.36384 33.19547
[2,] 29.56565 66.83616 53.10453

$out
[1]  51.12253 124.30000

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(12.4168031312412, 22.5504530375558, 25.6224319107185,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1q6m61258145479.ps tmp/1q6m61258145479.png")
> system("convert tmp/2kvc91258145479.ps tmp/2kvc91258145479.png")
> system("convert tmp/34t931258145479.ps tmp/34t931258145479.png")
> system("convert tmp/4wvom1258145479.ps tmp/4wvom1258145479.png")
> system("convert tmp/5d4l71258145479.ps tmp/5d4l71258145479.png")
> system("convert tmp/63t481258145479.ps tmp/63t481258145479.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.052   0.810   1.385