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(101.3,106.3,94,102.8,102,105.1,92.4,81.4,105.8,120.3,100.7,88.8,94.3,99.9,103.4,103.3,98.8,104.2,91.2,74.7,108.5,114.5,96.9,89.6,97.1,100.3,122.6,115.4,109,129.1,102.8,96.2,127.7,128.9,126.5,119.8,113.2,114.1,134.1,130,121.8,132.1,105.3,103,117.1,126.3,138.1,119.5,138,135.5,178.6,162.2,176.9,204.9,132.2,142.5,164.3,174.9,175.4,143)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> ylab = ''
> xlab = ''
> main = ''
> #'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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,] 101.3  94.3  97.1 113.2 138.0   NA
 [2,] 106.3  99.9 100.3 114.1 135.5   NA
 [3,]  94.0 103.4 122.6 134.1 178.6   NA
 [4,] 102.8 103.3 115.4 130.0 162.2   NA
 [5,] 102.0  98.8 109.0 121.8 176.9   NA
 [6,] 105.1 104.2 129.1 132.1 204.9   NA
 [7,]  92.4  91.2 102.8 105.3 132.2   NA
 [8,]  81.4  74.7  96.2 103.0 142.5   NA
 [9,] 105.8 108.5 127.7 117.1 164.3   NA
[10,] 120.3 114.5 128.9 126.3 174.9   NA
[11,] 100.7  96.9 126.5 138.1 175.4   NA
[12,]  88.8  89.6 119.8 119.5 143.0   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/15b831257864162.ps",horizontal=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/2a2jv1257864162.ps",horizontal=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/36fut1257864162.ps",horizontal=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/4mido1257864162.ps",horizontal=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,]  94.3  99.9  94.0 102.8  98.8 104.2  91.2  74.7 105.8 114.5  96.9  88.8
[2,]  97.1 100.3 103.4 103.3 102.0 105.1  92.4  81.4 108.5 120.3 100.7  89.6
[3,] 101.3 106.3 122.6 115.4 109.0 129.1 102.8  96.2 117.1 126.3 126.5 119.5
[4,] 113.2 114.1 134.1 130.0 121.8 132.1 105.3 103.0 127.7 128.9 138.1 119.8
[5,] 113.2 114.1 178.6 162.2 121.8 132.1 105.3 103.0 127.7 128.9 175.4 143.0

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
          [,1]      [,2]     [,3]      [,4]      [,5]     [,6]     [,7]
[1,]  89.92378  96.54895 100.9075  96.53385  95.00937 110.0219  93.6849
[2,] 112.67622 116.05105 144.2925 134.26615 122.99063 148.1781 111.9151
         [,8]     [,9]    [,10]    [,11]     [,12]
[1,]  80.9375 103.5333 120.2233 100.0733  98.16076
[2,] 111.4625 130.6667 132.3767 152.9267 140.83924

$out
[1] 138.0 135.5 176.9 204.9 132.2 142.5 164.3 174.9

$group
[1]  1  2  5  6  7  8  9 10

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

Warning message:
In bxp(list(stats = c(94.3, 97.1, 101.3, 113.2, 113.2, 99.9, 100.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5r96y1257864162.ps",horizontal=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]
[1,]  81.40  89.60  96.20 103.00 132.20   NA
[2,]  93.20  92.75 101.55 113.65 140.25   NA
[3,] 101.65  99.35 117.60 120.65 163.25   NA
[4,] 105.45 103.80 127.10 131.05 176.15   NA
[5,] 120.30 114.50 129.10 138.10 204.90   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]      [,2]     [,3]     [,4]     [,5] [,6]
[1,]  96.0627  94.31002 105.9465 112.7137 146.8758   NA
[2,] 107.2373 104.38998 129.2535 128.5863 179.6242   NA

$out
[1] 74.7

$group
[1] 2

$names
[1] "1" "2" "3" "4" "5" NA 

Warning message:
In bxp(list(stats = c(81.4, 93.2, 101.65, 105.45, 120.3, 89.6, 92.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6felg1257864162.ps",horizontal=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.74286  35.60  8.60
[2,] 20.43356  48.95 14.95
[3,] 24.39084  59.90 20.70
[4,] 31.98262  78.30 28.60
[5,] 41.14562 100.70 37.40

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 19.12323 46.51327 14.47414
[2,] 29.65844 73.28673 26.92586

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(14.7428626799547, 20.4335554741599, 24.3908392636398,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/15b831257864162.ps tmp/15b831257864162.png")
> system("convert tmp/2a2jv1257864162.ps tmp/2a2jv1257864162.png")
> system("convert tmp/36fut1257864162.ps tmp/36fut1257864162.png")
> system("convert tmp/4mido1257864162.ps tmp/4mido1257864162.png")
> system("convert tmp/5r96y1257864162.ps tmp/5r96y1257864162.png")
> system("convert tmp/6felg1257864162.ps tmp/6felg1257864162.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.087   0.834   1.536