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(18288.3,16049,16764.5,17880.2,16555.9,16087.1,16373.5,17842.2,22321.5,22786.7,18274.1,22392.9,23899.3,21343.5,22952.3,21374.4,21164.1,20906.5,17877.4,20664.3,22160,19813.6,17735.4,19640.2,20844.4,19823.1,18594.6,21350.6,18574.1,18924.2,17343.4,19961.2,19932.1,19464.6,16165.4,17574.9,19795.4,19439.5,17170,21072.4,17751.8,17515.5,18040.3,19090.1,17746.5,19202.1,15141.6,16258.1,18586.5,17209.4,17838.7,19123.5,16583.6,15991.2,16704.5,17422,17872,17823.2,13866.5,15912.8,17870.5,15420.3,16379.4,17903.9,15305.8,14583.3,14861,14968.9,16726.5,16283.6,11703.7,15101.8,15469.7,14956.9,15370.6,15998.1,14725.1,14768.9,13659.6,15070.3,16942.6,15761.3,12083,15023.6,15106.5)
> 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] 85
> (np <- floor(n / par1))
[1] 7
> 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] 8 7 7 7 7 7 7 7 7 7 7 7
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
 [1,] 18288.3 23899.3 20844.4 19795.4 18586.5 17870.5 15469.7 15106.5
 [2,] 16049.0 21343.5 19823.1 19439.5 17209.4 15420.3 14956.9      NA
 [3,] 16764.5 22952.3 18594.6 17170.0 17838.7 16379.4 15370.6      NA
 [4,] 17880.2 21374.4 21350.6 21072.4 19123.5 17903.9 15998.1      NA
 [5,] 16555.9 21164.1 18574.1 17751.8 16583.6 15305.8 14725.1      NA
 [6,] 16087.1 20906.5 18924.2 17515.5 15991.2 14583.3 14768.9      NA
 [7,] 16373.5 17877.4 17343.4 18040.3 16704.5 14861.0 13659.6      NA
 [8,] 17842.2 20664.3 19961.2 19090.1 17422.0 14968.9 15070.3      NA
 [9,] 22321.5 22160.0 19932.1 17746.5 17872.0 16726.5 16942.6      NA
[10,] 22786.7 19813.6 19464.6 19202.1 17823.2 16283.6 15761.3      NA
[11,] 18274.1 17735.4 16165.4 15141.6 13866.5 11703.7 12083.0      NA
[12,] 22392.9 19640.2 17574.9 16258.1 15912.8 15101.8 15023.6      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/1765y1274033038.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/2zx401274033038.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/3aomm1274033038.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/4aomm1274033038.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]
[1,] 15106.5 14956.90 15370.60 15998.10 14725.10 14583.30 13659.60 14968.90
[2,] 16670.1 15734.65 16571.95 17892.05 15930.85 15380.05 15617.25 16246.15
[3,] 18437.4 17209.40 17170.00 19123.50 16583.60 16087.10 16704.50 17842.20
[4,] 20319.9 19631.30 18216.65 21211.50 18162.95 18219.85 17610.40 19525.65
[5,] 23899.3 21343.50 18594.60 21374.40 21164.10 20906.50 18040.30 20664.30
         [,9]   [,10]    [,11]    [,12]
[1,] 16726.50 15761.3 11703.70 15023.60
[2,] 17344.55 17053.4 12974.75 15507.30
[3,] 17872.00 19202.1 15141.60 16258.10
[4,] 21046.05 19639.1 16950.40 18607.55
[5,] 22321.50 22786.7 18274.10 22392.90

$n
 [1] 8 7 7 7 7 7 7 7 7 7 7 7

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 16398.57 14882.38 16187.81 17141.18 15250.63 14391.22 15514.22 15883.74
[2,] 20476.23 19536.42 18152.19 21105.82 17916.57 17782.98 17894.78 19800.66
         [,9]    [,10]    [,11]    [,12]
[1,] 15661.52 17657.96 12767.41 14406.68
[2,] 20082.48 20746.24 17515.79 18109.52

$out
[1] 22952.3

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(15106.5, 16670.1, 18437.4, 20319.9, 23899.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5aomm1274033038.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]     [,7]    [,8]
[1,] 16049.0 17735.4 16165.40 15141.60 15912.80 14583.30 13659.60 15106.5
[2,] 16464.7 19726.9 18074.50 17342.75 16287.40 14914.95 14747.00 15106.5
[3,] 17861.2 21035.3 19194.40 17896.05 17315.70 15363.05 15046.95 15106.5
[4,] 20304.9 21767.2 19946.65 19320.80 17855.35 16552.95 15615.50 15106.5
[5,] 22786.7 23899.3 21350.60 21072.40 19123.50 17903.90 15998.10 15106.5

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

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 16109.66 20104.71 18340.5 16993.85 16600.55 14615.95 14650.82 15106.5
[2,] 19612.74 21965.89 20048.3 18798.25 18030.85 16110.15 15443.08 15106.5

$out
[1] 13866.5 11703.7 16942.6 12083.0

$group
[1] 5 6 7 7

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

Warning message:
In bxp(list(stats = c(16049, 16464.7, 17861.2, 20304.9, 22786.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6aomm1274033038.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,] 2102.602 4380.70 1644.700
[2,] 2208.273 5645.20 2408.900
[3,] 2381.800 6412.80 2970.025
[4,] 2528.696 7197.35 3510.475
[5,] 2854.598 8792.80 3975.650

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2235.653 5704.854 2467.589
[2,] 2527.946 7120.746 3472.461

$out
[1] 1619.389

$group
[1] 1

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

> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1765y1274033038.ps tmp/1765y1274033038.png",intern=TRUE))
character(0)
> try(system("convert tmp/2zx401274033038.ps tmp/2zx401274033038.png",intern=TRUE))
character(0)
> try(system("convert tmp/3aomm1274033038.ps tmp/3aomm1274033038.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aomm1274033038.ps tmp/4aomm1274033038.png",intern=TRUE))
character(0)
> try(system("convert tmp/5aomm1274033038.ps tmp/5aomm1274033038.png",intern=TRUE))
character(0)
> try(system("convert tmp/6aomm1274033038.ps tmp/6aomm1274033038.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.068   0.824   1.582