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(100.00,99.97,101.03,101.00,101.30,101.43,101.49,102.14,102.58,102.34,102.07,101.83,102.14,102.58,102.63,102.74,103.32,103.27,102.48,102.14,102.07,101.69,101.15,100.90,101.77,101.93,102.27,102.49,102.80,102.82,102.83,102.89,102.87,102.67,102.96,103.22,103.53,104.63,104.63,104.17,103.93,104.01,104.16,105.22,105.85,106.21,105.77,105.63,106.49,107.51,110.43,111.42,111.58,111.34,111.08,111.66,112.36,112.31,111.52,110.87,111.13,112.71,113.25,113.09,112.55,112.87,113.59,115.14,116.38,116.50,116.25,116.73)
> 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,] 100.00 102.14 101.77 103.53 106.49 111.13   NA
 [2,]  99.97 102.58 101.93 104.63 107.51 112.71   NA
 [3,] 101.03 102.63 102.27 104.63 110.43 113.25   NA
 [4,] 101.00 102.74 102.49 104.17 111.42 113.09   NA
 [5,] 101.30 103.32 102.80 103.93 111.58 112.55   NA
 [6,] 101.43 103.27 102.82 104.01 111.34 112.87   NA
 [7,] 101.49 102.48 102.83 104.16 111.08 113.59   NA
 [8,] 102.14 102.14 102.89 105.22 111.66 115.14   NA
 [9,] 102.58 102.07 102.87 105.85 112.36 116.38   NA
[10,] 102.34 101.69 102.67 106.21 112.31 116.50   NA
[11,] 102.07 101.15 102.96 105.77 111.52 116.25   NA
[12,] 101.83 100.90 103.22 105.63 110.87 116.73   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/107fu1258048117.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/2xjak1258048117.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/39d1d1258048117.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/4psdj1258048117.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]
[1,] 100.000  99.970 101.03 101.000 101.300 101.43 101.490 102.140 102.07
[2,] 101.770 101.930 102.27 102.490 102.800 102.82 102.480 102.140 102.58
[3,] 102.835 103.605 103.63 103.455 103.625 103.64 103.495 104.055 104.36
[4,] 106.490 107.510 110.43 111.420 111.580 111.34 111.080 111.660 112.36
[5,] 111.130 112.710 113.25 113.090 112.550 112.87 113.590 115.140 116.38
      [,10]   [,11]   [,12]
[1,] 101.69 101.150 100.900
[2,] 102.34 102.070 101.830
[3,] 104.44 104.365 104.425
[4,] 112.31 111.520 110.870
[5,] 116.50 116.250 116.730

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

$conf
          [,1]     [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  99.79045 100.0057  98.36654  97.69486  97.96162  98.14432  97.94772
[2,] 105.87955 107.2043 108.89346 109.21514 109.28838 109.13568 109.04228
         [,8]      [,9]     [,10]     [,11]    [,12]
[1,]  97.9143  98.05158  98.00903  98.26944  98.5939
[2,] 110.1957 110.66842 110.87097 110.46056 110.2561

$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(100, 101.77, 102.835, 106.49, 111.13, 99.97,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5t7mf1258048117.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,]  99.970 100.900 101.77 103.530 110.43 111.130   NA
[2,] 101.015 101.880 102.38 104.085 110.65 112.790   NA
[3,] 101.460 102.310 102.81 104.630 111.38 113.420   NA
[4,] 102.105 102.685 102.88 105.700 111.62 116.315   NA
[5,] 102.580 103.320 103.22 106.210 112.36 116.730   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 100.9628 101.9428 102.5819 103.8934 110.9376 111.8122   NA
[2,] 101.9572 102.6772 103.0381 105.3666 111.8224 115.0278   NA

$out
[1] 106.49 107.51

$group
[1] 5 5

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

Warning message:
In bxp(list(stats = c(99.97, 101.015, 101.46, 102.105, 102.58, 100.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6h3pg1258048117.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,] 4.037403 11.130 6.57500
[2,] 4.856980 11.765 6.59750
[3,] 5.101945 12.480 6.91875
[4,] 6.003470 14.560 7.75625
[5,] 6.138870 15.830 8.36250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4.579023 11.20518 6.390236
[2,] 5.624867 13.75482 7.447264

$out
[1] 3.8875 4.6975

$group
[1] 3 3

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

Warning message:
In bxp(list(stats = c(4.03740345601807, 4.85698010986061, 5.10194487239567,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/107fu1258048117.ps tmp/107fu1258048117.png")
> system("convert tmp/2xjak1258048117.ps tmp/2xjak1258048117.png")
> system("convert tmp/39d1d1258048117.ps tmp/39d1d1258048117.png")
> system("convert tmp/4psdj1258048117.ps tmp/4psdj1258048117.png")
> system("convert tmp/5t7mf1258048117.ps tmp/5t7mf1258048117.png")
> system("convert tmp/6h3pg1258048117.ps tmp/6h3pg1258048117.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.059   0.816   1.315