R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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.

  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(6827,6178,7084,8162,8462,9644,10466,10748,9963,8194,6848,7027,7269,6775,7819,8371,9069,10248,11030,10882,10333,9109,7685,7602,8350,7829,8829,9948,10638,11253,11424,11391,10665,9396,7775,7933,8186,7444,8484,9948,10252,12282,11637,11577,12417,9637,8094,9280,8334,7899,9994,10078,10801,12950,12222,12246,13281,10366,8730,9614,8639,8772,10894,10455,11179,10588,10794,12770,13812,10857,9290,10925,9491,8919,11607,8852,12537,14759,13667,13731,15110,12185,10645,12161,10840,10436,13589,13402,13103,14933,14147,14057,16234,12389,11595,12772)
> 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] 96
> (np <- floor(n / par1))
[1] 8
> 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 8 8 8 8 8 8 8 8 8 8 8
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8] [,9]
 [1,]  6827  7269  8350  8186  8334  8639  9491 10840   NA
 [2,]  6178  6775  7829  7444  7899  8772  8919 10436   NA
 [3,]  7084  7819  8829  8484  9994 10894 11607 13589   NA
 [4,]  8162  8371  9948  9948 10078 10455  8852 13402   NA
 [5,]  8462  9069 10638 10252 10801 11179 12537 13103   NA
 [6,]  9644 10248 11253 12282 12950 10588 14759 14933   NA
 [7,] 10466 11030 11424 11637 12222 10794 13667 14147   NA
 [8,] 10748 10882 11391 11577 12246 12770 13731 14057   NA
 [9,]  9963 10333 10665 12417 13281 13812 15110 16234   NA
[10,]  8194  9109  9396  9637 10366 10857 12185 12389   NA
[11,]  6848  7685  7775  8094  8730  9290 10645 11595   NA
[12,]  7027  7602  7933  9280  9614 10925 12161 12772   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/1cx021304954020.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/2ulx31304954020.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/3nz441304954020.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/4o2w91304954020.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]
[1,]  6827.0  6178.0  7084.0  8162.0  8462.0  9644.0 10466.0 10748.0  9963
[2,]  7727.5  7109.5  8151.5  8611.5  9660.5 10418.0 10912.0 11136.5 10499
[3,]  8342.0  7864.0  9411.5  9948.0 10719.5 11767.5 11530.5 11911.5 12849
[4,]  9065.0  8845.5 11250.5 10266.5 11858.0 13854.5 12944.5 13250.5 14461
[5,] 10840.0 10436.0 13589.0 10455.0 13103.0 14933.0 14147.0 14057.0 16234
       [,10]   [,11]   [,12]
[1,]  8194.0  6848.0  7027.0
[2,]  9252.5  7730.0  7767.5
[3,] 10001.5  8412.0  9447.0
[4,] 11521.0  9967.5 11543.0
[5,] 12389.0 11595.0 12772.0

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

$conf
         [,1]     [,2]      [,3]      [,4]      [,5]      [,6]     [,7]
[1,] 7594.853 6894.245  7680.354  9023.493  9491.945  9847.822 10395.12
[2,] 9089.147 8833.755 11142.646 10872.507 11947.055 13687.178 12665.88
         [,8]     [,9]     [,10]  [,11]     [,12]
[1,] 10730.59 10635.77  8734.283 7162.1  7337.951
[2,] 13092.41 15062.23 11268.717 9661.9 11556.049

$out
[1] 13402

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(6827, 7727.5, 8342, 9065, 10840, 6178, 7109.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5u7oy1304954020.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]
[1,]  6178.0  6775.0  7775.0  7444.0  7899  8639.0  8852 10436.0   NA
[2,]  6937.5  7643.5  8141.5  8335.0  9172  9872.5 10068 11992.0   NA
[3,]  8178.0  8720.0  9672.0  9792.5 10222 10825.5 12173 13252.5   NA
[4,]  9803.5 10290.5 10959.0 11607.0 12234 11052.0 13699 14102.0   NA
[5,] 10748.0 11030.0 11424.0 12417.0 13281 12770.0 15110 16234.0   NA

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

$conf
         [,1]     [,2]     [,3]      [,4]      [,5]     [,6]     [,7]     [,8]
[1,] 6870.798 7512.686  8386.92  8300.119  8825.401 10287.52 10516.88 12290.11
[2,] 9485.202 9927.314 10957.08 11284.881 11618.599 11363.48 13829.12 14214.89
     [,9]
[1,]   NA
[2,]   NA

$out
[1] 13812

$group
[1] 6

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

Warning message:
In bxp(list(stats = c(6178, 6937.5, 8178, 9803.5, 10748, 6775, 7643.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/665ry1304954020.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,] 1250.110 3309  895.250
[2,] 1340.000 4104 1547.125
[3,] 1587.631 4694 1805.625
[4,] 2071.917 5517 2826.875
[5,] 2304.677 6505 3554.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1253.798 4049.521 1221.922
[2,] 1921.463 5338.479 2389.328

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1250.11039512517, 1339.99981044292, 1587.63051014139,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1cx021304954020.ps tmp/1cx021304954020.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ulx31304954020.ps tmp/2ulx31304954020.png",intern=TRUE))
character(0)
> try(system("convert tmp/3nz441304954020.ps tmp/3nz441304954020.png",intern=TRUE))
character(0)
> try(system("convert tmp/4o2w91304954020.ps tmp/4o2w91304954020.png",intern=TRUE))
character(0)
> try(system("convert tmp/5u7oy1304954020.ps tmp/5u7oy1304954020.png",intern=TRUE))
character(0)
> try(system("convert tmp/665ry1304954020.ps tmp/665ry1304954020.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.650   0.178   1.949