R version 2.8.0 (2008-10-20)
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(103.34,103.38,103.64,104.04,104.11,104.11,104.11,104.17,105.16,105.86,106.11,106.11,106.11,106.13,106.67,106.85,106.97,107.02,107.02,107.07,107.76,108.10,108.18,108.22,108.22,108.17,108.31,108.31,108.36,108.46,108.46,108.46,109.43,109.55,109.62,109.70,109.70,109.56,109.92,109.81,109.78,109.80,109.80,109.79,110.40,110.95,111.07,111.09,111.10,111.01,111.01,111.35,111.42,111.24,111.24,111.47,111.57,111.96,112.02,112.02)
> 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,] 103.34 106.11 108.22 109.70 111.10   NA
 [2,] 103.38 106.13 108.17 109.56 111.01   NA
 [3,] 103.64 106.67 108.31 109.92 111.01   NA
 [4,] 104.04 106.85 108.31 109.81 111.35   NA
 [5,] 104.11 106.97 108.36 109.78 111.42   NA
 [6,] 104.11 107.02 108.46 109.80 111.24   NA
 [7,] 104.11 107.02 108.46 109.80 111.24   NA
 [8,] 104.17 107.07 108.46 109.79 111.47   NA
 [9,] 105.16 107.76 109.43 110.40 111.57   NA
[10,] 105.86 108.10 109.55 110.95 111.96   NA
[11,] 106.11 108.18 109.62 111.07 112.02   NA
[12,] 106.11 108.22 109.70 111.09 112.02   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/1z6ea1257374946.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/2yos91257374946.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/3ypr11257374946.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/4ekvl1257374946.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,] 103.34 103.38 103.64 104.04 104.11 104.11 104.11 104.17 105.16 105.86
[2,] 106.11 106.13 106.67 106.85 106.97 107.02 107.02 107.07 107.76 108.10
[3,] 108.22 108.17 108.31 108.31 108.36 108.46 108.46 108.46 109.43 109.55
[4,] 109.70 109.56 109.92 109.81 109.78 109.80 109.80 109.79 110.40 110.95
[5,] 111.10 111.01 111.01 111.35 111.42 111.24 111.24 111.47 111.57 111.96
      [,11]  [,12]
[1,] 106.11 106.11
[2,] 108.18 108.22
[3,] 109.62 109.70
[4,] 111.07 111.09
[5,] 112.02 112.02

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 105.6833 105.7464 106.0136 106.2185 106.3745 106.4957 106.4957 106.5381
[2,] 110.7567 110.5936 110.6064 110.4015 110.3455 110.4243 110.4243 110.3819
         [,9]    [,10]    [,11]    [,12]
[1,] 107.5646 107.5362 107.5779 107.6721
[2,] 111.2954 111.5638 111.6621 111.7279

$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(103.34, 106.11, 108.22, 109.7, 111.1, 103.38,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5l08h1257374946.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,] 103.34 106.11 108.17 109.560 111.010   NA
[2,] 103.84 106.76 108.31 109.785 111.170   NA
[3,] 104.11 107.02 108.46 109.805 111.385   NA
[4,] 105.51 107.93 109.49 110.675 111.765   NA
[5,] 106.11 108.22 109.70 111.090 112.020   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 103.3483 106.4864 107.9218 109.3991 111.1136   NA
[2,] 104.8717 107.5536 108.9982 110.2109 111.6564   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(103.34, 103.84, 104.11, 105.51, 106.11, 106.11,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6o5t81257374946.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,] 2.345965 5.910 2.640
[2,] 2.450929 6.255 2.780
[3,] 2.756027 7.215 2.860
[4,] 2.853463 7.340 3.105
[5,] 3.056465 7.760 3.590

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2.572429 6.720124 2.711765
[2,] 2.939625 7.709876 3.008235

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2.34596462036408, 2.45092947130616, 2.75602682744962,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1z6ea1257374946.ps tmp/1z6ea1257374946.png")
> system("convert tmp/2yos91257374946.ps tmp/2yos91257374946.png")
> system("convert tmp/3ypr11257374946.ps tmp/3ypr11257374946.png")
> system("convert tmp/4ekvl1257374946.ps tmp/4ekvl1257374946.png")
> system("convert tmp/5l08h1257374946.ps tmp/5l08h1257374946.png")
> system("convert tmp/6o5t81257374946.ps tmp/6o5t81257374946.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.677   1.302   1.915