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.

  Natural language support but running in an English locale

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(24.65,25.24,25.56,25.9,25.87,25.78,25.78,25.74,25.78,25.73,24.67,24.31,24.56,25,25.38,25.99,26.22,26.19,26.22,26.22,26.61,26.72,25.46,25.48,25.59,25.88,26,26.97,27.2,27.19,27.19,27.19,27.26,26.9,26.11,25.87,26.02,26.31,26.37,26.52,26.86,26.92,26.98,26.98,27.03,26.75,26.39,26.3,26.3,26.52,26.53,26.98,27.22,27.34,27.41,27.47,27.46,27.53,27.21,26.91,26.95,26.91,27.39,27.62,27.79,27.88,27.9,28.09,28.46,28.73,27.93,27.61)
> 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,] 24.65 24.56 25.59 26.02 26.30 26.95   NA
 [2,] 25.24 25.00 25.88 26.31 26.52 26.91   NA
 [3,] 25.56 25.38 26.00 26.37 26.53 27.39   NA
 [4,] 25.90 25.99 26.97 26.52 26.98 27.62   NA
 [5,] 25.87 26.22 27.20 26.86 27.22 27.79   NA
 [6,] 25.78 26.19 27.19 26.92 27.34 27.88   NA
 [7,] 25.78 26.22 27.19 26.98 27.41 27.90   NA
 [8,] 25.74 26.22 27.19 26.98 27.47 28.09   NA
 [9,] 25.78 26.61 27.26 27.03 27.46 28.46   NA
[10,] 25.73 26.72 26.90 26.75 27.53 28.73   NA
[11,] 24.67 25.46 26.11 26.39 27.21 27.93   NA
[12,] 24.31 25.48 25.87 26.30 26.91 27.61   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/1j6tg1210618504.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/2new91210618504.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/304is1210618504.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/404g91210618504.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,] 24.560 25.000 25.380 25.900 25.87 25.780 25.780 25.740 25.780 25.730 24.67
[2,] 24.650 25.240 25.560 25.990 26.22 26.190 26.220 26.220 26.610 26.720 25.46
[3,] 25.805 26.095 26.185 26.745 27.03 27.055 27.085 27.085 27.145 26.825 26.25
[4,] 26.300 26.520 26.530 26.980 27.22 27.340 27.410 27.470 27.460 27.530 27.21
[5,] 26.950 26.910 27.390 27.620 27.79 27.880 27.900 28.090 28.460 28.730 27.93
      [,12]
[1,] 24.310
[2,] 25.480
[3,] 26.085
[4,] 26.910
[5,] 27.610

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 24.74070 25.26936 25.55932 26.10642 26.38497 26.31321 26.31741 26.27871
[2,] 26.86930 26.92064 26.81068 27.38358 27.67503 27.79679 27.85259 27.89129
         [,9]    [,10]    [,11]    [,12]
[1,] 26.59672 26.30252 25.12119 25.16260
[2,] 27.69328 27.34748 27.37881 27.00740

$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(24.56, 24.65, 25.805, 26.3, 26.95, 25, 25.24,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/573wd1210618504.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,] 24.310 24.56 25.590 26.020 26.300 26.910   NA
[2,] 24.955 25.42 25.940 26.340 26.720 27.500   NA
[3,] 25.735 26.09 26.935 26.635 27.215 27.835   NA
[4,] 25.780 26.22 27.190 26.950 27.435 28.010   NA
[5,] 25.900 26.72 27.260 27.030 27.530 28.730   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 25.35871 25.72511 26.36487 26.35677 26.88888 27.60239   NA
[2,] 26.11129 26.45489 27.50513 26.91323 27.54112 28.06761   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(24.31, 24.955, 25.735, 25.78, 25.9, 24.56,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6adbn1210618504.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.6581996 1.720 0.64500
[2,] 0.7388202 1.965 0.82750
[3,] 0.8173198 2.235 0.93750
[4,] 0.9717709 2.840 1.12375
[5,] 1.1740145 3.300 1.38250

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] 0.7110694 1.835907 0.8023784
[2,] 0.9235701 2.634093 1.0726216

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.65819956446861, 0.738820224402931, 0.81731976267018,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1j6tg1210618504.ps tmp/1j6tg1210618504.png")
> system("convert tmp/2new91210618504.ps tmp/2new91210618504.png")
> system("convert tmp/304is1210618504.ps tmp/304is1210618504.png")
> system("convert tmp/404g91210618504.ps tmp/404g91210618504.png")
> system("convert tmp/573wd1210618504.ps tmp/573wd1210618504.png")
> system("convert tmp/6adbn1210618504.ps tmp/6adbn1210618504.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.369   1.356   2.584