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(4.23,4.38,4.43,4.44,4.44,4.44,4.44,4.44,4.45,4.45,4.45,4.45,4.45,4.45,4.45,4.45,4.46,4.46,4.46,4.48,4.58,4.67,4.68,4.68,4.69,4.69,4.69,4.69,4.69,4.69,4.69,4.73,4.78,4.79,4.79,4.8,4.8,4.81,5.16,5.26,5.29,5.29,5.29,5.3,5.3,5.3,5.3,5.3,5.3,5.3,5.3,5.35,5.44,5.47,5.47,5.48,5.48,5.48,5.48,5.48,5.48,5.48,5.5,5.55,5.57,5.58,5.58,5.58,5.59,5.59,5.59,5.55)
> 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,] 4.23 4.45 4.69 4.80 5.30 5.48   NA
 [2,] 4.38 4.45 4.69 4.81 5.30 5.48   NA
 [3,] 4.43 4.45 4.69 5.16 5.30 5.50   NA
 [4,] 4.44 4.45 4.69 5.26 5.35 5.55   NA
 [5,] 4.44 4.46 4.69 5.29 5.44 5.57   NA
 [6,] 4.44 4.46 4.69 5.29 5.47 5.58   NA
 [7,] 4.44 4.46 4.69 5.29 5.47 5.58   NA
 [8,] 4.44 4.48 4.73 5.30 5.48 5.58   NA
 [9,] 4.45 4.58 4.78 5.30 5.48 5.59   NA
[10,] 4.45 4.67 4.79 5.30 5.48 5.59   NA
[11,] 4.45 4.68 4.79 5.30 5.48 5.59   NA
[12,] 4.45 4.68 4.80 5.30 5.48 5.55   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/11f3e1291127301.ps",horizontal=F,onefile=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/21f3e1291127301.ps",horizontal=F,onefile=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/31f3e1291127301.ps",horizontal=F,onefile=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/4u6lh1291127301.ps",horizontal=F,onefile=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] [,12]
[1,] 4.230 4.38 4.430 4.440 4.44 4.44 4.44 4.440 4.45 4.450 4.450  4.45
[2,] 4.450 4.45 4.450 4.450 4.46 4.46 4.46 4.480 4.58 4.670 4.680  4.68
[3,] 4.745 4.75 4.925 4.975 4.99 4.99 4.99 5.015 5.04 5.045 5.045  5.05
[4,] 5.300 5.30 5.300 5.350 5.44 5.47 5.47 5.480 5.48 5.480 5.480  5.48
[5,] 5.480 5.48 5.500 5.550 5.57 5.58 5.58 5.580 5.59 5.590 5.590  5.55

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 4.196723 4.201723 4.376723 4.394471 4.357868 4.338517 4.338517 4.369968
[2,] 5.293277 5.298277 5.473277 5.555529 5.622132 5.641483 5.641483 5.660032
         [,9]    [,10]    [,11]    [,12]
[1,] 4.459471 4.522524 4.528974 4.533974
[2,] 5.620529 5.567476 5.561026 5.566026

$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(4.23, 4.45, 4.745, 5.3, 5.48, 4.38, 4.45,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5u6lh1291127301.ps",horizontal=F,onefile=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,] 4.430 4.450 4.690 5.16 5.300 5.480   NA
[2,] 4.435 4.450 4.690 5.21 5.325 5.525   NA
[3,] 4.440 4.460 4.690 5.29 5.470 5.575   NA
[4,] 4.450 4.625 4.785 5.30 5.480 5.585   NA
[5,] 4.450 4.680 4.800 5.30 5.480 5.590   NA

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

$conf
         [,1]     [,2]    [,3]    [,4]     [,5]     [,6] [,7]
[1,] 4.433158 4.380181 4.64667 5.24895 5.399303 5.547634   NA
[2,] 4.446842 4.539819 4.73333 5.33105 5.540697 5.602366   NA

$out
[1] 4.23 4.38 4.80 4.81

$group
[1] 1 1 4 4

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

Warning message:
In bxp(list(stats = c(4.43, 4.435, 4.44, 4.45, 4.45, 4.45, 4.45,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/65gk21291127301.ps",horizontal=F,onefile=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.4488838 1.070 0.66500
[2,] 0.4647323 1.105 0.72625
[3,] 0.4856436 1.140 0.78000
[4,] 0.5114248 1.140 0.88875
[5,] 0.5180508 1.140 0.90750

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] 0.4643469 1.124036 0.7058827
[2,] 0.5069404 1.155964 0.8541173

$out
[1] 1.25

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(0.448883800851252, 0.464732343493579, 0.485643647495754,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/11f3e1291127301.ps tmp/11f3e1291127301.png",intern=TRUE))
character(0)
> try(system("convert tmp/21f3e1291127301.ps tmp/21f3e1291127301.png",intern=TRUE))
character(0)
> try(system("convert tmp/31f3e1291127301.ps tmp/31f3e1291127301.png",intern=TRUE))
character(0)
> try(system("convert tmp/4u6lh1291127301.ps tmp/4u6lh1291127301.png",intern=TRUE))
character(0)
> try(system("convert tmp/5u6lh1291127301.ps tmp/5u6lh1291127301.png",intern=TRUE))
character(0)
> try(system("convert tmp/65gk21291127301.ps tmp/65gk21291127301.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.066   0.829   2.994