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(-9,-13,-18,-11,-9,-10,-13,-11,-5,-15,-6,-6,-3,-1,-3,-4,-6,0,-4,-2,-2,-6,-7,-6,-6,-3,-2,-5,-11,-11,-11,-10,-14,-8,-9,-5,-1,-2,-5,-4,-6,-2,-2,-2,-2,2,1,-8,-1,1,-1,2,2,1,-1,-2,-2,-1,-8,-4,-6,-3,-3,-7,-9,-11,-13,-11,-9,-17,-22,-25,-20)
> 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] 73
> (np <- floor(n / par1))
[1] 6
> 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] 7 6 6 6 6 6 6 6 6 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   -9   -3   -6   -1   -1   -6  -20
 [2,]  -13   -1   -3   -2    1   -3   NA
 [3,]  -18   -3   -2   -5   -1   -3   NA
 [4,]  -11   -4   -5   -4    2   -7   NA
 [5,]   -9   -6  -11   -6    2   -9   NA
 [6,]  -10    0  -11   -2    1  -11   NA
 [7,]  -13   -4  -11   -2   -1  -13   NA
 [8,]  -11   -2  -10   -2   -2  -11   NA
 [9,]   -5   -2  -14   -2   -2   -9   NA
[10,]  -15   -6   -8    2   -1  -17   NA
[11,]   -6   -7   -9    1   -8  -22   NA
[12,]   -6   -6   -5   -8   -4  -25   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/1xrwk1243892189.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/2se3b1243892189.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/39ofm1243892189.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/4muxl1243892189.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,] -9.0 -3.0   -5 -11.0 -11.0  -11 -13.0  -11 -14.0   -17  -9.0    -8
[2,] -7.5 -3.0   -5  -7.0  -9.0  -11 -13.0  -11  -9.0   -15  -9.0    -8
[3,] -6.0 -2.5   -3  -4.5  -7.5   -6  -7.5   -6  -3.5    -7  -7.5    -6
[4,] -2.0 -1.0   -2  -4.0  -6.0    0  -2.0   -2  -2.0    -1  -6.0    -5
[5,] -1.0  1.0   -1  -4.0  -6.0    1  -1.0   -2  -2.0     2  -6.0    -4

$n
 [1] 7 6 6 6 6 6 6 6 6 6 6 6

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]       [,6]        [,7]
[1,] -9.284511 -3.790065 -4.935097 -6.435097 -9.435097 -13.095355 -14.5953553
[2,] -2.715489 -1.209935 -1.064903 -2.564903 -5.564903   1.095355  -0.4046447
            [,8]      [,9]      [,10]     [,11]     [,12]
[1,] -11.8052907 -8.015226 -16.030452 -9.435097 -7.935097
[2,]  -0.1947093  1.015226   2.030452 -5.564903 -4.064903

$out
[1] -20 -13 -18   2   2   1 -22 -25

$group
[1]  1  2  3  4  5 11 11 12

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

Warning message:
In bxp(list(stats = c(-9, -7.5, -6, -2, -1, -3, -3, -2.5, -1, 1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5lwok1243892189.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]
[1,] -18.0 -7.0 -14.0 -8.0   -4 -25.0  -20
[2,] -13.0 -6.0 -11.0 -4.5   -2 -15.0  -20
[3,] -10.5 -3.5  -8.5 -2.0   -1 -10.0  -20
[4,]  -7.5 -2.0  -5.0 -1.5    1  -6.5  -20
[5,]  -5.0  0.0  -2.0  2.0    2  -3.0  -20

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

$conf
           [,1]      [,2]      [,3]       [,4]       [,5]       [,6] [,7]
[1,] -13.008587 -5.324427 -11.23664 -3.3683201 -2.3683201 -13.876907  -20
[2,]  -7.991413 -1.675573  -5.76336 -0.6316799  0.3683201  -6.123093  -20

$out
[1] -8

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(-18, -13, -10.5, -7.5, -5, -7, -6, -3.5, -2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/63tfc1243892189.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,] 4.262237  9.0  1.750
[2,] 4.824857 12.0  2.375
[3,] 5.681388 13.5  4.250
[4,] 7.053902 19.0  9.375
[5,] 7.949843 23.0 11.000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4.664706 10.30725 1.057253
[2,] 6.698070 16.69275 7.442747

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(4.26223728418147, 4.82485745516368, 5.68138783512352,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1xrwk1243892189.ps tmp/1xrwk1243892189.png")
> system("convert tmp/2se3b1243892189.ps tmp/2se3b1243892189.png")
> system("convert tmp/39ofm1243892189.ps tmp/39ofm1243892189.png")
> system("convert tmp/4muxl1243892189.ps tmp/4muxl1243892189.png")
> system("convert tmp/5lwok1243892189.ps tmp/5lwok1243892189.png")
> system("convert tmp/63tfc1243892189.ps tmp/63tfc1243892189.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.041   0.794   1.298