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(71.42,71.47,71.64,71.55,71.62,71.65,71.65,71.65,71.6,71.7,71.73,71.83,71.83,71.87,71.91,71.99,72.1,72.12,72.12,72.12,72.25,72.59,72.72,72.76,72.76,72.91,73,73.16,73.16,73.11,73.11,73.33,73.51,73.66,73.65,73.65,73.65,73.65,73.71,73.73,73.85,73.77,73.77,73.78,73.88,74.3,74.53,74.71,74.71,74.78,74.9,74.65,74.65,74.53,74.53,74.53,74.65,74.85,74.96,74.96,74.96,75.19,74.98,75.54,75.61,75.59,75.58,75.44,75.37,75.22,75.33,75.33)
> 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,] 71.42 71.83 72.76 73.65 74.71 74.96   NA
 [2,] 71.47 71.87 72.91 73.65 74.78 75.19   NA
 [3,] 71.64 71.91 73.00 73.71 74.90 74.98   NA
 [4,] 71.55 71.99 73.16 73.73 74.65 75.54   NA
 [5,] 71.62 72.10 73.16 73.85 74.65 75.61   NA
 [6,] 71.65 72.12 73.11 73.77 74.53 75.59   NA
 [7,] 71.65 72.12 73.11 73.77 74.53 75.58   NA
 [8,] 71.65 72.12 73.33 73.78 74.53 75.44   NA
 [9,] 71.60 72.25 73.51 73.88 74.65 75.37   NA
[10,] 71.70 72.59 73.66 74.30 74.85 75.22   NA
[11,] 71.73 72.72 73.65 74.53 74.96 75.33   NA
[12,] 71.83 72.76 73.65 74.71 74.96 75.33   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/155jb1210536935.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/20dcs1210536935.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/3gfnb1210536935.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/4mza51210536935.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,] 71.420 71.47 71.640 71.550 71.620 71.65 71.65 71.650 71.600 71.70 71.73
[2,] 71.830 71.87 71.910 71.990 72.100 72.12 72.12 72.120 72.250 72.59 72.72
[3,] 73.205 73.28 73.355 73.445 73.505 73.44 73.44 73.555 73.695 73.98 74.09
[4,] 74.710 74.78 74.900 74.650 74.650 74.53 74.53 74.530 74.650 74.85 74.96
[5,] 74.960 75.19 74.980 75.540 75.610 75.59 75.58 75.440 75.370 75.22 75.33
     [,12]
[1,] 71.83
[2,] 72.76
[3,] 74.18
[4,] 74.96
[5,] 75.33

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 71.3473 71.40296 71.42635 71.72921 71.86017 71.88547 71.88547 72.00047
[2,] 75.0627 75.15704 75.28365 75.16079 75.14983 74.99453 74.99453 75.10953
         [,9]    [,10]    [,11]    [,12]
[1,] 72.14692 72.52223 72.64513 72.76093
[2,] 75.24308 75.43777 75.53487 75.59907

$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(71.42, 71.83, 73.205, 74.71, 74.96, 71.47,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/598hm1210536935.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,] 71.470 71.83 72.760 73.650 74.530 74.960   NA
[2,] 71.575 71.95 73.055 73.720 74.590 75.205   NA
[3,] 71.645 72.12 73.160 73.775 74.680 75.350   NA
[4,] 71.675 72.42 73.580 74.090 74.875 75.560   NA
[5,] 71.730 72.76 73.660 74.530 74.960 75.610   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 71.59939 71.90563 72.92054 73.60624 74.55001 75.18808   NA
[2,] 71.69061 72.33437 73.39946 73.94376 74.80999 75.51192   NA

$out
[1] 71.42 71.83 74.71

$group
[1] 1 1 4

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

Warning message:
In bxp(list(stats = c(71.47, 71.575, 71.645, 71.675, 71.73, 71.83,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6h8gc1210536935.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.358249 3.340 1.8550
[2,] 1.407606 3.530 1.9075
[3,] 1.453347 3.745 1.9725
[4,] 1.496657 3.935 2.2525
[5,] 1.530433 3.990 2.4200

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.412730 3.560277 1.815143
[2,] 1.493964 3.929723 2.129857

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.35824887262975, 1.40760623313809, 1.45334682911849,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/155jb1210536935.ps tmp/155jb1210536935.png")
> system("convert tmp/20dcs1210536935.ps tmp/20dcs1210536935.png")
> system("convert tmp/3gfnb1210536935.ps tmp/3gfnb1210536935.png")
> system("convert tmp/4mza51210536935.ps tmp/4mza51210536935.png")
> system("convert tmp/598hm1210536935.ps tmp/598hm1210536935.png")
> system("convert tmp/6h8gc1210536935.ps tmp/6h8gc1210536935.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.380   1.347   2.593