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(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.9,3867.8,3872.3,3873.3,3876.7,3882.6,3883.5,3882.2,3888.1,3893.7,3901.9,3914.3,3930.3,3948.3,3971.5,3990.1,3993,3998,4015.8,4041.2,4060.7,4076.7,4103,4125.3,4139.7,4146.7,4158,4155.1,4144.8,4148.2,4142.5,4142.1,4145.4,4146.3,4143.5,4149.2,4158.9,4166.1,4179.1,4194.4,4211.7,4226.3,4235.8,4243.6,4258.7,4278.2,4298,4315.1,4334.3,4356,4374,4395.5,4417.8,4432.8,4446.3)
> 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] 73
> (np <- floor(n / par1))
[1] 6
> arr <- array(NA,dim=c(par1,np+1))
> darr <- array(NA,dim=c(par1,np+1))
> ari <- array(0,dim=par1)
> dx <- diff(x)
> j <- 0
> for (i in 1:n)
+ {
+ j = j + 1
+ ari[j] = ari[j] + 1
+ arr[j,ari[j]] <- x[i]
+ darr[j,ari[j]] <- dx[i]
+ if (j == par1) j = 0
+ }
> ari
 [1] 7 6 6 6 6 6 6 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 3851.3 3826.8 3883.5 4015.8 4142.5 4235.8 4446.3
 [2,] 3851.8 3818.2 3882.2 4041.2 4142.1 4243.6     NA
 [3,] 3854.1 3811.4 3888.1 4060.7 4145.4 4258.7     NA
 [4,] 3858.4 3806.8 3893.7 4076.7 4146.3 4278.2     NA
 [5,] 3861.6 3810.3 3901.9 4103.0 4143.5 4298.0     NA
 [6,] 3856.3 3818.2 3914.3 4125.3 4149.2 4315.1     NA
 [7,] 3855.8 3858.9 3930.3 4139.7 4158.9 4334.3     NA
 [8,] 3860.4 3867.8 3948.3 4146.7 4166.1 4356.0     NA
 [9,] 3855.1 3872.3 3971.5 4158.0 4179.1 4374.0     NA
[10,] 3839.5 3873.3 3990.1 4155.1 4194.4 4395.5     NA
[11,] 3833.0 3876.7 3993.0 4144.8 4211.7 4417.8     NA
[12,] 3833.6 3882.6 3998.0 4148.2 4226.3 4432.8     NA
> darr
       [,1] [,2] [,3]  [,4] [,5] [,6] [,7]
 [1,]   0.5 -8.6 -1.3  25.4 -0.4  7.8   NA
 [2,]   2.3 -6.8  5.9  19.5  3.3 15.1   NA
 [3,]   4.3 -4.6  5.6  16.0  0.9 19.5   NA
 [4,]   3.2  3.5  8.2  26.3 -2.8 19.8   NA
 [5,]  -5.3  7.9 12.4  22.3  5.7 17.1   NA
 [6,]  -0.5 40.7 16.0  14.4  9.7 19.2   NA
 [7,]   4.6  8.9 18.0   7.0  7.2 21.7   NA
 [8,]  -5.3  4.5 23.2  11.3 13.0 18.0   NA
 [9,] -15.6  1.0 18.6  -2.9 15.3 21.5   NA
[10,]  -6.5  3.4  2.9 -10.3 17.3 22.3   NA
[11,]   0.6  5.9  5.0   3.4 14.6 15.0   NA
[12,]  -6.8  0.9 17.8  -5.7  9.5 13.5   NA
> arr.mean <- array(NA,dim=par1)
> arr.median <- array(NA,dim=par1)
> arr.midrange <- array(NA,dim=par1)
> for (j in 1:par1)
+ {
+ arr.mean[j] <- mean(arr[j,],na.rm=TRUE)
+ arr.median[j] <- median(arr[j,],na.rm=TRUE)
+ arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2
+ }
> overall.mean <- mean(x)
> overall.median <- median(x)
> overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2
> postscript(file="/var/www/html/rcomp/tmp/1ngmu1242325091.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.mean,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2w1g61242325091.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.median,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/3do8f1242325091.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.midrange,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/4oiy71242325091.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,] 3826.80 3818.2 3811.4 3806.8 3810.30 3818.2 3855.8 3860.4 3855.10 3839.5
[2,] 3867.40 3851.8 3854.1 3858.4 3861.60 3856.3 3858.9 3867.8 3872.30 3873.3
[3,] 4015.80 3961.7 3974.4 3985.2 4002.45 4019.8 4035.0 4047.5 4064.75 4072.6
[4,] 4189.15 4142.1 4145.4 4146.3 4143.50 4149.2 4158.9 4166.1 4179.10 4194.4
[5,] 4446.30 4243.6 4258.7 4278.2 4298.00 4315.1 4334.3 4356.0 4374.00 4395.5
      [,11]  [,12]
[1,] 3833.0 3833.6
[2,] 3876.7 3882.6
[3,] 4068.9 4073.1
[4,] 4211.7 4226.3
[5,] 4417.8 4432.8

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]    [,7]     [,8]
[1,] 3823.656 3774.447 3786.502 3799.495 3820.615 3830.87 3841.49 3855.087
[2,] 4207.944 4148.953 4162.298 4170.905 4184.285 4208.73 4228.51 4239.913
         [,9]   [,10]    [,11]    [,12]
[1,] 3866.854 3865.48 3852.814 3851.402
[2,] 4262.646 4279.72 4284.986 4294.798

$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(3826.8, 3867.4, 4015.8, 4189.15, 4446.3, 3818.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5u6iu1242325091.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> z <- data.frame(t(darr))
> names(z) <- c(1:par1)
> (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries'))
$stats
      [,1] [,2]  [,3]  [,4]  [,5] [,6]  [,7]  [,8]   [,9]  [,10] [,11] [,12]
[1,] -8.60 -6.8 -4.60 -2.80 -5.30 -0.5  4.60 -5.30 -15.60 -10.30  0.60  -6.8
[2,] -1.30  2.3  0.90  3.20  5.70  9.7  7.00  4.50  -2.90  -6.50  3.40  -5.7
[3,]  0.05  4.6  4.95  5.85 10.15 15.2  8.05 12.15   8.15   3.15  5.45   5.2
[4,]  7.80 15.1 16.00 19.80 17.10 19.2 18.00 18.00  18.60  17.30 14.60  13.5
[5,]  7.80 19.5 19.50 26.30 22.30 19.2 21.70 23.20  21.50  22.30 15.00  17.8

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]       [,7]
[1,] -5.819794 -3.656413 -4.789988 -4.857536  2.796632  9.072193  0.9546447
[2,]  5.919794 12.856413 14.689988 16.557536 17.503368 21.327807 15.1453553
          [,8]      [,9]     [,10]     [,11]    [,12]
[1,]  3.442064 -5.718194 -12.20177 -1.774362 -7.18462
[2,] 20.857936 22.018194  18.50177 12.674362 17.58462

$out
[1] 25.4 40.7

$group
[1] 1 6

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

Warning message:
In bxp(list(stats = c(-8.60000000000036, -1.30000000000018, 0.0500000000001819,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6azhr1242325091.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]
[1,] 3833.00 3806.80 3882.2 4015.80 4142.10 4235.80 4446.3
[2,] 3845.40 3814.80 3890.9 4068.70 4144.45 4268.45 4446.3
[3,] 3854.60 3842.85 3922.3 4132.50 4154.05 4324.70 4446.3
[4,] 3857.35 3872.80 3980.8 4147.45 4186.75 4384.75 4446.3
[5,] 3861.60 3882.60 3998.0 4158.00 4226.30 4432.80 4446.3

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]   [,7]
[1,] 3849.150 3816.396 3881.296 4096.582 4134.757 4271.655 4446.3
[2,] 3860.050 3869.304 3963.304 4168.418 4173.343 4377.745 4446.3

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(3833, 3845.4, 3854.6, 3857.35, 3861.6, 3806.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/795u11242325091.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> z <- data.frame(cbind(arr.mean,arr.median,arr.midrange))
> names(z) <- list('mean','median','midrange')
> (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots'))
$stats
         [,1]     [,2]     [,3]
[1,] 3996.517 3961.700 3988.138
[2,] 4014.867 3993.825 4000.294
[3,] 4051.873 4027.400 4020.006
[4,] 4071.492 4066.825 4039.500
[5,] 4086.917 4073.100 4059.113

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4026.046 3994.104 4002.124
[2,] 4077.700 4060.696 4037.888

$out
numeric(0)

$group
numeric(0)

$names
[1] "mean"     "median"   "midrange"

Warning message:
In bxp(list(stats = c(3996.51666666667, 4014.86666666667, 4051.87261904762,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ngmu1242325091.ps tmp/1ngmu1242325091.png")
> system("convert tmp/2w1g61242325091.ps tmp/2w1g61242325091.png")
> system("convert tmp/3do8f1242325091.ps tmp/3do8f1242325091.png")
> system("convert tmp/4oiy71242325091.ps tmp/4oiy71242325091.png")
> system("convert tmp/5u6iu1242325091.ps tmp/5u6iu1242325091.png")
> system("convert tmp/6azhr1242325091.ps tmp/6azhr1242325091.png")
> system("convert tmp/795u11242325091.ps tmp/795u11242325091.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.210   0.948   1.564