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 = '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))
> 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] 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
> darr
       [,1] [,2] [,3]  [,4] [,5] [,6] [,7] [,8] [,9]
 [1,] -15.6  3.0 13.3  11.5  0.6  6.8  4.2  3.7 14.5
 [2,]  -9.7  0.8  8.5   6.5 -6.8  0.2 19.0 -6.6  7.9
 [3,]  -6.0  3.7 31.5   0.5 -8.6 -1.6 25.5 -1.4  5.3
 [4,] -12.0  0.5 30.0   2.3 -6.8  5.3 18.8  5.1 16.1
 [5,] -18.8  6.5 -0.4   4.3 -4.6  6.2 16.2  1.5 17.5
 [6,] -13.1 11.1 12.3   3.2  3.5  9.5 26.0 -4.9 15.2
 [7,]  -8.5 11.1 17.8  -5.3  7.9 11.9 23.1  5.4 17.1
 [8,] -15.5  6.2 14.2  -0.5 40.5 15.1 14.3 10.9 18.7
 [9,] -11.3  0.5 15.8   4.6 10.0 17.5  6.6  8.0 18.8
[10,]  -8.0  2.2 17.6  -5.3  3.9 24.5 11.2 11.5 16.2
[11,]  -7.3  0.2 13.8 -15.6 -0.2 19.5 -3.2 14.0   NA
[12,]   7.0  2.8  6.4  -6.5  4.0  1.7 -9.4 16.7   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/freestat/rcomp/tmp/16iv81212178353.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/freestat/rcomp/tmp/2qxhj1212178353.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/freestat/rcomp/tmp/3cnf61212178353.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/freestat/rcomp/tmp/4myr71212178353.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/freestat/rcomp/tmp/58nxl1212178353.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,]  0.6 -9.7 -8.6 -12.0 -4.6 -4.9 -8.5 -0.5 -11.3  -8.0 -15.60  -9.4
[2,]  3.0 -6.6 -1.6   0.5 -0.4  3.2  5.4  6.2   4.6   2.2  -5.25  -2.4
[3,]  4.2  0.8  0.5   5.1  4.3  9.5 11.1 14.2   8.0  11.2   0.00   3.4
[4,] 11.5  7.9  5.3  16.1  6.5 12.3 17.1 15.1  15.8  16.2  13.90   6.7
[5,] 14.5 19.0  5.3  30.0 16.2 15.2 23.1 18.7  18.8  24.5  19.50  16.7

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

$conf
           [,1]      [,2]   [,3]   [,4]  [,5]      [,6]   [,7]      [,8]
[1,] -0.2766667 -6.836667 -3.134 -3.116 0.666  4.707333  4.938  9.512667
[2,]  8.6766667  8.436667  4.134 13.316 7.934 14.292667 17.262 18.887333
          [,9]     [,10]     [,11]     [,12]
[1,]  2.101333  3.826667 -10.69746 -1.683391
[2,] 13.898667 18.573333  10.69746  8.483391

$out
[1] -15.6  31.5  25.5 -18.8  17.5 -13.1  26.0 -15.5  40.5

$group
[1] 1 3 3 5 5 6 6 8 8

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

Warning message:
In bxp(list(stats = c(0.599999999999909, 3, 4.19999999999982, 11.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6x6df1212178354.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/freestat/rcomp/tmp/73jmd1212178354.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,] 3874.367 3844.800 3852.150
[2,] 3883.950 3851.550 3865.800
[3,] 3894.833 3856.125 3906.713
[4,] 3918.106 3860.150 3952.775
[5,] 3939.911 3872.600 3985.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 3879.255 3852.202 3867.043
[2,] 3910.412 3860.048 3946.382

$out
[1] 3833.3

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(3874.36666666667, 3883.95, 3894.83333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/16iv81212178353.ps tmp/16iv81212178353.png")
> system("convert tmp/2qxhj1212178353.ps tmp/2qxhj1212178353.png")
> system("convert tmp/3cnf61212178353.ps tmp/3cnf61212178353.png")
> system("convert tmp/4myr71212178353.ps tmp/4myr71212178353.png")
> system("convert tmp/58nxl1212178353.ps tmp/58nxl1212178353.png")
> system("convert tmp/6x6df1212178354.ps tmp/6x6df1212178354.png")
> system("convert tmp/73jmd1212178354.ps tmp/73jmd1212178354.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.041   1.397   2.228