R version 2.8.0 (2008-10-20)
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.

  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(68897,38683,44720,39525,45315,50380,40600,36279,42438,38064,31879,11379,70249,39253,47060,41697,38708,49267,39018,32228,40870,39383,34571,12066,70938,34077,45409,40809,37013,44953,37848,32745,43412,34931,33008,8620,68906,39556,50669,36432,40891,48428,36222,33425,39401,37967,34801,12657,69116,41519,51321,38529,41547,52073,38401,40898,40439,41888,37898,8771,68184,50530,47221,41756,45633,48138,39486,39341,41117,41629,29722,7054,56676,34870,35117,30169,30936,35699,33228,27733,33666,35429,27438,8170)
> 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] 84
> (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] 7 7 7 7 7 7 7 7 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7] [,8]
 [1,] 68897 70249 70938 68906 69116 68184 56676   NA
 [2,] 38683 39253 34077 39556 41519 50530 34870   NA
 [3,] 44720 47060 45409 50669 51321 47221 35117   NA
 [4,] 39525 41697 40809 36432 38529 41756 30169   NA
 [5,] 45315 38708 37013 40891 41547 45633 30936   NA
 [6,] 50380 49267 44953 48428 52073 48138 35699   NA
 [7,] 40600 39018 37848 36222 38401 39486 33228   NA
 [8,] 36279 32228 32745 33425 40898 39341 27733   NA
 [9,] 42438 40870 43412 39401 40439 41117 33666   NA
[10,] 38064 39383 34931 37967 41888 41629 35429   NA
[11,] 31879 34571 33008 34801 37898 29722 27438   NA
[12,] 11379 12066  8620 12657  8771  7054  8170   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/1fq0w1291551338.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/27zhh1291551338.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/37zhh1291551338.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/47zhh1291551338.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,] 68184.0 34077.0 44720.0 36432.0 30936.0 44953.0 36222 27733.0 39401.0
[2,] 68540.5 36776.5 45064.5 37480.5 37860.5 46545.5 37035 32486.5 39920.0
[3,] 68906.0 39253.0 47060.0 39525.0 40891.0 48428.0 38401 33425.0 40870.0
[4,] 69682.5 40537.5 48945.0 41253.0 43431.0 49823.5 39252 37810.0 41777.5
[5,] 70938.0 41519.0 51321.0 41756.0 45633.0 52073.0 40600 40898.0 43412.0
     [,10]   [,11]   [,12]
[1,] 34931 27438.0  7054.0
[2,] 36698 30800.5  8395.0
[3,] 38064 33008.0  8771.0
[4,] 40506 34686.0 11722.5
[5,] 41888 37898.0 12657.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 68224.02 37006.99 44742.63 37272.12 37564.39 46470.43 37077.04 30245.89
[2,] 69587.98 41499.01 49377.37 41777.88 44217.61 50385.57 39724.96 36604.11
         [,9]    [,10]    [,11]     [,12]
[1,] 39760.73 35789.92 30687.64  6783.871
[2,] 41979.27 40338.08 35328.36 10758.129

$out
[1] 56676 50530 35117 30169 35699 33228 33666

$group
[1] 1 2 3 4 6 7 9

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

Warning message:
In bxp(list(stats = c(68184, 68540.5, 68906, 69682.5, 70938, 34077,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/50rh21291551338.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]
[1,] 31879.0 32228.0 32745.0 33425.0 37898.0 29722.0 27438   NA
[2,] 37171.5 36639.5 33542.5 35511.5 38465.0 39413.5 28951   NA
[3,] 40062.5 39318.0 37430.5 38684.0 41208.5 41692.5 33447   NA
[4,] 45017.5 44378.5 44182.5 44659.5 46604.5 47679.5 35273   NA
[5,] 50380.0 49267.0 45409.0 50669.0 52073.0 50530.0 35699   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7] [,8]
[1,] 36483.89 35788.19 32577.52 34511.54 37496.02 37922.32 30563.49   NA
[2,] 43641.11 42847.81 42283.48 42856.46 44920.98 45462.68 36330.51   NA

$out
 [1] 68897 11379 70249 12066 70938  8620 68906 12657 69116  8771 68184  7054
[13] 56676  8170

$group
 [1] 1 1 2 2 3 3 4 4 5 5 6 6 7 7

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

Warning message:
In bxp(list(stats = c(31879, 37171.5, 40062.5, 45017.5, 50380, 32228,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/60rh21291551338.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,] 2177.045  5603.0 1142.00
[2,] 2945.701  8559.0 2747.50
[3,] 4309.357 12376.0 3766.75
[4,] 5234.582 15450.5 3883.00
[5,] 5439.600 16453.0 5570.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 3265.382  9232.74 3248.841
[2,] 5353.331 15519.26 4284.659

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2177.04491718214, 2945.70102721643, 4309.35661304251,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1fq0w1291551338.ps tmp/1fq0w1291551338.png",intern=TRUE))
character(0)
> try(system("convert tmp/27zhh1291551338.ps tmp/27zhh1291551338.png",intern=TRUE))
character(0)
> try(system("convert tmp/37zhh1291551338.ps tmp/37zhh1291551338.png",intern=TRUE))
character(0)
> try(system("convert tmp/47zhh1291551338.ps tmp/47zhh1291551338.png",intern=TRUE))
character(0)
> try(system("convert tmp/50rh21291551338.ps tmp/50rh21291551338.png",intern=TRUE))
character(0)
> try(system("convert tmp/60rh21291551338.ps tmp/60rh21291551338.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.668   1.210   1.869