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(12.23,12.25,12.33,12.37,12.38,12.41,12.41,12.41,12.43,12.44,12.46,12.46,12.46,12.47,12.55,12.58,12.6,12.62,12.62,12.66,12.7,12.76,12.79,12.81,12.86,13.03,13.09,13.15,13.17,13.18,13.18,13.19,13.19,13.19,13.25,13.24,13.3,13.35,13.35,13.35,13.36,13.37,13.38,13.39,13.4,13.4,13.44,13.44,13.44,13.47,13.52,13.58,13.65,13.66,13.68,13.68,13.77,13.78,13.79,13.79,13.82,13.86,13.94,14,14.02,14.03,14.07,14.1,14.1,14.16,14.17,14.17)
> 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,] 12.23 12.46 12.86 13.30 13.44 13.82   NA
 [2,] 12.25 12.47 13.03 13.35 13.47 13.86   NA
 [3,] 12.33 12.55 13.09 13.35 13.52 13.94   NA
 [4,] 12.37 12.58 13.15 13.35 13.58 14.00   NA
 [5,] 12.38 12.60 13.17 13.36 13.65 14.02   NA
 [6,] 12.41 12.62 13.18 13.37 13.66 14.03   NA
 [7,] 12.41 12.62 13.18 13.38 13.68 14.07   NA
 [8,] 12.41 12.66 13.19 13.39 13.68 14.10   NA
 [9,] 12.43 12.70 13.19 13.40 13.77 14.10   NA
[10,] 12.44 12.76 13.19 13.40 13.78 14.16   NA
[11,] 12.46 12.79 13.25 13.44 13.79 14.17   NA
[12,] 12.46 12.81 13.24 13.44 13.79 14.17   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/167kx1210952020.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/2utpi1210952020.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/3tje01210952020.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/4nxaa1210952020.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,] 12.23 12.25 12.33 12.37 12.380 12.410 12.41 12.41 12.430 12.440 12.460
[2,] 12.46 12.47 12.55 12.58 12.600 12.620 12.62 12.66 12.700 12.760 12.790
[3,] 13.08 13.19 13.22 13.25 13.265 13.275 13.28 13.29 13.295 13.295 13.345
[4,] 13.44 13.47 13.52 13.58 13.650 13.660 13.68 13.68 13.770 13.780 13.790
[5,] 13.82 13.86 13.94 14.00 14.020 14.030 14.07 14.10 14.100 14.160 14.170
     [,12]
[1,] 12.46
[2,] 12.81
[3,] 13.34
[4,] 13.79
[5,] 14.17

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 12.44787 12.54497 12.59432 12.60497 12.58772 12.60417 12.59627 12.63207
[2,] 13.71213 13.83503 13.84568 13.89503 13.94228 13.94583 13.96373 13.94793
         [,9]    [,10]    [,11]    [,12]
[1,] 12.60482 12.63707 12.69997 12.70787
[2,] 13.98518 13.95293 13.99003 13.97213

$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(12.23, 12.46, 13.08, 13.44, 13.82, 12.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/53le61210952020.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,] 12.230 12.460 13.03 13.300 13.440 13.82   NA
[2,] 12.350 12.565 13.12 13.350 13.550 13.97   NA
[3,] 12.410 12.620 13.18 13.375 13.670 14.05   NA
[4,] 12.435 12.730 13.19 13.400 13.775 14.13   NA
[5,] 12.460 12.810 13.25 13.440 13.790 14.17   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 12.37123 12.54474 13.14807 13.35219 13.56738 13.97702   NA
[2,] 12.44877 12.69526 13.21193 13.39781 13.77262 14.12298   NA

$out
[1] 12.86

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(12.23, 12.35, 12.41, 12.435, 12.46, 12.46,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6zlpn1210952020.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,] 0.6057392 1.590 0.78500
[2,] 0.6137881 1.615 0.79875
[3,] 0.6243043 1.650 0.82250
[4,] 0.6298108 1.700 0.84000
[5,] 0.6360005 1.720 0.85500

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] 0.6169962 1.611231 0.8036856
[2,] 0.6316124 1.688769 0.8413144

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.60573921781572, 0.613788052778113, 0.624304281539835,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/167kx1210952020.ps tmp/167kx1210952020.png")
> system("convert tmp/2utpi1210952020.ps tmp/2utpi1210952020.png")
> system("convert tmp/3tje01210952020.ps tmp/3tje01210952020.png")
> system("convert tmp/4nxaa1210952020.ps tmp/4nxaa1210952020.png")
> system("convert tmp/53le61210952020.ps tmp/53le61210952020.png")
> system("convert tmp/6zlpn1210952020.ps tmp/6zlpn1210952020.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.270   0.818   3.091