R version 2.10.1 (2009-12-14)
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(89507,87562,85209,82360,79054,79069,107551,115759,115585,110260,103444,102303,101397,97994,94044,91159,87239,89235,118647,125620,125154,117529,109459,108483,107137,104699,100804,96066,91971,93228,120144,127233,127166,118194,109940,106683,102834,99882,96666,92540,88744,89321,115870,122401,122030,113802,105791,103076,98658,96945,92497,90687,88796,90015,113228,118711,117460,106556,97347,92657,93118,89037,83570,81693,75956,73993,97088,102394,96549,89727,82336,82653,82303,79596,74472,73562,66618,69029,89899,93774,90305,83799,80320,82497)
> 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,]  89507 101397 107137 102834  98658  93118 82303   NA
 [2,]  87562  97994 104699  99882  96945  89037 79596   NA
 [3,]  85209  94044 100804  96666  92497  83570 74472   NA
 [4,]  82360  91159  96066  92540  90687  81693 73562   NA
 [5,]  79054  87239  91971  88744  88796  75956 66618   NA
 [6,]  79069  89235  93228  89321  90015  73993 69029   NA
 [7,] 107551 118647 120144 115870 113228  97088 89899   NA
 [8,] 115759 125620 127233 122401 118711 102394 93774   NA
 [9,] 115585 125154 127166 122030 117460  96549 90305   NA
[10,] 110260 117529 118194 113802 106556  89727 83799   NA
[11,] 103444 109459 109940 105791  97347  82336 80320   NA
[12,] 102303 108483 106683 103076  92657  82653 82497   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/yougetitorg/rcomp/tmp/19d2c1304973564.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/yougetitorg/rcomp/tmp/237mj1304973564.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/yougetitorg/rcomp/tmp/3tem01304973564.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/yougetitorg/rcomp/tmp/4awn41304973564.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,]  82303.0  79596.0  74472.0 73562.0 66618 69029  89899.0  93774.0  90305
[2,]  91312.5  88299.5  84389.5 82026.5 77505 76531 102319.5 109076.5 106067
[3,]  98658.0  96945.0  92497.0 90687.0 87239 89235 113228.0 118711.0 117460
[4,] 102115.5  98938.0  95355.0 91849.5 88770 89668 117258.5 124010.5 123592
[5,] 107137.0 104699.0 100804.0 96066.0 91971 93228 120144.0 127233.0 127166
        [,10]    [,11]    [,12]
[1,]  83799.0  80320.0  82497.0
[2,]  98141.5  89841.5  87655.0
[3,] 110260.0 103444.0 102303.0
[4,] 115665.5 107625.0 104879.5
[5,] 118194.0 109940.0 108483.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,]  92206.62  90591.86 85948.58 84820.86 80511.72 81389.8 104306.7 109792.7
[2,] 105109.38 103298.14 99045.42 96553.14 93966.28 97080.2 122149.3 127629.3
         [,9]     [,10]     [,11]    [,12]
[1,] 106994.4  99794.95  92823.98  92016.8
[2,] 127925.6 120725.05 114064.02 112589.2

$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(82303, 91312.5, 98658, 102115.5, 107137, 79596,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/5mrux1304973564.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,]  79054.0  87239.0  91971  88744  88796  73993.0 66618.0   NA
[2,]  83784.5  92601.5  98435  94603  91592  82014.5 74017.0   NA
[3,]  95905.0 104940.0 106910 102955  97146  86303.5 81311.5   NA
[4,] 108905.5 118088.0 119169 114836 109892  94833.5 86849.0   NA
[5,] 115759.0 125620.0 127233 122401 118711 102394.0 93774.0   NA

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

$conf
          [,1]      [,2]      [,3]     [,4]      [,5]     [,6]     [,7] [,8]
[1,]  84447.14  93315.44  97453.08  93726.6  88799.25 80456.67 75458.74   NA
[2,] 107362.86 116564.56 116366.92 112183.4 105492.75 92150.33 87164.26   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(79054, 83784.5, 95905, 108905.5, 115759, 87239,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/6w3na1304973564.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,]  7894.350 22504.0  9823.00
[2,]  8822.955 24968.5 10884.25
[3,] 10167.782 26159.0 14035.50
[4,] 12458.027 31852.0 17374.25
[5,] 14379.965 36861.0 17783.50

$n
[1] 12 12 12

$conf
          [,1]     [,2]     [,3]
[1,]  8509.802 23019.39 11075.37
[2,] 11825.763 29298.61 16995.63

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(7894.34970330649, 8822.95479322016, 10167.7824884760,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/19d2c1304973564.ps tmp/19d2c1304973564.png",intern=TRUE))
character(0)
> try(system("convert tmp/237mj1304973564.ps tmp/237mj1304973564.png",intern=TRUE))
character(0)
> try(system("convert tmp/3tem01304973564.ps tmp/3tem01304973564.png",intern=TRUE))
character(0)
> try(system("convert tmp/4awn41304973564.ps tmp/4awn41304973564.png",intern=TRUE))
character(0)
> try(system("convert tmp/5mrux1304973564.ps tmp/5mrux1304973564.png",intern=TRUE))
character(0)
> try(system("convert tmp/6w3na1304973564.ps tmp/6w3na1304973564.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.310   1.010   1.801