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.

  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(10236,10893,10756,10940,10997,10827,10166,10186,10457,10368,10244,10511,10812,10738,10171,9721,9897,9828,9924,10371,10846,10413,10709,10662,10570,10297,10635,10872,10296,10383,10431,10574,10653,10805,10872,10625,10407,10463,10556,10646,10702,11353,11346,11451,11964,12574,13031,13812,14544,14931,14886,16005,17064,15168,16050,15839,15137,14954,15648,15305,15579,16348,15928,16171,15937,15713,15594,15683,16438,17032,17696,17745,19394)
> 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] 73
> (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 6 6 6 6 6 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 10236 10812 10570 10407 14544 15579 19394
 [2,] 10893 10738 10297 10463 14931 16348    NA
 [3,] 10756 10171 10635 10556 14886 15928    NA
 [4,] 10940  9721 10872 10646 16005 16171    NA
 [5,] 10997  9897 10296 10702 17064 15937    NA
 [6,] 10827  9828 10383 11353 15168 15713    NA
 [7,] 10166  9924 10431 11346 16050 15594    NA
 [8,] 10186 10371 10574 11451 15839 15683    NA
 [9,] 10457 10846 10653 11964 15137 16438    NA
[10,] 10368 10413 10805 12574 14954 17032    NA
[11,] 10244 10709 10872 13031 15648 17696    NA
[12,] 10511 10662 10625 13812 15305 17745    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/17jzm1211844486.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/2h60l1211844486.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/3zp291211844486.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/4rg6y1211844486.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]
[1,] 10236.0 10297.0 10171.0  9721  9897.0  9828  9924.0 10186.0 10457 10368.0
[2,] 10488.5 10463.0 10556.0 10646 10296.0 10383 10166.0 10371.0 10653 10413.0
[3,] 10812.0 10815.5 10695.5 10906 10849.5 11090 10888.5 11012.5 11405 11689.5
[4,] 15061.5 14931.0 14886.0 16005 15937.0 15168 15594.0 15683.0 15137 14954.0
[5,] 19394.0 16348.0 15928.0 16171 17064.0 15713 16050.0 15839.0 16438 17032.0
       [,11] [,12]
[1,] 10244.0 10511
[2,] 10709.0 10625
[3,] 11951.5 12237
[4,] 15648.0 15305
[5,] 17696.0 17745

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

$conf
          [,1]      [,2]     [,3]      [,4]      [,5]     [,6]      [,7]
[1,]  8081.078  7933.496  7902.51  7449.272  7210.873  8003.52  7387.265
[2,] 13542.922 13697.504 13488.49 14362.728 14488.127 14176.48 14389.735
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,]  7586.088  8512.675  8760.408  8765.685  9218.249
[2,] 14438.912 14297.325 14618.592 15137.315 15255.751

$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(10236, 10488.5, 10812, 15061.5, 19394, 10297,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5gs9j1211844486.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,] 10166  9721.0 10296.0 10407.0 14544.0 15579 19394
[2,] 10240  9910.5 10407.0 10601.0 14942.5 15698 19394
[3,] 10484 10392.0 10599.5 11349.5 15236.5 16054 19394
[4,] 10860 10723.5 10729.0 12269.0 15922.0 16735 19394
[5,] 10997 10846.0 10872.0 13812.0 17064.0 17745 19394

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]  [,7]
[1,] 10201.21 10021.19 10452.63 10588.71 14789.74 15581.02 19394
[2,] 10766.79 10762.81 10746.37 12110.29 15683.26 16526.98 19394

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(10166, 10240, 10484, 10860, 10997, 9721, 9910.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wtby1211844486.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,] 2547.726 5653.0 3277.750
[2,] 2610.074 5933.0 3681.375
[3,] 2785.987 6288.0 4119.750
[4,] 3035.566 7200.5 4298.625
[5,] 3533.932 7452.0 4573.000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2591.917 5709.885 3838.218
[2,] 2980.057 6866.115 4401.282

$out
[1] 9158

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(2547.72633276549, 2610.07370116771, 2785.98707788957,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/17jzm1211844486.ps tmp/17jzm1211844486.png")
> system("convert tmp/2h60l1211844486.ps tmp/2h60l1211844486.png")
> system("convert tmp/3zp291211844486.ps tmp/3zp291211844486.png")
> system("convert tmp/4rg6y1211844486.ps tmp/4rg6y1211844486.png")
> system("convert tmp/5gs9j1211844486.ps tmp/5gs9j1211844486.png")
> system("convert tmp/6wtby1211844486.ps tmp/6wtby1211844486.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.365   1.360   2.608