R version 2.7.0 (2008-04-22)
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 = '4'
> #'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] 26
> 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] 27 27 27 26
> arr
       [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
[1,] 3722.7 3679.4 3623.5 3603.9 3611.9 3646.8 3652.5 3735.8 3779.7 3833.3
[2,] 3707.1 3660.6 3612.2 3606.9 3618.4 3647.3 3665.8 3735.4 3795.5 3844.8
[3,] 3697.4 3647.5 3604.2 3607.7 3629.5 3649.5 3674.3 3747.7 3813.1 3851.3
[4,] 3691.4 3639.0 3596.9 3611.4 3640.6 3649.7 3705.8 3765.5 3826.9 3851.8
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 3854.1 3855.8 3833.0 3811.4 3858.7 3876.4 3887.1 3929.8 3993.0 4060.5
[2,] 3858.4 3860.4 3833.6 3806.8 3868.7 3883.2 3893.3 3947.3 3997.2 4076.7
[3,] 3861.6 3855.1 3826.8 3810.3 3872.6 3883.4 3902.8 3971.8 4016.2 4102.7
[4,] 3856.3 3839.5 3818.2 3818.2 3872.4 3881.8 3914.7 3991.3 4041.7 4125.8
      [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]
[1,] 4140.1 4145.3 4146.1 4159.0 4209.2 4253.0 4321.5
[2,] 4146.7 4149.0 4147.6 4167.0 4223.7 4270.5 4340.3
[3,] 4157.9 4142.4 4142.7 4178.5 4231.6 4285.7 4356.5
[4,] 4154.7 4141.0 4148.1 4192.5 4236.9 4302.8     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/1ybta1212317134.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/2x97g1212317134.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/38l7z1212317134.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/4ei6v1212317134.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]
[1,] 3603.90 3606.90 3604.20 3596.90
[2,] 3729.25 3721.25 3722.55 3705.80
[3,] 3855.80 3860.40 3861.60 3854.05
[4,] 4100.30 4111.70 4122.55 4125.80
[5,] 4321.50 4340.30 4356.50 4302.80

$n
[1] 27 27 27 26

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 3742.974 3741.675 3739.972 3723.907
[2,] 3968.626 3979.125 3983.228 3984.193

$out
numeric(0)

$group
numeric(0)

$names
[1] "1" "2" "3" "4"

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5s85l1212317134.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,] 3691.40 3639.00 3596.90 3603.90 3611.90 3646.80 3652.50 3735.40 3779.7
[2,] 3694.40 3643.25 3600.55 3605.40 3615.15 3647.05 3659.15 3735.60 3787.6
[3,] 3702.25 3654.05 3608.20 3607.30 3623.95 3648.40 3670.05 3741.75 3804.3
[4,] 3714.90 3670.00 3617.85 3609.55 3635.05 3649.60 3690.05 3756.60 3820.0
[5,] 3722.70 3679.40 3623.50 3611.40 3640.60 3649.70 3705.80 3765.50 3826.9
       [,10]   [,11]   [,12]  [,13]   [,14]   [,15]  [,16]   [,17]   [,18]
[1,] 3833.30 3854.10 3839.50 3818.2 3806.80 3858.70 3876.4 3887.10 3929.80
[2,] 3839.05 3855.20 3847.30 3822.5 3808.55 3863.70 3879.1 3890.20 3938.55
[3,] 3848.05 3857.35 3855.45 3829.9 3810.85 3870.55 3882.5 3898.05 3959.55
[4,] 3851.55 3860.00 3858.10 3833.3 3814.80 3872.50 3883.3 3908.75 3981.55
[5,] 3851.80 3861.60 3860.40 3833.6 3818.20 3872.60 3883.4 3914.70 3991.30
       [,19]   [,20]  [,21]   [,22]   [,23]   [,24]   [,25]   [,26]  [,27]
[1,] 3993.00 4060.50 4140.1 4141.00 4142.70 4159.00 4209.20 4253.00 4321.5
[2,] 3995.10 4068.60 4143.4 4141.70 4144.40 4163.00 4216.45 4261.75 4330.9
[3,] 4006.70 4089.70 4150.7 4143.85 4146.85 4172.75 4227.65 4278.10 4340.3
[4,] 4028.95 4114.25 4156.3 4147.15 4147.85 4185.50 4234.25 4294.25 4348.4
[5,] 4041.70 4125.80 4157.9 4149.00 4148.10 4192.50 4236.90 4302.80 4356.5

$n
 [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 3

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 3686.055 3632.918 3594.533 3604.022 3608.229 3646.386 3645.639 3725.16
[2,] 3718.445 3675.183 3621.867 3610.579 3639.671 3650.414 3694.461 3758.34
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 3778.704 3838.175 3853.558 3846.918 3821.368 3805.913 3863.598 3879.182
[2,] 3829.896 3857.925 3861.142 3863.982 3838.432 3815.788 3877.502 3885.818
        [,17]   [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 3883.396 3925.58 3979.958 4053.636 4140.509 4139.545 4144.124 4154.975
[2,] 3912.705 3993.52 4033.441 4125.764 4160.891 4148.155 4149.576 4190.525
        [,25]    [,26]    [,27]
[1,] 4213.588 4252.425 4324.336
[2,] 4241.712 4303.775 4356.264

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" NA  

Warning message:
In bxp(list(stats = c(3691.4, 3694.4, 3702.25, 3714.9, 3722.7, 3639,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6qzwq1212317134.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,] 212.1697 705.90 371.050
[2,] 213.3107 711.75 377.550
[3,] 217.0401 725.50 387.250
[4,] 221.7224 742.85 395.225
[5,] 223.8165 752.30 400.000

$n
[1] 4 4 4

$conf
         [,1]    [,2]     [,3]
[1,] 210.3949 700.931 373.2868
[2,] 223.6854 750.069 401.2132

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(212.169678432508, 213.310743807377, 217.040126823757,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/1ybta1212317134.ps tmp/1ybta1212317134.png")
> system("convert tmp/2x97g1212317134.ps tmp/2x97g1212317134.png")
> system("convert tmp/38l7z1212317134.ps tmp/38l7z1212317134.png")
> system("convert tmp/4ei6v1212317134.ps tmp/4ei6v1212317134.png")
> system("convert tmp/5s85l1212317134.ps tmp/5s85l1212317134.png")
> system("convert tmp/6qzwq1212317134.ps tmp/6qzwq1212317134.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.392   0.880   6.013