R version 2.8.0 (2008-10-20)
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(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/freestat/rcomp/tmp/1p9ej1292349837.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/freestat/rcomp/tmp/2i0e41292349837.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/freestat/rcomp/tmp/3i0e41292349837.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/freestat/rcomp/tmp/4b9vp1292349837.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/freestat/rcomp/tmp/5b9vp1292349837.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/freestat/rcomp/tmp/6b9vp1292349837.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/1p9ej1292349837.ps tmp/1p9ej1292349837.png",intern=TRUE))
character(0)
> try(system("convert tmp/2i0e41292349837.ps tmp/2i0e41292349837.png",intern=TRUE))
character(0)
> try(system("convert tmp/3i0e41292349837.ps tmp/3i0e41292349837.png",intern=TRUE))
character(0)
> try(system("convert tmp/4b9vp1292349837.ps tmp/4b9vp1292349837.png",intern=TRUE))
character(0)
> try(system("convert tmp/5b9vp1292349837.ps tmp/5b9vp1292349837.png",intern=TRUE))
character(0)
> try(system("convert tmp/6b9vp1292349837.ps tmp/6b9vp1292349837.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.683   1.246   2.109