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,8.6,8.8,8.5,8.3,8.2,8,7.9,8,9.3,9.6,9,8.7,8.3,8.4,7.8,7.8,7.6,7.7,7.6,7.6,8.6,8.6,8.2,7.5,7.1,7,6.9,6.6,6.3,6.1,5.9,6,7.2,7.2,6.4,6.1,5.9,6.1,5.9,5.8,5.7,5.6,5.3,5.5,6.5,6.5,6.1,5.9,5.8,6.2,6.5,6.6,6.7,6.6,6.5,6.8,7.8,7.9,7.4)
> 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] 60
> (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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,]  9.0  8.7  7.5  6.1  5.9   NA
 [2,]  8.6  8.3  7.1  5.9  5.8   NA
 [3,]  8.8  8.4  7.0  6.1  6.2   NA
 [4,]  8.5  7.8  6.9  5.9  6.5   NA
 [5,]  8.3  7.8  6.6  5.8  6.6   NA
 [6,]  8.2  7.6  6.3  5.7  6.7   NA
 [7,]  8.0  7.7  6.1  5.6  6.6   NA
 [8,]  7.9  7.6  5.9  5.3  6.5   NA
 [9,]  8.0  7.6  6.0  5.5  6.8   NA
[10,]  9.3  8.6  7.2  6.5  7.8   NA
[11,]  9.6  8.6  7.2  6.5  7.9   NA
[12,]  9.0  8.2  6.4  6.1  7.4   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/15zy81259775186.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/2d2wl1259775186.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/32bor1259775186.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/4rs3u1259775186.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,]  5.9  5.8  6.1  5.9  5.8  5.7  5.6  5.3  5.5   6.5   6.5   6.1
[2,]  6.1  5.9  6.2  6.5  6.6  6.3  6.1  5.9  6.0   7.2   7.2   6.4
[3,]  7.5  7.1  7.0  6.9  6.6  6.7  6.6  6.5  6.8   7.8   7.9   7.4
[4,]  8.7  8.3  8.4  7.8  7.8  7.6  7.7  7.6  7.6   8.6   8.6   8.2
[5,]  9.0  8.6  8.8  8.5  8.3  8.2  8.0  7.9  8.0   9.3   9.6   9.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 5.662847 5.404166 5.445486 5.981423 5.752083 5.781423 5.469444 5.298784
[2,] 9.337153 8.795834 8.554514 7.818577 7.447917 7.618577 7.730556 7.701216
         [,9]    [,10]    [,11]    [,12]
[1,] 5.669444 6.810764 6.910764 6.128125
[2,] 7.930556 8.789236 8.889236 8.671875

$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(5.9, 6.1, 7.5, 8.7, 9, 5.8, 5.9, 7.1, 8.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5uzi51259775186.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,] 7.90 7.60 5.90 5.30 5.80   NA
[2,] 8.10 7.65 6.20 5.65 6.35   NA
[3,] 8.55 8.00 6.75 5.90 6.60   NA
[4,] 9.00 8.50 7.15 6.10 7.10   NA
[5,] 9.60 8.70 7.50 6.50 7.90   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]    [,2]     [,3]     [,4]    [,5] [,6]
[1,] 8.139504 7.61231 6.316699 5.694752 6.25792   NA
[2,] 8.960496 8.38769 7.183301 6.105248 6.94208   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(7.9, 8.1, 8.55, 9, 9.6, 7.6, 7.65, 8, 8.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6emm71259775186.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,] 1.002497 2.40 1.20
[2,] 1.031253 2.50 1.35
[3,] 1.105665 2.65 1.60
[4,] 1.229125 2.85 2.00
[5,] 1.431084 3.10 2.60

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.015415 2.490363 1.303531
[2,] 1.195916 2.809637 1.896469

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.00249688278817, 1.03125312882614, 1.10566534669478,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/15zy81259775186.ps tmp/15zy81259775186.png")
> system("convert tmp/2d2wl1259775186.ps tmp/2d2wl1259775186.png")
> system("convert tmp/32bor1259775186.ps tmp/32bor1259775186.png")
> system("convert tmp/4rs3u1259775186.ps tmp/4rs3u1259775186.png")
> system("convert tmp/5uzi51259775186.ps tmp/5uzi51259775186.png")
> system("convert tmp/6emm71259775186.ps tmp/6emm71259775186.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.067   0.831   1.307