R version 2.6.2 (2008-02-08)
Copyright (C) 2008 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(29.90,28.77,15.64,23.73,25.65,21.81,28.97,24.29,25.33,28.84,19.99,19.75,22.70,23.28,24.15,20.38,27.75,27.31,25.61,22.64,26.05,28.07,21.02,25.00,17.93,35.45,17.70,28.53,26.55,26.51,30.78,26.83,27.49,25.89,20.44,19.79,18.14,27.98,35.90,34.38,21.58,21.53,31.14,28.25,25.16,20.51,30.05,20.17,32.37,22.46,25.40,19.82,18.14,20.10,20.25,19.73,24.74,26.17,20.14,31.71,26.66,20.75,20.01,26.67,23.91,26.81,29.31,31.76,22.99,23.94,27.04,20.28,23.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] 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,] 29.90 22.70 17.93 18.14 32.37 26.66 23.32
 [2,] 28.77 23.28 35.45 27.98 22.46 20.75    NA
 [3,] 15.64 24.15 17.70 35.90 25.40 20.01    NA
 [4,] 23.73 20.38 28.53 34.38 19.82 26.67    NA
 [5,] 25.65 27.75 26.55 21.58 18.14 23.91    NA
 [6,] 21.81 27.31 26.51 21.53 20.10 26.81    NA
 [7,] 28.97 25.61 30.78 31.14 20.25 29.31    NA
 [8,] 24.29 22.64 26.83 28.25 19.73 31.76    NA
 [9,] 25.33 26.05 27.49 25.16 24.74 22.99    NA
[10,] 28.84 28.07 25.89 20.51 26.17 23.94    NA
[11,] 19.99 21.02 20.44 30.05 20.14 27.04    NA
[12,] 19.75 25.00 19.79 20.17 31.71 20.28    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/1ym2e1210622023.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/2zzqj1210622023.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/32rk81210622023.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/42wby1210622023.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,] 17.93 20.75 15.64 19.82 18.14 20.10 20.25 19.73 22.990 20.51 19.99 19.750
[2,] 20.42 22.46 17.70 20.38 21.58 21.53 25.61 22.64 24.740 23.94 20.14 19.790
[3,] 23.32 25.63 22.08 25.20 24.78 24.16 29.14 25.56 25.245 26.03 20.73 20.225
[4,] 28.28 28.77 25.40 28.53 26.55 26.81 30.78 28.25 26.050 28.07 27.04 25.000
[5,] 32.37 35.45 35.90 34.38 27.75 27.31 31.14 31.76 27.490 28.84 30.05 31.710

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 18.62613 21.55985 17.11325 19.94299 21.57419 20.75423 25.80518 21.94137
[2,] 28.01387 29.70015 27.04675 30.45701 27.98581 27.56577 32.47482 29.17863
         [,9]    [,10]    [,11]    [,12]
[1,] 24.40001 23.36602 16.27928 16.86438
[2,] 26.08999 28.69398 25.18072 23.58562

$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(17.93, 20.42, 23.32, 28.28, 32.37, 20.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5m0yp1210622023.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,] 15.640 20.380 17.700 18.140 18.140 20.010 23.32
[2,] 20.900 22.670 20.115 21.020 19.960 21.870 23.32
[3,] 24.810 24.575 26.530 26.570 21.355 25.300 23.32
[4,] 28.805 26.680 28.010 30.595 25.785 26.925 23.32
[5,] 29.900 28.070 35.450 35.900 32.370 31.760 23.32

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]  [,7]
[1,] 21.20448 22.74601 22.92904 22.20278 18.69818 22.99438 23.32
[2,] 28.41552 26.40399 30.13096 30.93722 24.01182 27.60562 23.32

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15.64, 20.9, 24.81, 28.805, 29.9, 20.38, 22.67,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6v7631210622023.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.485162  4.500 1.02500
[2,] 3.382236  8.970 3.94875
[3,] 4.301951 11.425 4.98875
[4,] 5.460227 14.500 6.35875
[5,] 7.274415 20.260 7.86000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 3.354166  8.90273 3.889533
[2,] 5.249737 13.94727 6.087967

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.48516216847409, 3.38223608366654, 4.30195127119066,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ym2e1210622023.ps tmp/1ym2e1210622023.png")
> system("convert tmp/2zzqj1210622023.ps tmp/2zzqj1210622023.png")
> system("convert tmp/32rk81210622023.ps tmp/32rk81210622023.png")
> system("convert tmp/42wby1210622023.ps tmp/42wby1210622023.png")
> system("convert tmp/5m0yp1210622023.ps tmp/5m0yp1210622023.png")
> system("convert tmp/6v7631210622023.ps tmp/6v7631210622023.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.366   0.893   2.220