R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)

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(32819,32700,32242,32810,33865,32226,31077,31293,30236,30160,32436,30695,27525,26434,25739,25204,24977,24320,22680,22052,21467,21383,21777,21928,21814,22937,23595,20830,19650,19195,19644,18483,18079,19178,18391,18441,18584,20108,20148,19394,17745,17696,17032,16438,15683,15594,15713,15937,16171,15928,16348,15579,15305,15648,14954,15137,15839,16050,15168,17064,16005,14886,14931,14544,13812,13031,12574,11964,11451,11346,11353,10702,10646,10556,10463,10407,10625,10872,10805,10653,10574,10431,10383,10296,10872,10635,10297,10570,10662,10709,10413,10846,10371,9924,9828)
> 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] 95
> (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] 8 8 8 8 8 8 8 8 8 8 8 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 32819 27525 21814 18584 16171 16005 10646 10872
 [2,] 32700 26434 22937 20108 15928 14886 10556 10635
 [3,] 32242 25739 23595 20148 16348 14931 10463 10297
 [4,] 32810 25204 20830 19394 15579 14544 10407 10570
 [5,] 33865 24977 19650 17745 15305 13812 10625 10662
 [6,] 32226 24320 19195 17696 15648 13031 10872 10709
 [7,] 31077 22680 19644 17032 14954 12574 10805 10413
 [8,] 31293 22052 18483 16438 15137 11964 10653 10846
 [9,] 30236 21467 18079 15683 15839 11451 10574 10371
[10,] 30160 21383 19178 15594 16050 11346 10431  9924
[11,] 32436 21777 18391 15713 15168 11353 10383  9828
[12,] 30695 21928 18441 15937 17064 10702 10296    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/10z6q1305814441.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/2ca4e1305814441.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/30jd41305814441.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/4l4d31305814441.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,] 10646.0 10556.0 10297 10407.0 10625.0 10709.0 10413.0 10653.0 10371.0
[2,] 13438.5 12760.5 12697 12557.0 12237.0 11951.5 11689.5 11405.0 11012.5
[3,] 17377.5 18018.0 18248 17486.5 16525.0 16672.0 15993.0 15787.5 15761.0
[4,] 24669.5 24685.5 24667 23017.0 22313.5 21757.5 21162.0 20267.5 19773.0
[5,] 32819.0 32700.0 32242 32810.0 33865.0 32226.0 31077.0 31293.0 30236.0
       [,10]   [,11]   [,12]
[1,]  9924.0  9828.0 10296.0
[2,] 10888.5 10868.0 13319.5
[3,] 15822.0 15440.5 17064.0
[4,] 20280.5 20084.0 20184.5
[5,] 30160.0 32436.0 21928.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 11103.70 11356.52 11561.39 11643.39 10896.12 11194.23 10701.53 10836.78
[2,] 23651.30 24679.48 24934.61 23329.61 22153.88 22149.77 21284.47 20738.22
         [,9]    [,10]    [,11]    [,12]
[1,] 10867.26 10575.49 10292.31 12964.33
[2,] 20654.74 21068.51 20588.69 21163.67

$out
[1] 30695

$group
[1] 12

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

Warning message:
In bxp(list(stats = c(10646, 13438.5, 17377.5, 24669.5, 32819, 10556,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5jix61305814441.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,] 30160 21383.0 18079.0 15594 14954.0 10702.0 10296.0  9828.0
[2,] 30886 21852.5 18462.0 15825 15236.5 11402.0 10419.0 10334.0
[3,] 32234 23500.0 19419.5 17364 15743.5 12802.5 10565.0 10570.0
[4,] 32755 25471.5 21322.0 18989 16110.5 14715.0 10649.5 10685.5
[5,] 33865 27525.0 23595.0 20148 17064.0 16005.0 10872.0 10872.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 31381.54 21849.35 18115.03 15920.88 15344.86 11291.42 10459.87 10402.55
[2,] 33086.46 25150.65 20723.97 18807.12 16142.14 14313.58 10670.13 10737.45

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(30160, 30886, 32234, 32755, 33865, 21383,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6zjd31305814441.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,] 6705.645 19865.0 6865.000
[2,] 6974.007 20519.5 7825.750
[3,] 7430.448 21731.0 8199.125
[4,] 7740.500 22288.0 8566.000
[5,] 7886.079 23240.0 8612.000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 7080.846 20924.38 7861.492
[2,] 7780.051 22537.62 8536.758

$out
[1]  9988 10317

$group
[1] 3 3

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

Warning message:
In bxp(list(stats = c(6705.64478629759, 6974.00738543284, 7430.44830288666,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/10z6q1305814441.ps tmp/10z6q1305814441.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ca4e1305814441.ps tmp/2ca4e1305814441.png",intern=TRUE))
character(0)
> try(system("convert tmp/30jd41305814441.ps tmp/30jd41305814441.png",intern=TRUE))
character(0)
> try(system("convert tmp/4l4d31305814441.ps tmp/4l4d31305814441.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jix61305814441.ps tmp/5jix61305814441.png",intern=TRUE))
character(0)
> try(system("convert tmp/6zjd31305814441.ps tmp/6zjd31305814441.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.626   0.168   1.801