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(8357.00
+ ,7454.00
+ ,8076.00
+ ,7248.00
+ ,7339.00
+ ,7292.00
+ ,7359.00
+ ,7537.00
+ ,7441.00
+ ,8057.00
+ ,8037.00
+ ,8257.00
+ ,8692.00
+ ,8119.00
+ ,8236.00
+ ,7432.00
+ ,7669.00
+ ,7453.00
+ ,7566.00
+ ,7731.00
+ ,7657.00
+ ,8130.00
+ ,8401.00
+ ,8737.00
+ ,9009.00
+ ,7919.00
+ ,8228.00
+ ,7903.00
+ ,7912.00
+ ,7857.00
+ ,7965.00
+ ,8091.00
+ ,8024.00
+ ,8772.00
+ ,8656.00
+ ,8953.00
+ ,9014.00
+ ,8103.00
+ ,8876.00
+ ,8231.00
+ ,8173.00
+ ,8087.00
+ ,8296.00
+ ,8007.00
+ ,8382.00
+ ,9168.00
+ ,9137.00
+ ,9321.00
+ ,9234.00
+ ,8451.00
+ ,9101.00
+ ,8279.00
+ ,8284.00
+ ,8225.00
+ ,8597.00
+ ,8305.00
+ ,8620.00
+ ,9102.00
+ ,9258.00
+ ,9652.00
+ ,9522.00
+ ,8874.00
+ ,9415.00
+ ,8525.00
+ ,8862.00
+ ,8421.00
+ ,8626.00
+ ,8750.00
+ ,8852.00
+ ,9412.00
+ ,9570.00
+ ,9513.00
+ ,9986.00
+ ,8907.00
+ ,9663.00
+ ,8799.00
+ ,8931.00
+ ,8732.00
+ ,8936.00
+ ,9127.00
+ ,9070.00
+ ,9773.00
+ ,9670.00
+ ,9929.00
+ ,10095.00
+ ,9025.00
+ ,9659.00
+ ,8954.00
+ ,9022.00
+ ,8855.00
+ ,9034.00
+ ,9196.00
+ ,9038.00
+ ,9650.00
+ ,9715.00
+ ,10052.00
+ ,10436.00
+ ,9314.00
+ ,9717.00
+ ,8997.00
+ ,9062.00
+ ,8885.00
+ ,9058.00
+ ,9095.00
+ ,9149.00
+ ,9857.00
+ ,9848.00
+ ,10269.00
+ ,10341.00
+ ,9690.00
+ ,10125.00
+ ,9349.00
+ ,9224.00
+ ,9224.00
+ ,9454.00
+ ,9347.00
+ ,9430.00
+ ,9933.00
+ ,10148.00
+ ,10677.00
+ ,10735.00
+ ,9760.00
+ ,10567.00
+ ,9333.00
+ ,9409.00
+ ,9502.00
+ ,9348.00
+ ,9319.00
+ ,9594.00
+ ,10160.00
+ ,10182.00
+ ,10810.00
+ ,11105.00
+ ,9874.00
+ ,10958.00
+ ,9311.00
+ ,9610.00
+ ,9398.00
+ ,9784.00
+ ,9425.00
+ ,9557.00
+ ,10166.00
+ ,10337.00
+ ,10770.00
+ ,11265.00
+ ,10183.00
+ ,10941.00
+ ,9628.00
+ ,9709.00
+ ,9637.00
+ ,9579.00
+ ,9741.00
+ ,9754.00
+ ,10508.00
+ ,10749.00
+ ,11079.00
+ ,11608.00
+ ,10668.00
+ ,10933.00
+ ,9703.00
+ ,9799.00
+ ,9656.00
+ ,9648.00
+ ,9712.00
+ ,9766.00
+ ,10540.00
+ ,10564.00
+ ,10911.00
+ ,11218.00
+ ,10230.00
+ ,10410.00
+ ,9227.00
+ ,9378.00
+ ,9105.00
+ ,9128.00)
> 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] 175
> (np <- floor(n / par1))
[1] 14
> 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] 15 15 15 15 15 15 15 14 14 14 14 14
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]  [,8]  [,9] [,10] [,11] [,12] [,13]
 [1,] 8357 8692 9009 9014 9234 9522 9986 10095 10436 10341 10735 11105 11265
 [2,] 7454 8119 7919 8103 8451 8874 8907  9025  9314  9690  9760  9874 10183
 [3,] 8076 8236 8228 8876 9101 9415 9663  9659  9717 10125 10567 10958 10941
 [4,] 7248 7432 7903 8231 8279 8525 8799  8954  8997  9349  9333  9311  9628
 [5,] 7339 7669 7912 8173 8284 8862 8931  9022  9062  9224  9409  9610  9709
 [6,] 7292 7453 7857 8087 8225 8421 8732  8855  8885  9224  9502  9398  9637
 [7,] 7359 7566 7965 8296 8597 8626 8936  9034  9058  9454  9348  9784  9579
 [8,] 7537 7731 8091 8007 8305 8750 9127  9196  9095  9347  9319  9425  9741
 [9,] 7441 7657 8024 8382 8620 8852 9070  9038  9149  9430  9594  9557  9754
[10,] 8057 8130 8772 9168 9102 9412 9773  9650  9857  9933 10160 10166 10508
[11,] 8037 8401 8656 9137 9258 9570 9670  9715  9848 10148 10182 10337 10749
[12,] 8257 8737 8953 9321 9652 9513 9929 10052 10269 10677 10810 10770 11079
      [,14] [,15]
 [1,] 11608 11218
 [2,] 10668 10230
 [3,] 10933 10410
 [4,]  9703  9227
 [5,]  9799  9378
 [6,]  9656  9105
 [7,]  9648  9128
 [8,]  9712    NA
 [9,]  9766    NA
[10,] 10540    NA
[11,] 10564    NA
[12,] 10911    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/1hx9i1262554486.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/239471262554486.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/3dyk41262554486.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/4cltd1262554486.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]
[1,]  8357  7454  8076.0 7248 7339.0 7292 7359.0 7537 7441  8057.0  8037.0
[2,]  9124  8285  8988.5 8255 8228.5 8156 8446.5 8091 8382  9102.0  9137.0
[3,] 10095  9025  9663.0 8954 9022.0 8855 9034.0 9111 9054  9711.5  9692.5
[4,] 10920  9817 10488.5 9322 9393.5 9311 9401.0 9347 9557 10160.0 10182.0
[5,] 11608 10668 10958.0 9703 9799.0 9656 9784.0 9741 9766 10540.0 10749.0
       [,12]
[1,]  8257.0
[2,]  9321.0
[3,]  9990.5
[4,] 10770.0
[5,] 11079.0

$n
 [1] 15 15 15 15 15 15 15 14 14 14 14 14

$conf
          [,1]     [,2]      [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,]  9362.314 8400.014  9051.069 8518.713 8546.733 8383.813 8644.608 8580.625
[2,] 10827.686 9649.986 10274.931 9389.287 9497.267 9326.187 9423.392 9641.375
        [,9]     [,10]     [,11]     [,12]
[1,] 8557.83  9264.735  9251.225  9378.627
[2,] 9550.17 10158.265 10133.775 10602.373

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(8357, 9124, 10095, 10920, 11608, 7454, 8285,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5axun1262554486.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]    [,9]   [,10]
[1,] 7248.0 7432.0 7857.0 8007.0 8225.0 8421 8732.0  8855.0  8885.0  9224.0
[2,] 7349.0 7611.5 7915.5 8138.0 8294.5 8688 8919.0  9023.5  9060.0  9348.0
[3,] 7495.5 7925.0 8057.5 8339.0 8608.5 8868 9098.5  9117.0  9231.5  9572.0
[4,] 8066.5 8318.5 8714.0 9075.5 9168.0 9464 9721.5  9687.0  9852.5 10136.5
[5,] 8357.0 8737.0 9009.0 9321.0 9652.0 9570 9986.0 10095.0 10436.0 10677.0
       [,11]   [,12]   [,13]   [,14]   [,15]
[1,]  9319.0  9311.0  9579.0  9648.0  9105.0
[2,]  9378.5  9491.0  9673.0  9707.5  9177.5
[3,]  9677.0  9829.0  9968.5 10169.5  9378.0
[4,] 10374.5 10553.5 10845.0 10789.5 10320.0
[5,] 10810.0 11105.0 11265.0 11608.0 11218.0

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

$conf
         [,1]     [,2]     [,3]   [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 7168.243 7602.533 7693.299 7911.4 8210.09 8514.061 8732.474 8814.373
[2,] 7822.757 8247.467 8421.701 8766.6 9006.91 9221.939 9464.526 9419.627
         [,9]   [,10]     [,11]     [,12]     [,13]     [,14]     [,15]
[1,] 8870.035 9212.36  9222.718  9344.387  9433.943  9675.993  8695.717
[2,] 9592.965 9931.64 10131.282 10313.613 10503.057 10663.007 10060.283

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(7248, 7349, 7495.5, 8066.5, 8357, 7432, 7611.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6iomk1262554486.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,]  741.1936 2204.0  954.500
[2,]  765.7996 2394.5 1036.625
[3,]  782.1865 2471.5 1160.000
[4,]  924.1287 2852.0 1438.875
[5,] 1031.6669 3251.0 1796.000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 709.9715 2262.831  976.531
[2,] 854.4014 2680.169 1343.469

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(741.193576322959, 765.799643669446, 782.186484806354,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1hx9i1262554486.ps tmp/1hx9i1262554486.png",intern=TRUE))
character(0)
> try(system("convert tmp/239471262554486.ps tmp/239471262554486.png",intern=TRUE))
character(0)
> try(system("convert tmp/3dyk41262554486.ps tmp/3dyk41262554486.png",intern=TRUE))
character(0)
> try(system("convert tmp/4cltd1262554486.ps tmp/4cltd1262554486.png",intern=TRUE))
character(0)
> try(system("convert tmp/5axun1262554486.ps tmp/5axun1262554486.png",intern=TRUE))
character(0)
> try(system("convert tmp/6iomk1262554486.ps tmp/6iomk1262554486.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.083   0.801   1.849