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(57.7,63.6,78,77.4,74.1,85.9,82,78.4,68.1,70.9,85.2,149.6,57.9,63.7,85,66.1,80.2,83.4,85.7,81.8,69.4,76.4,90.3,157.3,65.3,68.4,72.7,86.6,82.6,84.8,93.4,82.2,75.2,83.9,85.4,166.3,70.4,73.9,82.4,92.3,82.7,95.8,105.8,84.2,82.7,88.4,90.2,176.6,69.5,77.3,98.6,86.4,90.8,101.5,112.2,93.6,93.8,90.8,98.1,187.6,75,83.7,99.7,104.9,98.9,117.3,115.7,102.2,101.9,96.6,110,203.7,82.3,93.3,121.9,100.9,107.7,130,123.2,116.1,105.3,107.7,123.9,205.2,90.3,106.9,122.4,111.3,122.6,124.8,139.5,118.8,111,121.2,120.6,219.1,101.3,105,113.4,133.6,123.9,136.2,151.7,121.9,120.2,132.2,125.2,233.8)
> 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] 108
> (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] 9 9 9 9 9 9 9 9 9 9 9 9
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10]
 [1,]  57.7  57.9  65.3  70.4  69.5  75.0  82.3  90.3 101.3    NA
 [2,]  63.6  63.7  68.4  73.9  77.3  83.7  93.3 106.9 105.0    NA
 [3,]  78.0  85.0  72.7  82.4  98.6  99.7 121.9 122.4 113.4    NA
 [4,]  77.4  66.1  86.6  92.3  86.4 104.9 100.9 111.3 133.6    NA
 [5,]  74.1  80.2  82.6  82.7  90.8  98.9 107.7 122.6 123.9    NA
 [6,]  85.9  83.4  84.8  95.8 101.5 117.3 130.0 124.8 136.2    NA
 [7,]  82.0  85.7  93.4 105.8 112.2 115.7 123.2 139.5 151.7    NA
 [8,]  78.4  81.8  82.2  84.2  93.6 102.2 116.1 118.8 121.9    NA
 [9,]  68.1  69.4  75.2  82.7  93.8 101.9 105.3 111.0 120.2    NA
[10,]  70.9  76.4  83.9  88.4  90.8  96.6 107.7 121.2 132.2    NA
[11,]  85.2  90.3  85.4  90.2  98.1 110.0 123.9 120.6 125.2    NA
[12,] 149.6 157.3 166.3 176.6 187.6 203.7 205.2 219.1 233.8    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/1xzhn1291465609.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/html/rcomp/tmp/2q9g81291465609.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/html/rcomp/tmp/3iift1291465609.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/html/rcomp/tmp/4iift1291465609.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] [,11] [,12]
[1,]  57.7  63.6  72.7  66.1  74.1  83.4  82.0  78.4  68.1  70.9  85.2 149.6
[2,]  65.3  68.4  82.4  86.4  82.6  85.9  93.4  82.2  75.2  83.9  90.2 166.3
[3,]  70.4  77.3  98.6  92.3  90.8 101.5 112.2  93.6  93.8  90.8  98.1 187.6
[4,]  82.3  93.3 113.4 104.9 107.7 124.8 123.2 116.1 105.3 107.7 120.6 205.2
[5,] 101.3 106.9 122.4 111.3 123.9 136.2 151.7 121.9 120.2 132.2 125.2 233.8

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

$conf
         [,1]   [,2]      [,3]      [,4]      [,5]      [,6]      [,7]    [,8]
[1,] 61.44667 64.186  82.27333  82.55667  77.58067  81.01267  96.50533  75.746
[2,] 79.35333 90.414 114.92667 102.04333 104.01933 121.98733 127.89467 111.454
          [,9]     [,10]     [,11]    [,12]
[1,]  77.94733  78.26533  82.08933 167.1127
[2,] 109.65267 103.33467 114.11067 208.0873

$out
[1] 133.6

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(57.7, 65.3, 70.4, 82.3, 101.3, 63.6, 68.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5iift1291465609.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] [,10]
[1,] 57.7 57.90 65.30  70.40  77.30  83.70  82.30 106.90 101.30    NA
[2,] 69.5 67.75 73.95  82.55  88.60  97.75 103.10 111.15 116.80    NA
[3,] 77.7 81.00 83.25  86.30  93.70 102.05 111.90 120.90 124.55    NA
[4,] 83.6 85.35 86.00  94.05 100.05 112.85 123.55 123.70 134.90    NA
[5,] 85.9 90.30 93.40 105.80 112.20 117.30 130.00 139.50 151.70    NA

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]      [,6]     [,7]     [,8]
[1,] 71.2689 72.97252 77.75391 81.05477 88.47758  95.16279 102.5726 115.1759
[2,] 84.1311 89.02748 88.74609 91.54523 98.92242 108.93721 121.2274 126.6241
         [,9] [,10]
[1,] 116.2945    NA
[2,] 132.8055    NA

$out
 [1] 149.6 157.3 166.3 176.6  69.5 187.6  75.0 203.7 205.2  90.3 219.1 233.8

$group
 [1] 1 2 3 4 5 5 6 6 7 8 8 9

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

Warning message:
In bxp(list(stats = c(57.7, 69.5, 77.7, 83.6, 85.9, 57.9, 67.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6iift1291465609.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,] 14.63365 40.00 17.00
[2,] 17.16846 43.55 24.35
[3,] 18.93659 50.95 29.95
[4,] 20.58473 64.40 32.45
[5,] 23.57817 84.20 38.90

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 17.37840 41.44018 26.25554
[2,] 20.49478 60.45982 33.64446

$out
[1] 28.76517

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(14.6336465418265, 17.1684551369228, 18.9365908443413,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1xzhn1291465609.ps tmp/1xzhn1291465609.png",intern=TRUE))
character(0)
> try(system("convert tmp/2q9g81291465609.ps tmp/2q9g81291465609.png",intern=TRUE))
character(0)
> try(system("convert tmp/3iift1291465609.ps tmp/3iift1291465609.png",intern=TRUE))
character(0)
> try(system("convert tmp/4iift1291465609.ps tmp/4iift1291465609.png",intern=TRUE))
character(0)
> try(system("convert tmp/5iift1291465609.ps tmp/5iift1291465609.png",intern=TRUE))
character(0)
> try(system("convert tmp/6iift1291465609.ps tmp/6iift1291465609.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.078   0.843   2.476