R version 2.12.2 (2011-02-25)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-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(3893.9,3799.2,3769.6,3768.6,3854.9,3778.5,3779.7,3803.2,3900.3,3792.6,3767.4,3752.6,3829.6,3722.6,3692.9,3681,3762.9,3661.7,3633.1,3621.5,3710,3619.4,3595.2,3573.2,3650.1,3554.2,3537,3528.6,3597.1,3521.9,3516.5,3515.7,3600.2,3517.1,3513.7,3528.2,3608.3,3502.5,3502.5,3495.3,3543.8,3425.3,3418.4,3406.4,3446.1,3341.1,3347,3354.9,3399,3288.9,3279,3275.2,3314,3227.1,3225.3,3228.6,3287.1,3210.1,3213.1,3228,3287,3211,3199.8,3166.3,3164,3156.7,3156,3165.5,3179.2,3182.5,3179.5,3193.5,3219.6,3221.9,3210.1)
> 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] 75
> (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 7 7 6 6 6 6 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 3893.9 3829.6 3650.1 3608.3 3399.0 3287.0 3219.6
 [2,] 3799.2 3722.6 3554.2 3502.5 3288.9 3211.0 3221.9
 [3,] 3769.6 3692.9 3537.0 3502.5 3279.0 3199.8 3210.1
 [4,] 3768.6 3681.0 3528.6 3495.3 3275.2 3166.3     NA
 [5,] 3854.9 3762.9 3597.1 3543.8 3314.0 3164.0     NA
 [6,] 3778.5 3661.7 3521.9 3425.3 3227.1 3156.7     NA
 [7,] 3779.7 3633.1 3516.5 3418.4 3225.3 3156.0     NA
 [8,] 3803.2 3621.5 3515.7 3406.4 3228.6 3165.5     NA
 [9,] 3900.3 3710.0 3600.2 3446.1 3287.1 3179.2     NA
[10,] 3792.6 3619.4 3517.1 3341.1 3210.1 3182.5     NA
[11,] 3767.4 3595.2 3513.7 3347.0 3213.1 3179.5     NA
[12,] 3752.6 3573.2 3528.2 3354.9 3228.0 3193.5     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/wessaorg/rcomp/tmp/1qx281305809226.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/wessaorg/rcomp/tmp/2fcop1305809226.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/wessaorg/rcomp/tmp/353lh1305809226.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/wessaorg/rcomp/tmp/4413x1305809226.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,] 3219.60 3211.0 3199.80 3166.30 3164.00 3156.7 3156.00 3165.50 3179.20
[2,] 3343.00 3255.4 3244.55 3275.20 3314.00 3227.1 3225.30 3228.60 3287.10
[3,] 3608.30 3502.5 3502.50 3511.95 3570.45 3473.6 3467.45 3461.05 3523.15
[4,] 3739.85 3638.4 3614.95 3681.00 3762.90 3661.7 3633.10 3621.50 3710.00
[5,] 3893.90 3799.2 3769.60 3768.60 3854.90 3778.5 3779.70 3803.20 3900.30
      [,10]   [,11]   [,12]
[1,] 3182.5 3179.50 3193.50
[2,] 3210.1 3213.10 3228.00
[3,] 3429.1 3430.35 3441.55
[4,] 3619.4 3595.20 3573.20
[5,] 3792.6 3767.40 3752.60

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3371.308 3273.779 3281.303 3250.196 3280.895 3193.269 3204.406 3207.617
[2,] 3845.292 3731.221 3723.697 3773.704 3860.005 3753.931 3730.494 3714.483
         [,9]    [,10]    [,11]    [,12]
[1,] 3250.366 3165.088 3183.883 3218.885
[2,] 3795.934 3693.112 3676.817 3664.215

$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(3219.6, 3343, 3608.3, 3739.85, 3893.9, 3211,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/5qa0v1305809226.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]
[1,] 3752.60 3573.20 3513.70 3341.10 3210.1 3156.00 3210.10
[2,] 3769.10 3620.45 3516.80 3380.65 3226.2 3164.75 3214.85
[3,] 3786.15 3671.35 3528.40 3435.70 3251.9 3179.35 3219.60
[4,] 3829.05 3716.30 3575.65 3502.50 3288.0 3196.65 3220.75
[5,] 3900.30 3829.60 3650.10 3608.30 3314.0 3211.00 3221.90

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]   [,6]     [,7]
[1,] 3758.806 3627.632 3501.558 3380.123 3223.713 3164.8 3214.218
[2,] 3813.494 3715.068 3555.242 3491.277 3280.087 3193.9 3224.982

$out
[1] 3399 3287

$group
[1] 5 6

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

Warning message:
In bxp(list(stats = c(3752.6, 3769.1, 3786.15, 3829.05, 3900.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/6aftr1305809226.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,] 217.3353 559.10 302.225
[2,] 231.0021 588.05 325.125
[3,] 239.7566 615.95 350.050
[4,] 251.9442 656.00 363.050
[5,] 269.4465 721.10 396.850

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 230.2048 584.9575 332.7522
[2,] 249.3084 646.9425 367.3478

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(217.335289357251, 231.002140387832, 239.75658727635,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1qx281305809226.ps tmp/1qx281305809226.png",intern=TRUE))
character(0)
> try(system("convert tmp/2fcop1305809226.ps tmp/2fcop1305809226.png",intern=TRUE))
character(0)
> try(system("convert tmp/353lh1305809226.ps tmp/353lh1305809226.png",intern=TRUE))
character(0)
> try(system("convert tmp/4413x1305809226.ps tmp/4413x1305809226.png",intern=TRUE))
character(0)
> try(system("convert tmp/5qa0v1305809226.ps tmp/5qa0v1305809226.png",intern=TRUE))
character(0)
> try(system("convert tmp/6aftr1305809226.ps tmp/6aftr1305809226.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  4.030   0.370   4.531