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(36845,35338,35022,34777,26887,23970,22780,17351,21382,24561,17409,11514,31514,27071,29462,26105,22397,23843,21705,18089,20764,25316,17704,15548,28029,29383,36438,32034,22679,24319,18004,17537,20366,22782,19169,13807,29743,25591,29096,26482,22405,27044,17970,18730,19684,19785,18479,10698,31956,29506,34506,27165,26736,23691,18157,17328,18205,20995,17382,9367,31124,26551,30651,25859,25100,25778,20418,18688,20424,24776,19814,12738,42553)
> 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,] 36845 31514 28029 29743 31956 31124 42553
 [2,] 35338 27071 29383 25591 29506 26551    NA
 [3,] 35022 29462 36438 29096 34506 30651    NA
 [4,] 34777 26105 32034 26482 27165 25859    NA
 [5,] 26887 22397 22679 22405 26736 25100    NA
 [6,] 23970 23843 24319 27044 23691 25778    NA
 [7,] 22780 21705 18004 17970 18157 20418    NA
 [8,] 17351 18089 17537 18730 17328 18688    NA
 [9,] 21382 20764 20366 19684 18205 20424    NA
[10,] 24561 25316 22782 19785 20995 24776    NA
[11,] 17409 17704 19169 18479 17382 19814    NA
[12,] 11514 15548 13807 10698  9367 12738    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/1nl7v1210605423.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/2lynq1210605423.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/3zuzk1210605423.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/4qqgg1210605423.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,] 28029.0 25591 29096.0 25859.0 22397.0 23691.0 17970.0 17328 18205 19785.0
[2,] 30433.5 26551 29462.0 26105.0 22405.0 23843.0 18004.0 17351 19684 20995.0
[3,] 31514.0 28227 32578.5 26823.5 23889.5 24144.5 19287.5 17813 20395 23671.5
[4,] 34400.5 29506 35022.0 32034.0 26736.0 25778.0 21705.0 18688 20764 24776.0
[5,] 36845.0 29506 36438.0 34777.0 26887.0 27044.0 22780.0 18730 21382 25316.0
       [,11] [,12]
[1,] 17382.0  9367
[2,] 17409.0 10698
[3,] 18091.5 12126
[4,] 19169.0 13807
[5,] 19814.0 15548

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 29144.97 26320.93 28992.12 22999.10 21095.87 22896.36 16900.24 16950.59
[2,] 33883.03 30133.07 36164.88 30647.90 26683.13 25392.64 21674.76 18675.41
         [,9]    [,10]    [,11]    [,12]
[1,] 19698.37 21232.63 16956.24 10120.59
[2,] 21091.63 26110.37 19226.76 14131.41

$out
[1] 42553 35338

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(28029, 30433.5, 31514, 34400.5, 36845, 25591,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5qx1r1210605423.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,] 11514.0 15548.0 13807.0 10698.0  9367.0 12738 42553
[2,] 19395.5 19426.5 18586.5 18604.5 17769.5 20116 42553
[3,] 24265.5 23120.0 22730.5 21095.0 22343.0 24938 42553
[4,] 34899.5 26588.0 28706.0 26763.0 28335.5 26205 42553
[5,] 36845.0 31514.0 36438.0 29743.0 34506.0 31124 42553

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]  [,7]
[1,] 17194.02 19853.59 18114.93 17373.85 17523.78 22160.77 42553
[2,] 31336.98 26386.41 27346.07 24816.15 27162.22 27715.23 42553

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(11514, 19395.5, 24265.5, 34899.5, 36845, 15548,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6lbkl1210605423.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,]  646.3019  1402.0  824.500
[2,] 1221.2731  3265.0 1526.125
[3,] 2187.7042  5170.5 3037.375
[4,] 3341.7796  8130.0 3910.250
[5,] 4968.1795 14524.0 5133.750

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1220.527 2951.541 1949.960
[2,] 3154.881 7389.459 4124.790

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(646.301916032025, 1221.27314493111, 2187.70419722014,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1nl7v1210605423.ps tmp/1nl7v1210605423.png")
> system("convert tmp/2lynq1210605423.ps tmp/2lynq1210605423.png")
> system("convert tmp/3zuzk1210605423.ps tmp/3zuzk1210605423.png")
> system("convert tmp/4qqgg1210605423.ps tmp/4qqgg1210605423.png")
> system("convert tmp/5qx1r1210605423.ps tmp/5qx1r1210605423.png")
> system("convert tmp/6lbkl1210605423.ps tmp/6lbkl1210605423.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.392   1.359   2.588