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(100.00,98.86,96.87,103.18,104.66,103.74,103.87,98.10,98.91,109.43,104.89,88.63,97.03,93.79,97.61,98.17,96.31,97.82,87.52,85.71,89.15,90.65,92.74,80.00,84.33,93.26,115.89,105.21,94.23,106.12,94.41,101.30,98.43,107.94,102.82,92.33,97.48,88.37,92.50,81.48,92.55,98.96,77.28,91.01,81.21,88.41,96.02,91.59,86.30,86.37,105.41,76.73,93.26,95.02,84.32,93.24,89.81,111.66,103.23,103.21)
> 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,] 100.00 97.03  84.33 97.48  86.30   NA
 [2,]  98.86 93.79  93.26 88.37  86.37   NA
 [3,]  96.87 97.61 115.89 92.50 105.41   NA
 [4,] 103.18 98.17 105.21 81.48  76.73   NA
 [5,] 104.66 96.31  94.23 92.55  93.26   NA
 [6,] 103.74 97.82 106.12 98.96  95.02   NA
 [7,] 103.87 87.52  94.41 77.28  84.32   NA
 [8,]  98.10 85.71 101.30 91.01  93.24   NA
 [9,]  98.91 89.15  98.43 81.21  89.81   NA
[10,] 109.43 90.65 107.94 88.41 111.66   NA
[11,] 104.89 92.74 102.82 96.02 103.23   NA
[12,]  88.63 80.00  92.33 91.59 103.21   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/1sld61257951689.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/2ov9g1257951689.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/3mpn31257951689.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/4qx7b1257951689.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]
[1,]  84.33 86.37  92.50  76.73 92.55  95.02  77.28  85.71 81.21  88.41  92.74
[2,]  86.30 88.37  96.87  81.48 93.26  97.82  84.32  91.01 89.15  90.65  96.02
[3,]  97.03 93.26  97.61  98.17 94.23  98.96  87.52  93.24 89.81 107.94 102.82
[4,]  97.48 93.79 105.41 103.18 96.31 103.74  94.41  98.10 98.43 109.43 103.23
[5,] 100.00 98.86 115.89 105.21 96.31 106.12 103.87 101.30 98.91 111.66 104.89
     [,12]
[1,] 88.63
[2,] 88.63
[3,] 91.59
[4,] 92.33
[5,] 92.33

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

$conf
          [,1]     [,2]      [,3]      [,4]     [,5]      [,6]     [,7]
[1,]  89.13024 89.43024  91.57566  82.83683 92.07488  94.77694 80.39043
[2,] 104.92976 97.08976 103.64434 113.50317 96.38512 103.14306 94.64957
         [,8]     [,9]    [,10]     [,11]    [,12]
[1,] 88.23022 83.25278  94.6701  97.72543 88.97559
[2,] 98.24978 96.36722 121.2099 107.91457 94.20441

$out
[1] 104.66  80.00 103.21

$group
[1]  5 12 12

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

Warning message:
In bxp(list(stats = c(84.33, 86.3, 97.03, 97.48, 100, 86.37, 88.37,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5or8k1257951689.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.870 80.000  84.330 77.280  76.730   NA
[2,]  98.480 88.335  93.745 84.925  86.335   NA
[3,] 101.590 93.265  99.865 91.300  93.250   NA
[4,] 104.265 97.320 105.665 94.285 103.220   NA
[5,] 109.430 98.170 115.890 98.960 111.660   NA

$n
[1] 12 12 12 12 12  0

$conf
          [,1]     [,2]     [,3]     [,4]      [,5] [,6]
[1,]  98.95142 89.16688  94.4282 87.03084  85.54864   NA
[2,] 104.22858 97.36312 105.3018 95.56916 100.95136   NA

$out
[1] 88.63

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(96.87, 98.48, 101.59, 104.265, 109.43, 80,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6oxb61257951689.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.514568 11.100  3.050
[2,]  5.099539 12.320  5.670
[3,]  7.267042 16.685  7.875
[4,]  9.674803 23.320 10.635
[5,] 13.008249 28.480 11.180

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 5.180233 11.66783  5.61043
[2,] 9.353850 21.70217 10.13957

$out
[1] 21.70 18.78

$group
[1] 3 3

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

Warning message:
In bxp(list(stats = c(4.51456753189052, 5.0995393457279, 7.26704164631528,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1sld61257951689.ps tmp/1sld61257951689.png")
> system("convert tmp/2ov9g1257951689.ps tmp/2ov9g1257951689.png")
> system("convert tmp/3mpn31257951689.ps tmp/3mpn31257951689.png")
> system("convert tmp/4qx7b1257951689.ps tmp/4qx7b1257951689.png")
> system("convert tmp/5or8k1257951689.ps tmp/5or8k1257951689.png")
> system("convert tmp/6oxb61257951689.ps tmp/6oxb61257951689.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.051   0.815   1.331