R version 2.7.1 (2008-06-23)
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(16.8,16.91,16.91,17.16,17.02,17.23,17.22,17.29,17.3,17.22,17.19,17.23,17.36,17.39,17.29,17.28,17.4,17.51,17.54,17.64,17.65,17.5,17.37,17.56,17.49,17.61,17.79,17.83,17.56,17.95,18.09,18.38,18.38,18.44,18.84,19.01,19.06,19.06,18.97,18.98,19.41,19.55,19.64,19.71,19.48,19.48,19.41,19.25,19.14,19.21,19.3,19.53,19.14,19.16,19.24,19.38,19.27,19.27,19.07,19.15,19.24,19.36,19.57,19.59,19.36,19.46,19.65,19.46,19.51,19.64,19.64,19.69,19.28,19.67,19.65,19.6,19.53,19.64,19.67,19.81,19.73,19.87,19.97,20.12,19.94,20.31,20.13,20.22,20.38,20.44,20.34,20.14,19.97,19.82,19.98,20.12)
> 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] 96
> (np <- floor(n / par1))
[1] 8
> 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] 8 8 8 8 8 8 8 8 8 8 8 8
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8] [,9]
 [1,] 16.80 17.36 17.49 19.06 19.14 19.24 19.28 19.94   NA
 [2,] 16.91 17.39 17.61 19.06 19.21 19.36 19.67 20.31   NA
 [3,] 16.91 17.29 17.79 18.97 19.30 19.57 19.65 20.13   NA
 [4,] 17.16 17.28 17.83 18.98 19.53 19.59 19.60 20.22   NA
 [5,] 17.02 17.40 17.56 19.41 19.14 19.36 19.53 20.38   NA
 [6,] 17.23 17.51 17.95 19.55 19.16 19.46 19.64 20.44   NA
 [7,] 17.22 17.54 18.09 19.64 19.24 19.65 19.67 20.34   NA
 [8,] 17.29 17.64 18.38 19.71 19.38 19.46 19.81 20.14   NA
 [9,] 17.30 17.65 18.38 19.48 19.27 19.51 19.73 19.97   NA
[10,] 17.22 17.50 18.44 19.48 19.27 19.64 19.87 19.82   NA
[11,] 17.19 17.37 18.84 19.41 19.07 19.64 19.97 19.98   NA
[12,] 17.23 17.56 19.01 19.25 19.15 19.69 20.12 20.12   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/1pp7g1218631705.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/2jux11218631705.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/34y3j1218631705.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/4fy891218631705.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,] 16.800 16.910 16.910 17.160 17.02 17.230 17.220 17.29 17.300 17.220 17.190
[2,] 17.425 17.500 17.540 17.555 17.48 17.730 17.815 18.01 18.015 17.970 18.105
[3,] 19.100 19.135 19.135 19.255 19.25 19.310 19.440 19.42 19.375 19.375 19.240
[4,] 19.260 19.515 19.610 19.595 19.47 19.595 19.660 19.76 19.620 19.730 19.805
[5,] 19.940 20.310 20.130 20.220 20.38 20.440 20.340 20.14 19.970 19.870 19.980
      [,12]
[1,] 17.230
[2,] 18.285
[3,] 19.200
[4,] 19.905
[5,] 20.120

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 18.07494 18.00939 17.97867 18.11543 18.13836 18.26818 18.40936 18.44242
[2,] 20.12506 20.26061 20.29133 20.39457 20.36164 20.35182 20.47064 20.39758
         [,9]    [,10]    [,11]    [,12]
[1,] 18.47842 18.39184 18.29036 18.29504
[2,] 20.27158 20.35816 20.18964 20.10496

$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(16.8, 17.425, 19.1, 19.26, 19.94, 16.91, 17.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5r3du1218631706.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]   [,8] [,9]
[1,] 16.800 17.280 17.49 18.970 19.070 19.24 19.53 19.820   NA
[2,] 16.965 17.365 17.70 19.060 19.145 19.41 19.62 19.975   NA
[3,] 17.205 17.450 18.02 19.410 19.225 19.54 19.67 20.135   NA
[4,] 17.230 17.550 18.41 19.515 19.285 19.64 19.84 20.325   NA
[5,] 17.300 17.650 19.01 19.710 19.380 19.69 20.12 20.440   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 17.08413 17.36562 17.69616 19.20247 19.16115 19.43510 19.56966 19.97536
[2,] 17.32587 17.53438 18.34384 19.61753 19.28885 19.64490 19.77034 20.29464
     [,9]
[1,]   NA
[2,]   NA

$out
[1] 19.53 19.28

$group
[1] 5 7

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

Warning message:
In bxp(list(stats = c(16.8, 16.965, 17.205, 17.23, 17.3, 17.28,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ukwl1218631706.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.005065 2.650 1.15000
[2,] 1.077111 2.820 1.42375
[3,] 1.145370 3.090 1.71750
[4,] 1.193342 3.215 1.89125
[5,] 1.222749 3.400 1.92500

$n
[1] 12 12 12

$conf
         [,1]     [,2]    [,3]
[1,] 1.092356 2.909838 1.50427
[2,] 1.198384 3.270162 1.93073

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.00506485220749, 1.07711065791414, 1.14537008858105,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1pp7g1218631705.ps tmp/1pp7g1218631705.png")
> system("convert tmp/2jux11218631705.ps tmp/2jux11218631705.png")
> system("convert tmp/34y3j1218631705.ps tmp/34y3j1218631705.png")
> system("convert tmp/4fy891218631705.ps tmp/4fy891218631705.png")
> system("convert tmp/5r3du1218631706.ps tmp/5r3du1218631706.png")
> system("convert tmp/6ukwl1218631706.ps tmp/6ukwl1218631706.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.257   0.784   3.485