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(44.13,44.13,44.17,44.14,44.15,44.14,44.14,44.14,44.19,44.29,44.29,44.29,44.29,44.27,44.26,44.33,44.32,44.34,44.34,44.34,44.37,44.47,44.51,44.51,44.51,44.52,44.7,44.84,44.9,44.95,44.94,44.94,44.91,45.28,45.36,45.34,45.34,45.34,45.44,45.62,45.75,45.77,45.77,45.77,46.09,46.25,46.35,46.34,46.34,46.28,46.59,46.42,46.29,46.29,46.29,46.3,46.52,46.66,46.67,46.72,46.72,46.72,46.76,46.89,47.04,47.02,47.02,47.18,47.22,47.8,47.88,47.91)
> 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,] 44.13 44.29 44.51 45.34 46.34 46.72   NA
 [2,] 44.13 44.27 44.52 45.34 46.28 46.72   NA
 [3,] 44.17 44.26 44.70 45.44 46.59 46.76   NA
 [4,] 44.14 44.33 44.84 45.62 46.42 46.89   NA
 [5,] 44.15 44.32 44.90 45.75 46.29 47.04   NA
 [6,] 44.14 44.34 44.95 45.77 46.29 47.02   NA
 [7,] 44.14 44.34 44.94 45.77 46.29 47.02   NA
 [8,] 44.14 44.34 44.94 45.77 46.30 47.18   NA
 [9,] 44.19 44.37 44.91 46.09 46.52 47.22   NA
[10,] 44.29 44.47 45.28 46.25 46.66 47.80   NA
[11,] 44.29 44.51 45.36 46.35 46.67 47.88   NA
[12,] 44.29 44.51 45.34 46.34 46.72 47.91   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/12nkq1211293854.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/2lq361211293854.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/3whcs1211293854.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/4u0zl1211293854.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,] 44.130 44.13 44.17 44.14 44.150 44.14 44.140 44.140 44.19 44.290 44.290
[2,] 44.290 44.27 44.26 44.33 44.320 44.34 44.340 44.340 44.37 44.470 44.510
[3,] 44.925 44.93 45.07 45.23 45.325 45.36 45.355 45.355 45.50 45.765 45.855
[4,] 46.340 46.28 46.59 46.42 46.290 46.29 46.290 46.300 46.52 46.660 46.670
[5,] 46.720 46.72 46.76 46.89 47.040 47.02 47.020 47.180 47.22 47.800 47.880
     [,12]
[1,] 44.29
[2,] 44.51
[3,] 45.84
[4,] 46.72
[5,] 47.91

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 43.60268 43.63349 43.56707 43.88188 44.05429 44.10219 44.09719 44.09074
[2,] 46.24732 46.22651 46.57293 46.57812 46.59571 46.61781 46.61281 46.61926
         [,9]    [,10]    [,11]    [,12]
[1,] 44.11318 44.35238 44.46173 44.41448
[2,] 46.88682 47.17762 47.24827 47.26552

$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(44.13, 44.29, 44.925, 46.34, 46.72, 44.13,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5qzi61211293854.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,] 44.130 44.260 44.510 45.34 46.280 46.720   NA
[2,] 44.140 44.305 44.770 45.53 46.290 46.825   NA
[3,] 44.145 44.340 44.925 45.77 46.380 47.030   NA
[4,] 44.240 44.420 45.115 46.17 46.625 47.510   NA
[5,] 44.290 44.510 45.360 46.35 46.720 47.910   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6] [,7]
[1,] 44.09939 44.28755 44.76764 45.47809 46.2272 46.71757   NA
[2,] 44.19061 44.39245 45.08236 46.06191 46.5328 47.34243   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(44.13, 44.14, 44.145, 44.24, 44.29, 44.26,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6kgfe1211293854.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.092520 2.590 1.66750
[2,] 1.132058 2.670 1.68375
[3,] 1.145165 2.885 1.75375
[4,] 1.299876 3.275 1.89625
[5,] 1.394983 3.620 1.93250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.068623 2.609055 1.656827
[2,] 1.221708 3.160945 1.850673

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.09252002269981, 1.13205837250856, 1.14516519317426,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/12nkq1211293854.ps tmp/12nkq1211293854.png")
> system("convert tmp/2lq361211293854.ps tmp/2lq361211293854.png")
> system("convert tmp/3whcs1211293854.ps tmp/3whcs1211293854.png")
> system("convert tmp/4u0zl1211293854.ps tmp/4u0zl1211293854.png")
> system("convert tmp/5qzi61211293854.ps tmp/5qzi61211293854.png")
> system("convert tmp/6kgfe1211293854.ps tmp/6kgfe1211293854.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.281   0.806   1.993