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(96.96,93.11,95.62,98.30,96.38,100.82,99.06,94.03,102.07,99.31,98.64,101.82,99.14,97.63,100.06,101.32,101.49,105.43,105.09,99.48,108.53,104.34,106.10,107.35,103.00,104.50,105.17,104.84,106.18,108.86,107.77,102.74,112.63,106.26,108.86,111.38,106.85,107.86,107.94,111.38,111.29,113.72,111.88,109.87,113.72,111.71,114.81,112.05,111.54,110.87,110.87,115.48,111.63,116.24,113.56,106.01,110.45,107.77,108.61,108.19)
> 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,]  96.96  99.14 103.00 106.85 111.54   NA
 [2,]  93.11  97.63 104.50 107.86 110.87   NA
 [3,]  95.62 100.06 105.17 107.94 110.87   NA
 [4,]  98.30 101.32 104.84 111.38 115.48   NA
 [5,]  96.38 101.49 106.18 111.29 111.63   NA
 [6,] 100.82 105.43 108.86 113.72 116.24   NA
 [7,]  99.06 105.09 107.77 111.88 113.56   NA
 [8,]  94.03  99.48 102.74 109.87 106.01   NA
 [9,] 102.07 108.53 112.63 113.72 110.45   NA
[10,]  99.31 104.34 106.26 111.71 107.77   NA
[11,]  98.64 106.10 108.86 114.81 108.61   NA
[12,] 101.82 107.35 111.38 112.05 108.19   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/15o9p1258052592.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/2r5q81258052592.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/3tvek1258052592.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/4au5f1258052592.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]
[1,]  96.96  93.11  95.62  98.30  96.38 100.82  99.06  94.03 108.53  99.31
[2,]  99.14  97.63 100.06 101.32 101.49 105.43 105.09  99.48 108.53 104.34
[3,] 103.00 104.50 105.17 104.84 106.18 108.86 107.77 102.74 110.45 106.26
[4,] 106.85 107.86 107.94 111.38 111.29 113.72 111.88 106.01 112.63 107.77
[5,] 111.54 110.87 110.87 115.48 111.63 116.24 113.56 109.87 113.72 111.71
      [,11]  [,12]
[1,] 106.10 101.82
[2,] 106.10 107.35
[3,] 108.61 108.19
[4,] 108.86 111.38
[5,] 108.86 112.05

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

$conf
          [,1]     [,2]      [,3]      [,4]      [,5]     [,6]     [,7]
[1,]  97.55213  97.2715  99.60201  97.73163  99.25534 103.0023 102.9722
[2,] 108.44787 111.7285 110.73799 111.94837 113.10466 114.7177 112.5678
          [,8]     [,9]    [,10]    [,11]    [,12]
[1,]  98.12592 107.5530 103.8364 106.6598 105.3424
[2,] 107.35408 113.3470 108.6836 110.5602 111.0376

$out
[1] 102.07  98.64 114.81

$group
[1]  9 11 11

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

Warning message:
In bxp(list(stats = c(96.96, 99.14, 103, 106.85, 111.54, 93.11,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5w1h91258052592.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,]  93.110  97.630 102.74 106.850 106.010   NA
[2,]  96.000  99.770 104.67 108.905 108.400   NA
[3,]  98.470 102.915 106.22 111.545 110.870   NA
[4,] 100.065 105.765 108.86 112.885 112.595   NA
[5,] 102.070 108.530 112.63 114.810 116.240   NA

$n
[1] 12 12 12 12 12  0

$conf
          [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,]  96.61593 100.1806 104.3089 109.7297 108.9566   NA
[2,] 100.32407 105.6494 108.1311 113.3603 112.7834   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(93.11, 96, 98.47, 100.065, 102.07, 97.63,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6qey31258052592.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.072870 10.230  2.760
[2,] 5.183501 13.450  4.065
[3,] 5.972253 15.250  7.250
[4,] 6.373058 16.005  9.045
[5,] 7.318526 17.760 10.230

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 5.429688 14.08465 4.978589
[2,] 6.514819 16.41535 9.521411

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(4.07286999547003, 5.1835007150327, 5.97225349848516,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/15o9p1258052592.ps tmp/15o9p1258052592.png")
> system("convert tmp/2r5q81258052592.ps tmp/2r5q81258052592.png")
> system("convert tmp/3tvek1258052592.ps tmp/3tvek1258052592.png")
> system("convert tmp/4au5f1258052592.ps tmp/4au5f1258052592.png")
> system("convert tmp/5w1h91258052592.ps tmp/5w1h91258052592.png")
> system("convert tmp/6qey31258052592.ps tmp/6qey31258052592.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.069   0.804   1.301