R version 2.9.0 (2009-04-17)
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(11310,64305,15310,37299,21302,61308,72300,26303,18301,54305,66309,50301,31298,52291,87286,81288,14293,90302,50306,15310,44310,26314,98313,76310,25313,48309,95307,10320,87327,63328,34333,90333,81332,7342,30424,13344,88347,40339,23330,1339,10341,46342,81342,2342,76350,35368,93367,88377,39376,41366,77375,56382,79397,26385,73397,28404,98413,73414,47423,52431,24441,92439,90441,441,13448,18458,18459,69477,41491,10492,73508,82515,13525,55533,19550,85558,57563,60570,49568,51570,26561,61558,78548,77537,539,18540,47542,86542,81544,16543,22538,25538,99527,63518,95508,65496,5488,96475,81465,5463,81458,74445,21434,67427,27418,81407,82395,97359)
> 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] 108
> (np <- floor(n / par1))
[1] 9
> 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] 9 9 9 9 9 9 9 9 9 9 9 9
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10]
 [1,] 11310 31298 25313 88347 39376 24441 13525   539  5488    NA
 [2,] 64305 52291 48309 40339 41366 92439 55533 18540 96475    NA
 [3,] 15310 87286 95307 23330 77375 90441 19550 47542 81465    NA
 [4,] 37299 81288 10320  1339 56382   441 85558 86542  5463    NA
 [5,] 21302 14293 87327 10341 79397 13448 57563 81544 81458    NA
 [6,] 61308 90302 63328 46342 26385 18458 60570 16543 74445    NA
 [7,] 72300 50306 34333 81342 73397 18459 49568 22538 21434    NA
 [8,] 26303 15310 90333  2342 28404 69477 51570 25538 67427    NA
 [9,] 18301 44310 81332 76350 98413 41491 26561 99527 27418    NA
[10,] 54305 26314  7342 35368 73414 10492 61558 63518 81407    NA
[11,] 66309 98313 30424 93367 47423 73508 78548 95508 82395    NA
[12,] 50301 76310 13344 88377 52431 82515 77537 65496 97359    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/1p2n91244310470.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/2ik371244310470.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/3w52e1244310470.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/4etoz1244310470.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] [,11] [,12]
[1,]   539 18540 15310   441 10341 16543 18459  2342 18301  7342 30424 13344
[2,] 11310 41366 23330  5463 14293 26385 22538 25538 27418 26314 66309 52431
[3,] 24441 52291 77375 37299 57563 60570 49568 28404 44310 54305 78548 76310
[4,] 31298 64305 87286 81288 81458 63328 72300 67427 81332 63518 93367 82515
[5,] 39376 96475 95307 86542 87327 90302 81342 90333 99527 81407 98313 97359

$n
 [1] 9 9 9 9 9 9 9 9 9 9 9 9

$conf
         [,1]     [,2]      [,3]    [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 13913.99 40209.79  43691.51 -2635.5 22189.43 41113.35 23360.01  6342.46
[2,] 34968.01 64372.21 111058.49 77233.5 92936.57 80026.65 75775.99 50465.54
         [,9]    [,10]    [,11]    [,12]
[1,] 15915.29 34710.89 64297.45 60465.76
[2,] 72704.71 73899.11 92798.55 92154.24

$out
[1] 88347

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(539, 11310, 24441, 31298, 39376, 18540, 41366,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5u1d61244310470.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]  [,8]  [,9] [,10]
[1,] 11310.0 14293.0  7342.0  1339.0 26385.0   441.0 13525.0   539  5463    NA
[2,] 19801.5 28806.0 19328.5 16835.5 40371.0 15953.0 38064.5 20539 24426    NA
[3,] 43800.0 51298.5 41321.0 43340.5 54406.5 32966.0 56548.0 55530 77926    NA
[4,] 62806.5 84287.0 84329.5 84844.5 75394.5 78011.5 69547.5 84043 81930    NA
[5,] 72300.0 98313.0 95307.0 93367.0 98413.0 92439.0 85558.0 99527 97359    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]      [,6]     [,7]    [,8]
[1,] 24185.13 25993.24 11673.61 12321.14 38432.05  4660.702 42188.39 26565.4
[2,] 63414.87 76603.76 70968.39 74359.86 70380.95 61271.298 70907.61 84494.6
          [,9] [,10]
[1,]  51698.04    NA
[2,] 104153.96    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(11310, 19801.5, 43800, 62806.5, 72300, 14293,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6p4931244310470.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,] 22870.70 62883 19988.0
[2,] 25193.57 73912 28571.0
[3,] 26704.68 78966 39546.5
[4,] 32594.25 85058 58935.0
[5,] 37671.22 87991 75825.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 23329.19 73882.23 25697.28
[2,] 30080.18 84049.77 53395.72

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(22870.6980423520, 25193.574665704, 26704.6825922424,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1p2n91244310470.ps tmp/1p2n91244310470.png")
> system("convert tmp/2ik371244310470.ps tmp/2ik371244310470.png")
> system("convert tmp/3w52e1244310470.ps tmp/3w52e1244310470.png")
> system("convert tmp/4etoz1244310470.ps tmp/4etoz1244310470.png")
> system("convert tmp/5u1d61244310470.ps tmp/5u1d61244310470.png")
> system("convert tmp/6p4931244310470.ps tmp/6p4931244310470.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.073   0.806   1.346