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(90,69.3,87.3,57.4,56.2,61.6,77.7,177.2,97.6,81.6,96.8,191.3,106,75.1,72,63.5,57.4,62.3,79.4,178.1,109.3,85.2,102.7,193.7,108.4,73.4,85.9,58.5,58.6,62.7,77.5,180.5,102.2,82.6,97.8,197.8,93.8,72.4,77.7,58.7,53.1,64.3,76.4,188.4,105.5,79.8,96.1,202.5,97.3,89.5,64.7,61.2,57.8,62,76.3,195,110.9,81.4,101.7,202.2,97.4,68.5,86.8,59.1,62.4,66.2,68,198.5,120.4,90.2,103.2,207.3,106.4,75.5,97.3,60,67.5,71.2,73.7,213.3,114.6,96.1,117,229.2,105.6,99.9,79.3,72.5,67.4,78.3,85.7,177.4,113.6,94.1,105.7,228.3,100.3,70.3,94.2,66.5,64.4,73.7,87.9,152.2,97.3,89.3,107.6,228.4)
> 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,]  90.0 106.0 108.4  93.8  97.3  97.4 106.4 105.6 100.3    NA
 [2,]  69.3  75.1  73.4  72.4  89.5  68.5  75.5  99.9  70.3    NA
 [3,]  87.3  72.0  85.9  77.7  64.7  86.8  97.3  79.3  94.2    NA
 [4,]  57.4  63.5  58.5  58.7  61.2  59.1  60.0  72.5  66.5    NA
 [5,]  56.2  57.4  58.6  53.1  57.8  62.4  67.5  67.4  64.4    NA
 [6,]  61.6  62.3  62.7  64.3  62.0  66.2  71.2  78.3  73.7    NA
 [7,]  77.7  79.4  77.5  76.4  76.3  68.0  73.7  85.7  87.9    NA
 [8,] 177.2 178.1 180.5 188.4 195.0 198.5 213.3 177.4 152.2    NA
 [9,]  97.6 109.3 102.2 105.5 110.9 120.4 114.6 113.6  97.3    NA
[10,]  81.6  85.2  82.6  79.8  81.4  90.2  96.1  94.1  89.3    NA
[11,]  96.8 102.7  97.8  96.1 101.7 103.2 117.0 105.7 107.6    NA
[12,] 191.3 193.7 197.8 202.5 202.2 207.3 229.2 228.3 228.4    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/1x9wm1291633017.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/271vq1291633017.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/371vq1291633017.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/4iava1291633017.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,]  90.0 68.5 64.7 57.4 53.1 61.6 73.7 152.2  97.3  79.8  96.1 191.3
[2,]  97.3 70.3 77.7 58.7 57.4 62.3 76.3 177.4 102.2  81.6  97.8 197.8
[3,] 100.3 73.4 85.9 60.0 58.6 64.3 77.5 180.5 109.3  85.2 102.7 202.5
[4,] 106.0 75.5 87.3 63.5 64.4 71.2 79.4 195.0 113.6  90.2 105.7 228.3
[5,] 108.4 75.5 97.3 66.5 67.5 78.3 79.4 213.3 120.4  96.1 117.0 229.2

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

$conf
        [,1]     [,2]   [,3]   [,4]     [,5]     [,6]     [,7]     [,8]    [,9]
[1,]  95.718 70.66133 80.844 57.472 54.91333 59.61267 75.86733 171.2307 103.296
[2,] 104.882 76.13867 90.956 62.528 62.28667 68.98733 79.13267 189.7693 115.304
        [,10]     [,11]    [,12]
[1,] 80.67067  98.53933 186.4367
[2,] 89.72933 106.86067 218.5633

$out
[1] 89.5 99.9 72.5 68.0 85.7 87.9

$group
[1] 2 2 4 7 7 7

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

Warning message:
In bxp(list(stats = c(90, 97.3, 100.3, 106, 108.4, 68.5, 70.3, 73.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5iava1291633017.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,] 56.20  57.40  58.50  53.10  57.80  59.1  60.00  67.40  64.40    NA
[2,] 65.45  67.75  68.05  68.35  63.35  67.1  72.45  78.80  72.00    NA
[3,] 84.45  82.30  84.25  78.75  85.45  88.5  96.70  97.00  91.75    NA
[4,] 97.20 107.65 105.30 100.80 106.30 111.8 115.80 109.65 103.95    NA
[5,] 97.60 109.30 108.40 105.50 110.90 120.4 117.00 113.60 107.60    NA

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

$conf
         [,1]      [,2]      [,3]     [,4]      [,5]      [,6]      [,7]
[1,] 69.96861  64.10134  67.26002 63.94934  65.86022  68.11203  76.92777
[2,] 98.93139 100.49866 101.23998 93.55066 105.03978 108.88797 116.47223
         [,8]      [,9] [,10]
[1,]  82.9291  77.17739    NA
[2,] 111.0709 106.32261    NA

$out
 [1] 177.2 191.3 178.1 193.7 180.5 197.8 188.4 202.5 195.0 202.2 198.5 207.3
[13] 213.3 229.2 177.4 228.3 152.2 228.4

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

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

Warning message:
In bxp(list(stats = c(56.2, 65.45, 84.45, 97.2, 97.6, 57.4, 67.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6iava1291633017.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.874167 14.4  3.10
[2,]  5.949789 16.5  6.10
[3,]  6.458951 20.4  8.65
[4,] 10.504302 32.0 10.50
[5,] 17.156809 37.9 11.40

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4.381608 13.33035  6.64313
[2,] 8.536295 27.46965 10.65687

$out
[1] 61.1 17.6 30.5

$group
[1] 2 3 3

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

Warning message:
In bxp(list(stats = c(4.87416659542942, 5.94978907767772, 6.45895148736613,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1x9wm1291633017.ps tmp/1x9wm1291633017.png",intern=TRUE))
character(0)
> try(system("convert tmp/271vq1291633017.ps tmp/271vq1291633017.png",intern=TRUE))
character(0)
> try(system("convert tmp/371vq1291633017.ps tmp/371vq1291633017.png",intern=TRUE))
character(0)
> try(system("convert tmp/4iava1291633017.ps tmp/4iava1291633017.png",intern=TRUE))
character(0)
> try(system("convert tmp/5iava1291633017.ps tmp/5iava1291633017.png",intern=TRUE))
character(0)
> try(system("convert tmp/6iava1291633017.ps tmp/6iava1291633017.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.440   0.900   2.328