R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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(89.3,88.1,93.6,79.7,83.8,62.3,62.3,77.6,80.3,97,94,75.1,74,77.6,75.1,85,75.4,63.2,64.7,77,82.6,97.6,99,75.3,71.6,76.8,83.9,79.7,77.5,73.1,65.6,85.2,98.3,98,100.6,84.1,76.7,82.4,95.5,79.9,82.4,83.6,73.1,91.1,118.6,102.9,111.8,93.9,91.6,92,91.1,97.5,94.7,96.7,78.7,103.5,113.8,106.1,120.3,114.2,106.3,98.8,113.1,97.7,116.3,107.2,94.5,123.5,126.6,126.5,141.4,124.3,124.9,108.9,126.7,107.7,121.8,118.3,122.8,149.5,147,139.3,162.1,142.2,141.4,124.7,114,126.6,121.9,125.1,122.1,135.9,148.4,137.5,145.3,139.9,128.2,115.4,124.7,111.5,121.1,122.5,127.4,143.7,157.8,148.8,162.9,153.9)
> 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,] 89.3 74.0  71.6  76.7  91.6 106.3 124.9 141.4 128.2    NA
 [2,] 88.1 77.6  76.8  82.4  92.0  98.8 108.9 124.7 115.4    NA
 [3,] 93.6 75.1  83.9  95.5  91.1 113.1 126.7 114.0 124.7    NA
 [4,] 79.7 85.0  79.7  79.9  97.5  97.7 107.7 126.6 111.5    NA
 [5,] 83.8 75.4  77.5  82.4  94.7 116.3 121.8 121.9 121.1    NA
 [6,] 62.3 63.2  73.1  83.6  96.7 107.2 118.3 125.1 122.5    NA
 [7,] 62.3 64.7  65.6  73.1  78.7  94.5 122.8 122.1 127.4    NA
 [8,] 77.6 77.0  85.2  91.1 103.5 123.5 149.5 135.9 143.7    NA
 [9,] 80.3 82.6  98.3 118.6 113.8 126.6 147.0 148.4 157.8    NA
[10,] 97.0 97.6  98.0 102.9 106.1 126.5 139.3 137.5 148.8    NA
[11,] 94.0 99.0 100.6 111.8 120.3 141.4 162.1 145.3 162.9    NA
[12,] 75.1 75.3  84.1  93.9 114.2 124.3 142.2 139.9 153.9    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/rcomp/tmp/15dd31295180502.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/rcomp/tmp/2h0q31295180502.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/rcomp/tmp/36r2m1295180502.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/rcomp/tmp/47now1295180502.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,]  71.6  76.8  75.1  79.7  75.4  62.3  62.3  77.0  80.3  97.0  94.0  75.1
[2,]  76.7  82.4  91.1  79.9  82.4  73.1  65.6  85.2  98.3  98.0 100.6  84.1
[3,]  91.6  92.0  95.5  97.5  94.7  96.7  78.7 103.5 118.6 106.1 120.3 114.2
[4,] 124.9 108.9 114.0 107.7 121.1 118.3 122.1 135.9 147.0 137.5 145.3 139.9
[5,] 141.4 124.7 126.7 126.6 121.9 125.1 127.4 149.5 157.8 148.8 162.9 153.9

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

$conf
          [,1]      [,2]      [,3]      [,4]    [,5]      [,6]      [,7]
[1,]  66.21467  78.04333  83.43933  82.85867  74.318  72.89467  48.94333
[2,] 116.98533 105.95667 107.56067 112.14133 115.082 120.50533 108.45667
        [,8]      [,9]     [,10]   [,11]   [,12]
[1,]  76.798  92.95133  85.29667  96.758  84.812
[2,] 130.202 144.24867 126.90333 143.842 143.588

$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.6, 76.7, 91.6, 124.9, 141.4, 76.8, 82.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5gy7l1295180502.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,] 62.30 63.20  65.60  73.10  78.70  94.50 107.70 114.00 111.50    NA
[2,] 76.35 74.55  74.95  81.15  91.80 102.55 120.05 123.40 121.80    NA
[3,] 82.05 76.20  81.80  87.35  97.10 114.70 125.80 131.25 127.80    NA
[4,] 91.45 83.80  91.60  99.20 109.95 125.40 144.60 140.65 151.35    NA
[5,] 97.00 97.60 100.60 118.60 120.30 141.40 162.10 148.40 162.90    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]      [,5]     [,6]     [,7]     [,8]
[1,] 75.16279 71.98101 74.20582 79.11727  88.82166 104.2780 114.6026 123.3822
[2,] 88.93721 80.41899 89.39418 95.58273 105.37834 125.1220 136.9974 139.1178
         [,9] [,10]
[1,] 114.3220    NA
[2,] 141.2780    NA

$out
[1] 99

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(62.3, 76.35, 82.05, 91.45, 97, 63.2, 74.55,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6d0ts1295180502.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,] 16.69597 46.50 22.90
[2,] 19.41281 49.75 33.25
[3,] 25.52078 63.95 44.95
[4,] 27.86615 71.15 49.70
[5,] 30.44220 78.80 56.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 21.66516 54.18932 37.44704
[2,] 29.37641 73.71068 52.45296

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(16.6959659132911, 19.4128080037268, 25.5207802990254,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/15dd31295180502.ps tmp/15dd31295180502.png",intern=TRUE))
character(0)
> try(system("convert tmp/2h0q31295180502.ps tmp/2h0q31295180502.png",intern=TRUE))
character(0)
> try(system("convert tmp/36r2m1295180502.ps tmp/36r2m1295180502.png",intern=TRUE))
character(0)
> try(system("convert tmp/47now1295180502.ps tmp/47now1295180502.png",intern=TRUE))
character(0)
> try(system("convert tmp/5gy7l1295180502.ps tmp/5gy7l1295180502.png",intern=TRUE))
character(0)
> try(system("convert tmp/6d0ts1295180502.ps tmp/6d0ts1295180502.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.170   0.170   1.363