R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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.

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(47460,50104,61465,53726,39477,43895,31481,29896,33842,39120,33702,25094,51442,45594,52518,48564,41745,49585,32747,33379,35645,37034,35681,20972,58552,54955,65540,51570,51145,46641,35704,33253,35193,41668,34865,21210,56126,49231,59723,48103,47472,50497,40059,34149,36860,46356,36577,23872,57276,56389,57657,62300,48929,51168,39636,33213,38127,43291,30600,21956,48033,46148,50736,48114,38390,44112,36287,30333,35908,40005,35263,26591,49709,47840,64781,57802,48154,54353,39737,37732,37163,43782,40649,29412)
> 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,] 47460 51442 58552 56126 57276 48033 49709   NA
 [2,] 50104 45594 54955 49231 56389 46148 47840   NA
 [3,] 61465 52518 65540 59723 57657 50736 64781   NA
 [4,] 53726 48564 51570 48103 62300 48114 57802   NA
 [5,] 39477 41745 51145 47472 48929 38390 48154   NA
 [6,] 43895 49585 46641 50497 51168 44112 54353   NA
 [7,] 31481 32747 35704 40059 39636 36287 39737   NA
 [8,] 29896 33379 33253 34149 33213 30333 37732   NA
 [9,] 33842 35645 35193 36860 38127 35908 37163   NA
[10,] 39120 37034 41668 46356 43291 40005 43782   NA
[11,] 33702 35681 34865 36577 30600 35263 40649   NA
[12,] 25094 20972 21210 23872 21956 26591 29412   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/wessaorg/rcomp/tmp/13o781322593466.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/wessaorg/rcomp/tmp/22flq1322593466.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/wessaorg/rcomp/tmp/3bjch1322593466.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/wessaorg/rcomp/tmp/4nxhu1322593466.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]   [,10]
[1,] 47460 45594.0 50736.0 48103 38390.0 43895.0 31481.0 29896 33842.0 37034.0
[2,] 48871 46994.0 55087.5 48339 40611.0 45376.5 34225.5 31773 35419.0 39562.5
[3,] 51442 49231.0 59723.0 51570 47472.0 49585.0 36287.0 33253 35908.0 41668.0
[4,] 56701 52529.5 63123.0 55764 48541.5 50832.5 39686.5 33764 37011.5 43536.5
[5,] 58552 56389.0 65540.0 62300 51145.0 54353.0 40059.0 34149 38127.0 46356.0
       [,11]   [,12]
[1,] 33702.0 20972.0
[2,] 34283.5 21583.0
[3,] 35263.0 23872.0
[4,] 36129.0 25842.5
[5,] 36577.0 29412.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,] 46766.05 45925.29 54924.33 47135.91 42736.03 46326.76 33025.78 32064.01
[2,] 56117.95 52536.71 64521.67 56004.09 52207.97 52843.24 39548.22 34441.99
         [,9]    [,10]   [,11]   [,12]
[1,] 34956.98 39294.79 34160.9 21328.3
[2,] 36859.02 44041.21 36365.1 26415.7

$out
[1] 37732 30600 40649

$group
[1]  8 11 11

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

Warning message:
In bxp(list(stats = c(47460, 48871, 51442, 56701, 58552, 45594,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5anhw1322593466.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,] 25094.0 20972.0 21210.0 23872.0 21956.0 26591.0 29412.0   NA
[2,] 32591.5 34512.0 35029.0 36718.5 35670.0 35585.5 38734.5   NA
[3,] 39298.5 39389.5 44154.5 46914.0 46110.0 39197.5 45811.0   NA
[4,] 48782.0 49074.5 53262.5 49864.0 56832.5 47090.5 52031.0   NA
[5,] 61465.0 52518.0 65540.0 59723.0 62300.0 50736.0 64781.0   NA

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7] [,8]
[1,] 31913.9 32747.45 35838.08 40918.25 36457.64 33949.99 39746.38   NA
[2,] 46683.1 46031.55 52470.92 52909.75 55762.36 44445.01 51875.62   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(25094, 32591.5, 39298.5, 48782, 61465, 20972,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6t9zn1322593466.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,] 1413.105  4285.0 1592.5
[2,] 3073.632  8509.0 2982.5
[3,] 3683.747 10253.5 5458.5
[4,] 4830.826 11923.5 7627.5
[5,] 5698.774 14804.0 8035.5

$n
[1] 12 12 12

$conf
         [,1]      [,2]     [,3]
[1,] 2882.279  8696.124 3339.884
[2,] 4485.214 11810.876 7577.116

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1413.10495684065, 3073.63206224009, 3683.74671408204,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/13o781322593466.ps tmp/13o781322593466.png",intern=TRUE))
character(0)
> try(system("convert tmp/22flq1322593466.ps tmp/22flq1322593466.png",intern=TRUE))
character(0)
> try(system("convert tmp/3bjch1322593466.ps tmp/3bjch1322593466.png",intern=TRUE))
character(0)
> try(system("convert tmp/4nxhu1322593466.ps tmp/4nxhu1322593466.png",intern=TRUE))
character(0)
> try(system("convert tmp/5anhw1322593466.ps tmp/5anhw1322593466.png",intern=TRUE))
character(0)
> try(system("convert tmp/6t9zn1322593466.ps tmp/6t9zn1322593466.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.409   0.321   1.735