R version 2.6.2 (2008-02-08)
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.

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(11835.70,11542.20,13093.70,11180.20,12035.70,12112.00,10875.20,9897.30,11672.10,12385.70,11405.60,9830.90,11025.10,10853.80,12252.60,11839.40,11669.10,11601.40,11178.40,9516.40,12102.80,12989.00,11610.20,10205.50,11356.20,11307.10,12648.60,11947.20,11714.10,12192.50,11268.80,9097.40,12639.80,13040.10,11687.30,11191.70,11391.90,11793.10,13933.20,12778.10,11810.30,13698.40,11956.60,10723.80,13938.90,13979.80,13807.40,12973.90,12509.80,12934.10,14908.30,13772.10,13012.60,14049.90,11816.50,11593.20,14466.20,13615.90,14733.90,13880.70,13527.50,13584.00,16170.20,13260.60,14741.90,15486.50,13154.50,12621.20,15031.60,15452.40,15428.00,13105.90,14716.80,14180.00,16202.20,14392.40,15140.60,15960.10,14729.90,13705.20,15728.50,17315.60,16152.80)
> 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] 83
> (np <- floor(n / par1))
[1] 6
> 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] 7 7 7 7 7 7 7 7 7 7 7 6
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]
 [1,] 11835.7 11025.1 11356.2 11391.9 12509.8 13527.5 14716.8
 [2,] 11542.2 10853.8 11307.1 11793.1 12934.1 13584.0 14180.0
 [3,] 13093.7 12252.6 12648.6 13933.2 14908.3 16170.2 16202.2
 [4,] 11180.2 11839.4 11947.2 12778.1 13772.1 13260.6 14392.4
 [5,] 12035.7 11669.1 11714.1 11810.3 13012.6 14741.9 15140.6
 [6,] 12112.0 11601.4 12192.5 13698.4 14049.9 15486.5 15960.1
 [7,] 10875.2 11178.4 11268.8 11956.6 11816.5 13154.5 14729.9
 [8,]  9897.3  9516.4  9097.4 10723.8 11593.2 12621.2 13705.2
 [9,] 11672.1 12102.8 12639.8 13938.9 14466.2 15031.6 15728.5
[10,] 12385.7 12989.0 13040.1 13979.8 13615.9 15452.4 17315.6
[11,] 11405.6 11610.2 11687.3 13807.4 14733.9 15428.0 16152.8
[12,]  9830.9 10205.5 11191.7 12973.9 13880.7 13105.9      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/1v76k1210621904.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/2b0mk1210621904.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/3rjbu1210621904.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/4v0xw1210621904.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,] 11025.10 10853.80 12252.60 11180.20 11669.10 11601.40 10875.20  9097.40
[2,] 11374.05 11424.65 12871.15 11893.30 11762.20 12152.25 11223.60  9706.85
[3,] 11835.70 11793.10 13933.20 12778.10 12035.70 13698.40 11816.50 10723.80
[4,] 13018.65 13259.05 15539.25 13516.35 13877.25 14768.20 12555.55 12107.20
[5,] 14716.80 14180.00 16202.20 14392.40 15140.60 15960.10 13154.50 13705.20
        [,9]    [,10]    [,11]   [,12]
[1,] 11672.1 12385.70 11405.60  9830.9
[2,] 12371.3 13014.55 11648.75 10205.5
[3,] 13938.9 13615.90 13807.40 12082.8
[4,] 14748.9 14716.10 15080.95 13105.9
[5,] 15728.5 15452.40 16152.80 13880.7

$n
 [1] 7 7 7 7 7 7 7 7 7 7 7 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 10853.57 10697.63 12339.85 11808.84 10772.63 12136.20 11021.08  9290.35
[2,] 12817.83 12888.57 15526.55 13747.36 13298.77 15260.60 12611.92 12157.25
         [,9]    [,10]    [,11]    [,12]
[1,] 12519.04 12599.76 11757.75 10211.95
[2,] 15358.76 14632.04 15857.05 13953.65

$out
[1] 14729.9 17315.6

$group
[1]  7 10

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

Warning message:
In bxp(list(stats = c(11025.1, 11374.05, 11835.7, 13018.65, 14716.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5l87a1210621905.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]
[1,]  9830.90  9516.40 11191.70 10723.8 11593.20 12621.20 13705.20
[2,] 11027.70 10939.45 11287.95 11801.7 12721.95 13207.55 14554.60
[3,] 11607.15 11605.80 11700.70 12876.0 13694.00 14162.95 15140.60
[4,] 12073.85 11971.10 12416.15 13870.3 14258.05 15440.20 16056.45
[5,] 13093.70 12989.00 13040.10 13979.8 14908.30 16170.20 17315.60

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 11129.99 11135.26 11186.12 11932.50 12993.37 13144.62 14425.14
[2,] 12084.31 12076.34 12215.28 13819.50 14394.63 15181.28 15856.06

$out
[1] 9097.4

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(9830.9, 11027.7, 11607.15, 12073.85, 13093.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6frqg1210621905.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,] 1150.413 3212.20 1331.950
[2,] 1356.960 3581.60 1673.075
[3,] 1585.496 3999.70 2246.325
[4,] 1706.292 4483.25 2618.400
[5,] 1983.270 4929.90 3432.200

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1426.163 3588.451 1815.156
[2,] 1744.829 4410.949 2677.494

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1150.41341078513, 1356.95970811162, 1585.49613146784,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1v76k1210621904.ps tmp/1v76k1210621904.png")
> system("convert tmp/2b0mk1210621904.ps tmp/2b0mk1210621904.png")
> system("convert tmp/3rjbu1210621904.ps tmp/3rjbu1210621904.png")
> system("convert tmp/4v0xw1210621904.ps tmp/4v0xw1210621904.png")
> system("convert tmp/5l87a1210621905.ps tmp/5l87a1210621905.png")
> system("convert tmp/6frqg1210621905.ps tmp/6frqg1210621905.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.269   0.804   1.532