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(78.1,66.7,79,65.2,66.5,77.2,80.2,77.9,78,86.8,92.9,185.8,91,79.1,84.2,70.1,71.3,79.6,92.3,78.7,82.5,98.2,115.4,205.6,94,83.2,80.3,70.4,71.1,78.8,86.3,77.5,80.1,89.8,99.9,218,85.4,77.5,78.6,68.8,64.8,79.8,94.3,79.9,87.5,99.1,109.9,273.6,91.3,80.6,80.4,71.8,75.5,86.6,91.5,86.8,84.6,88.6,102.1,260.3,79,70.6,79.3,66.8,61.2,72.5,83.5,75.8,83.4,89.4,104.9,251.6,80,76.3,81.1,63.1,63.5,78.8,91.7,83.8,83.8,95.8,108.9,258.2,88.7,79.5,74.3,70.5,59.1,73.2,81.2,75,74.6,89.5,107,246.4,83.6,72.1,68.7,60.1,61.1,72.7,85.3,71.4,75.2,89.8,100.9,222.7)
> 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,]  78.1  91.0  94.0  85.4  91.3  79.0  80.0  88.7  83.6    NA
 [2,]  66.7  79.1  83.2  77.5  80.6  70.6  76.3  79.5  72.1    NA
 [3,]  79.0  84.2  80.3  78.6  80.4  79.3  81.1  74.3  68.7    NA
 [4,]  65.2  70.1  70.4  68.8  71.8  66.8  63.1  70.5  60.1    NA
 [5,]  66.5  71.3  71.1  64.8  75.5  61.2  63.5  59.1  61.1    NA
 [6,]  77.2  79.6  78.8  79.8  86.6  72.5  78.8  73.2  72.7    NA
 [7,]  80.2  92.3  86.3  94.3  91.5  83.5  91.7  81.2  85.3    NA
 [8,]  77.9  78.7  77.5  79.9  86.8  75.8  83.8  75.0  71.4    NA
 [9,]  78.0  82.5  80.1  87.5  84.6  83.4  83.8  74.6  75.2    NA
[10,]  86.8  98.2  89.8  99.1  88.6  89.4  95.8  89.5  89.8    NA
[11,]  92.9 115.4  99.9 109.9 102.1 104.9 108.9 107.0 100.9    NA
[12,] 185.8 205.6 218.0 273.6 260.3 251.6 258.2 246.4 222.7    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/18ey51291657082.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/2j5f91291657082.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/3j5f91291657082.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/4uwet1291657082.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,] 78.1 66.7 78.6 60.1 59.1 72.5 80.2 71.4 74.6  86.8  92.9 185.8
[2,] 80.0 72.1 78.6 65.2 61.2 73.2 83.5 75.8 78.0  89.4 100.9 218.0
[3,] 85.4 77.5 79.3 68.8 64.8 78.8 86.3 77.9 82.5  89.8 104.9 246.4
[4,] 91.0 79.5 80.4 70.4 71.1 79.6 91.7 79.9 83.8  95.8 108.9 258.2
[5,] 94.0 83.2 81.1 71.8 75.5 86.6 94.3 83.8 87.5  99.1 115.4 273.6

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

$conf
         [,1]     [,2]   [,3]     [,4]   [,5]     [,6]     [,7]     [,8]
[1,] 79.60667 73.60267 78.352 66.06133 59.586 75.42933 81.98133 75.74067
[2,] 91.19333 81.39733 80.248 71.53867 70.014 82.17067 90.61867 80.05933
         [,9]    [,10]    [,11]   [,12]
[1,] 79.44533 86.42933 100.6867 225.228
[2,] 85.55467 93.17067 109.1133 267.572

$out
[1] 84.2 74.3 68.7 86.8

$group
[1] 3 3 3 8

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

Warning message:
In bxp(list(stats = c(78.1, 80, 85.4, 91, 94, 66.7, 72.1, 77.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5uwet1291657082.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,] 65.20  70.10 70.40  64.80  71.8  61.20  63.10  59.10  60.10    NA
[2,] 71.95  78.90 78.15  78.05  80.5  71.55  77.55  73.75  70.05    NA
[3,] 78.05  83.35 81.75  82.65  86.7  79.15  82.45  77.25  73.95    NA
[4,] 83.50  95.25 91.90  96.70  91.4  86.45  93.75  89.10  87.55    NA
[5,] 92.90 115.40 99.90 109.90 102.1 104.90 108.90 107.00 100.90    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 72.78197 75.89266 75.47853 74.14361 81.72844 72.35401 75.06107 70.24876
[2,] 83.31803 90.80734 88.02147 91.15639 91.67156 85.94599 89.83893 84.25124
         [,9] [,10]
[1,] 65.96813    NA
[2,] 81.93187    NA

$out
[1] 185.8 205.6 218.0 273.6 260.3 251.6 258.2 246.4 222.7

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

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

Warning message:
In bxp(list(stats = c(65.2, 71.95, 78.05, 83.5, 92.9, 70.1, 78.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6uwet1291657082.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,] 3.940107 11.70  1.80
[2,] 4.491847 13.50  5.50
[3,] 4.916980 15.45  6.90
[4,] 5.719997 16.45  9.05
[5,] 6.595095 16.50 11.00

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4.356812 14.10449 5.280821
[2,] 5.477147 16.79551 8.519179

$out
[1] 29.14048 22.50000 87.80000 40.20000

$group
[1] 1 2 2 3

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

Warning message:
In bxp(list(stats = c(3.94010716154325, 4.49184699649609, 4.91697982765181,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/18ey51291657082.ps tmp/18ey51291657082.png",intern=TRUE))
character(0)
> try(system("convert tmp/2j5f91291657082.ps tmp/2j5f91291657082.png",intern=TRUE))
character(0)
> try(system("convert tmp/3j5f91291657082.ps tmp/3j5f91291657082.png",intern=TRUE))
character(0)
> try(system("convert tmp/4uwet1291657082.ps tmp/4uwet1291657082.png",intern=TRUE))
character(0)
> try(system("convert tmp/5uwet1291657082.ps tmp/5uwet1291657082.png",intern=TRUE))
character(0)
> try(system("convert tmp/6uwet1291657082.ps tmp/6uwet1291657082.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.669   1.232   1.870