R version 2.12.2 (2011-02-25)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-bit)

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(106.42,106.22,106.32,105.81,105.92,107.54,107.34,107.24,107.74,105.71,105.41,106.22,106.32,106.12,106.22,105.92,105.71,105.71,105.92,105.71,105.41,104.49,101.35,99.72,99.01,97.89,95.86,94.95,95.35,95.15,95.46,95.56,95.05,94.64,93.63,93.12,93.53,97.18,96.27,95.15,97.08,101.95,103.07,103.68,102.87,102.56,103.38,103.27,102.89,102.69,101.54,102.9,101.53,101.96,101.99,101.11,101.75,101.71,104.11,103.57,103.32,103.64,103.68,103.79,103.01,101.54,101.9,103.68,104.62,104.11,105.04,104.83,105.05,104.68,107.32,109.9,109.77,110.69,110.54,110.89,110.95,109.73,110.85,110.39,110.58,110.4,111.07,110.86,111.38,111.44,110.36,110.06,108.34,107.94,107.39,107.1,107.61,107.74,106.9,106.71,106.6,108.21,110.54,110.91,109.51,110.27,111.39,112.13,111.64)
> 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] 109
> (np <- floor(n / par1))
[1] 9
> 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] 10  9  9  9  9  9  9  9  9  9  9  9
> arr
        [,1]   [,2]  [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 106.42 106.32 99.01  93.53 102.89 103.32 105.05 110.58 107.61 111.64
 [2,] 106.22 106.12 97.89  97.18 102.69 103.64 104.68 110.40 107.74     NA
 [3,] 106.32 106.22 95.86  96.27 101.54 103.68 107.32 111.07 106.90     NA
 [4,] 105.81 105.92 94.95  95.15 102.90 103.79 109.90 110.86 106.71     NA
 [5,] 105.92 105.71 95.35  97.08 101.53 103.01 109.77 111.38 106.60     NA
 [6,] 107.54 105.71 95.15 101.95 101.96 101.54 110.69 111.44 108.21     NA
 [7,] 107.34 105.92 95.46 103.07 101.99 101.90 110.54 110.36 110.54     NA
 [8,] 107.24 105.71 95.56 103.68 101.11 103.68 110.89 110.06 110.91     NA
 [9,] 107.74 105.41 95.05 102.87 101.75 104.62 110.95 108.34 109.51     NA
[10,] 105.71 104.49 94.64 102.56 101.71 104.11 109.73 107.94 110.27     NA
[11,] 105.41 101.35 93.63 103.38 104.11 105.04 110.85 107.39 111.39     NA
[12,] 106.22  99.72 93.12 103.27 103.57 104.83 110.39 107.10 112.13     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/wessaorg/rcomp/tmp/150i61304686313.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2nrv81304686313.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3uacr1304686313.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/49cbu1304686313.ps",horizontal=F,onefile=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,]  99.010  97.89  95.86 102.90  95.35  95.15  95.46  95.56  95.05  94.64
[2,] 102.890 102.69 101.54 102.90 101.53 101.95 101.99 103.68 102.87 102.56
[3,] 105.685 104.68 106.22 105.81 105.71 105.71 105.92 105.71 105.41 104.49
[4,] 107.610 106.22 106.90 106.71 106.60 108.21 110.36 110.06 108.34 107.94
[5,] 111.640 110.40 111.07 110.86 111.38 111.44 110.54 110.91 110.95 110.27
      [,11]  [,12]
[1,] 101.35  99.72
[2,] 103.38 103.27
[3,] 105.04 104.83
[4,] 107.39 107.10
[5,] 111.39 112.13

$n
 [1] 10  9  9  9  9  9  9  9  9  9  9  9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 103.3267 102.8209 103.3971 103.8034 103.0398 102.4131 101.5118 102.3499
[2,] 108.0433 106.5391 109.0429 107.8166 108.3802 109.0069 110.3282 109.0701
         [,9]    [,10]    [,11]    [,12]
[1,] 102.5291 101.6565 102.9281 102.8129
[2,] 108.2909 107.3235 107.1519 106.8471

$out
[1] 93.53 97.18 94.95 95.15 93.63 93.12

$group
[1]  1  2  4  4 11 12

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(99.01, 102.89, 105.685, 107.61, 111.64, 97.89,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/5j8uc1304686313.ps",horizontal=F,onefile=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]    [,8]    [,9]
[1,] 105.410 104.49 93.630  93.530 101.110 101.540 107.320 107.100 106.600
[2,] 105.865 104.95 94.795  96.675 101.625 103.165 108.525 108.140 107.255
[3,] 106.270 105.71 95.250 102.255 101.975 103.680 110.145 110.380 108.860
[4,] 107.290 106.02 95.710 103.170 102.895 104.365 110.770 110.965 110.725
[5,] 107.740 106.32 95.860 103.680 104.110 105.040 110.950 111.440 112.130
      [,10]
[1,] 111.64
[2,] 111.64
[3,] 111.64
[4,] 111.64
[5,] 111.64

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

$conf
         [,1]     [,2]     [,3]      [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 105.6200 105.2220 94.83266  99.29259 101.3957 103.1327 109.1210 109.0915
[2,] 106.9200 106.1980 95.66734 105.21741 102.5543 104.2273 111.1690 111.6685
         [,9]  [,10]
[1,] 107.2773 111.64
[2,] 110.4427 111.64

$out
[1] 101.35  99.72  99.01  97.89  93.12 105.05 104.68

$group
[1] 2 2 3 3 3 7 7

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

Warning message:
In bxp(list(stats = c(105.41, 105.865, 106.27, 107.29, 107.74, 104.49,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/60mi61304686313.ps",horizontal=F,onefile=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.780358 13.220 3.530
[2,] 4.970398 15.280 3.920
[3,] 5.186329 15.905 5.215
[4,] 5.369482 17.025 5.865
[5,] 5.675113 19.010 8.370

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 5.004304 15.10909 4.327872
[2,] 5.368354 16.70091 6.102128

$out
[1] 4.337167

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(4.78035825017331, 4.97039817213143, 5.18632928475192,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/150i61304686313.ps tmp/150i61304686313.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nrv81304686313.ps tmp/2nrv81304686313.png",intern=TRUE))
character(0)
> try(system("convert tmp/3uacr1304686313.ps tmp/3uacr1304686313.png",intern=TRUE))
character(0)
> try(system("convert tmp/49cbu1304686313.ps tmp/49cbu1304686313.png",intern=TRUE))
character(0)
> try(system("convert tmp/5j8uc1304686313.ps tmp/5j8uc1304686313.png",intern=TRUE))
character(0)
> try(system("convert tmp/60mi61304686313.ps tmp/60mi61304686313.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.990   0.380   4.435