R version 2.6.2 (2008-02-08)
Copyright (C) 2008 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(3722.7,3707.1,3697.4,3691.4,3679.4,3660.6,3647.5,3639,3623.5,3612.2,3604.2,3596.9,3603.9,3606.9,3607.7,3611.4,3611.9,3618.4,3629.5,3640.6,3646.8,3647.3,3649.5,3649.7,3652.5,3665.8,3674.3,3705.8,3735.8,3735.4,3747.7,3765.5,3779.7,3795.5,3813.1,3826.9,3833.3,3844.8,3851.3,3851.8,3854.1,3858.4,3861.6,3856.3,3855.8,3860.4,3855.1,3839.5,3833,3833.6,3826.8,3818.2,3811.4,3806.8,3810.3,3818.2,3858.7,3868.7,3872.6,3872.4,3876.4,3883.2,3883.4,3881.8,3887.1,3893.3,3902.8,3914.7,3929.8,3947.3,3971.8,3991.3,3993,3997.2,4016.2,4041.7,4060.5,4076.7,4102.7,4125.8,4140.1,4146.7,4157.9,4154.7,4145.3,4149,4142.4,4141,4146.1,4147.6,4142.7,4148.1,4159,4167,4178.5,4192.5,4209.2,4223.7,4231.6,4236.9,4253,4270.5,4285.7,4302.8,4321.5,4340.3,4356.5)
> 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] 107
> (np <- floor(n / par1))
[1] 8
> 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 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,] 3722.7 3603.9 3652.5 3833.3 3833.0 3876.4 3993.0 4145.3 4209.2
 [2,] 3707.1 3606.9 3665.8 3844.8 3833.6 3883.2 3997.2 4149.0 4223.7
 [3,] 3697.4 3607.7 3674.3 3851.3 3826.8 3883.4 4016.2 4142.4 4231.6
 [4,] 3691.4 3611.4 3705.8 3851.8 3818.2 3881.8 4041.7 4141.0 4236.9
 [5,] 3679.4 3611.9 3735.8 3854.1 3811.4 3887.1 4060.5 4146.1 4253.0
 [6,] 3660.6 3618.4 3735.4 3858.4 3806.8 3893.3 4076.7 4147.6 4270.5
 [7,] 3647.5 3629.5 3747.7 3861.6 3810.3 3902.8 4102.7 4142.7 4285.7
 [8,] 3639.0 3640.6 3765.5 3856.3 3818.2 3914.7 4125.8 4148.1 4302.8
 [9,] 3623.5 3646.8 3779.7 3855.8 3858.7 3929.8 4140.1 4159.0 4321.5
[10,] 3612.2 3647.3 3795.5 3860.4 3868.7 3947.3 4146.7 4167.0 4340.3
[11,] 3604.2 3649.5 3813.1 3855.1 3872.6 3971.8 4157.9 4178.5 4356.5
[12,] 3596.9 3649.7 3826.9 3839.5 3872.4 3991.3 4154.7 4192.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/html/rcomp/tmp/153di1210335280.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/2t6lf1210335280.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/3tv7e1210335280.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/4bh281210335280.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]
[1,] 3603.9 3606.9 3607.7 3611.4 3611.9 3618.4 3629.5 3639.0 3623.5 3612.2
[2,] 3722.7 3707.1 3697.4 3705.8 3735.8 3735.4 3747.7 3765.5 3779.7 3795.5
[3,] 3833.3 3844.8 3851.3 3851.8 3854.1 3858.4 3861.6 3856.3 3858.7 3868.7
[4,] 3993.0 3997.2 4016.2 4041.7 4060.5 4076.7 4102.7 4125.8 4140.1 4146.7
[5,] 4209.2 4223.7 4231.6 4236.9 4253.0 4270.5 4285.7 4302.8 4321.5 4340.3
      [,11]   [,12]
[1,] 3604.2 3596.90
[2,] 3813.1 3738.30
[3,] 3872.6 3855.95
[4,] 4157.9 4073.00
[5,] 4356.5 4192.50

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3690.942 3692.014 3683.399 3674.893 3683.091 3678.649 3674.633 3666.542
[2,] 3975.658 3997.586 4019.201 4028.707 4025.109 4038.151 4048.567 4046.058
         [,9]    [,10]    [,11]    [,12]
[1,] 3668.889 3683.735 3691.005 3668.982
[2,] 4048.511 4053.665 4054.195 4042.918

$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(3603.9, 3722.7, 3833.3, 3993, 4209.2, 3606.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5x97n1210335280.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]
[1,] 3596.90 3603.90 3652.50 3839.50 3806.8 3876.40 3993.00 4141.00 4209.20
[2,] 3617.85 3609.55 3690.05 3848.05 3814.8 3883.30 4028.95 4144.00 4234.25
[3,] 3654.05 3623.95 3741.75 3854.60 3829.9 3898.05 4089.70 4147.85 4270.50
[4,] 3694.40 3647.05 3787.60 3857.35 3863.7 3938.55 4143.40 4163.00 4312.15
[5,] 3722.70 3649.70 3826.90 3861.60 3872.6 3991.30 4157.90 4178.50 4356.50

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]     [,7]     [,8]
[1,] 3619.135 3606.846 3697.257 3850.358 3807.596 3872.85 4037.499 4139.184
[2,] 3688.965 3641.054 3786.243 3858.842 3852.204 3923.25 4141.901 4156.516
         [,9]
[1,] 4233.389
[2,] 4307.611

$out
[1] 3833.3 4192.5

$group
[1] 4 8

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

Warning message:
In bxp(list(stats = c(3596.9, 3617.85, 3654.05, 3694.4, 3722.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ytud1210335280.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,] 208.6785 595.60 249.55
[2,] 213.3667 620.35 304.45
[3,] 221.2086 646.60 338.60
[4,] 236.1742 680.90 353.10
[5,] 250.9612 752.30 360.40

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 210.8060 618.9827 316.4104
[2,] 231.6112 674.2173 360.7896

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(208.678496736008, 213.366728571026, 221.208583921710,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/153di1210335280.ps tmp/153di1210335280.png")
> system("convert tmp/2t6lf1210335280.ps tmp/2t6lf1210335280.png")
> system("convert tmp/3tv7e1210335280.ps tmp/3tv7e1210335280.png")
> system("convert tmp/4bh281210335280.ps tmp/4bh281210335280.png")
> system("convert tmp/5x97n1210335280.ps tmp/5x97n1210335280.png")
> system("convert tmp/6ytud1210335280.ps tmp/6ytud1210335280.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.273   0.800   1.563