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(217.8,218.79,218.99,219.53,219.55,219.74,219.74,219.74,219.8,219.97,220.07,220.07,220.1,225.8,233.17,233.83,233.63,233.63,233.65,233.8,233.84,233.74,233.88,233.88,233.81,234.68,236.14,236.91,236.87,236.78,236.78,236.9,236.94,236.97,236.96,236.94,236.99,237.24,237.62,237.54,237.41,237.4,237.41,237.28,237.17,237.18,237.18,237.18,236.77,239.23,240.23,240.33,240.33,240.34,240.34,240.27,240.29,240.29,240.29,240.29,240.31,239.95,242.33,242.11,241.53,241.53,241.53,241.41,241.41,241.66,241.8,241.99)
> 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] 72
> (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] 6 6 6 6 6 6 6 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6] [,7]
 [1,] 217.80 220.10 233.81 236.99 236.77 240.31   NA
 [2,] 218.79 225.80 234.68 237.24 239.23 239.95   NA
 [3,] 218.99 233.17 236.14 237.62 240.23 242.33   NA
 [4,] 219.53 233.83 236.91 237.54 240.33 242.11   NA
 [5,] 219.55 233.63 236.87 237.41 240.33 241.53   NA
 [6,] 219.74 233.63 236.78 237.40 240.34 241.53   NA
 [7,] 219.74 233.65 236.78 237.41 240.34 241.53   NA
 [8,] 219.74 233.80 236.90 237.28 240.27 241.41   NA
 [9,] 219.80 233.84 236.94 237.17 240.29 241.41   NA
[10,] 219.97 233.74 236.97 237.18 240.29 241.66   NA
[11,] 220.07 233.88 236.96 237.18 240.29 241.80   NA
[12,] 220.07 233.88 236.94 237.18 240.29 241.99   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/freestat/rcomp/tmp/1d0hg1210681084.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/freestat/rcomp/tmp/2jyxv1210681084.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/freestat/rcomp/tmp/35mmd1210681084.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/freestat/rcomp/tmp/4qctf1210681084.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,] 217.80 218.79 233.17 233.830 233.63 233.63 233.650 233.80 233.840 233.740
[2,] 220.10 225.80 233.17 233.830 233.63 233.63 233.650 233.80 233.840 233.740
[3,] 235.29 235.96 236.88 237.225 237.14 237.09 237.095 237.09 237.055 237.075
[4,] 236.99 239.23 240.23 240.330 240.33 240.34 240.340 240.27 240.290 240.290
[5,] 240.31 239.95 242.33 242.110 241.53 241.53 241.530 241.41 241.410 241.660
      [,11]  [,12]
[1,] 233.88 233.88
[2,] 233.88 233.88
[3,] 237.07 237.06
[4,] 240.29 240.29
[5,] 241.80 241.99

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 224.3954 227.2972 232.3261 233.0323 232.8183 232.7618 232.7797 232.9166
[2,] 246.1846 244.6228 241.4339 241.4177 241.4617 241.4182 241.4103 241.2634
         [,9]    [,10]    [,11]    [,12]
[1,] 232.8945 232.8500 232.9353 232.9253
[2,] 241.2155 241.3000 241.2047 241.1947

$out
 [1] 218.99 219.53 219.55 219.74 219.74 219.74 219.80 219.97 220.07 220.07

$group
 [1]  3  4  5  6  7  8  9 10 11 12

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

Warning message:
In bxp(list(stats = c(217.8, 220.1, 235.29, 236.99, 240.31, 218.79,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5eon51210681084.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,] 218.790 233.170 236.140 236.99 240.23 241.410   NA
[2,] 219.260 233.400 236.460 237.18 240.25 241.410   NA
[3,] 219.740 233.695 236.885 237.26 240.29 241.530   NA
[4,] 219.885 233.835 236.940 237.41 240.33 241.895   NA
[5,] 220.070 233.880 236.970 237.62 240.34 242.330   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 219.4549 233.4966 236.6661 237.1551 240.2535 241.3088   NA
[2,] 220.0251 233.8934 237.1039 237.3649 240.3265 241.7512   NA

$out
[1] 217.80 220.10 225.80 233.81 234.68 236.77 239.23 240.31 239.95

$group
[1] 1 2 2 3 3 5 5 6 6

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

Warning message:
In bxp(list(stats = c(218.79, 219.26, 219.74, 219.885, 220.07, 233.17,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6w0vg1210681084.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,] 7.835458 21.160 4.86250
[2,] 7.870063 21.680 4.92125
[3,] 7.934394 21.790 5.09625
[4,] 8.234924 22.245 5.42625
[5,] 8.494588 22.580 5.66500

$n
[1] 12 12 12

$conf
         [,1]    [,2]     [,3]
[1,] 7.767978 21.5323 4.865916
[2,] 8.100809 22.0477 5.326584

$out
[1]  9.558108 23.340000 13.407500 10.712500

$group
[1] 1 2 3 3

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

Warning message:
In bxp(list(stats = c(7.83545786792323, 7.87006294852085, 7.93439359161429,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1d0hg1210681084.ps tmp/1d0hg1210681084.png")
> system("convert tmp/2jyxv1210681084.ps tmp/2jyxv1210681084.png")
> system("convert tmp/35mmd1210681084.ps tmp/35mmd1210681084.png")
> system("convert tmp/4qctf1210681084.ps tmp/4qctf1210681084.png")
> system("convert tmp/5eon51210681084.ps tmp/5eon51210681084.png")
> system("convert tmp/6w0vg1210681084.ps tmp/6w0vg1210681084.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.793   1.160   1.982