R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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.

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(4143,4429,5219,4929,5761,5592,4163,4962,5208,4755,4491,5732,5731,5040,6102,4904,5369,5578,4619,4731,5011,5299,4146,4625,4736,4219,5116,4205,4121,5103,4300,4578,3809,5657,4248,3830,4736,4839,4411,4570,4104,4801,3953,3828,4440,4026,4109,4785,3224,3552,3940,3913,3681,4309,3830,4143,4087,3818,3380,3430,3458,3970,5260,5024,5634,6549,4676)
> 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] 67
> (np <- floor(n / par1))
[1] 5
> 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] 6 6 6 6 6 6 6 5 5 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,] 4143 5731 4736 4736 3224 3458
 [2,] 4429 5040 4219 4839 3552 3970
 [3,] 5219 6102 5116 4411 3940 5260
 [4,] 4929 4904 4205 4570 3913 5024
 [5,] 5761 5369 4121 4104 3681 5634
 [6,] 5592 5578 5103 4801 4309 6549
 [7,] 4163 4619 4300 3953 3830 4676
 [8,] 4962 4731 4578 3828 4143   NA
 [9,] 5208 5011 3809 4440 4087   NA
[10,] 4755 5299 5657 4026 3818   NA
[11,] 4491 4146 4248 4109 3380   NA
[12,] 5732 4625 3830 4785 3430   NA
> darr
       [,1]  [,2]  [,3]  [,4] [,5]  [,6]
 [1,]   286  -691  -517   103  328   512
 [2,]   790  1062   897  -428  388  1290
 [3,]  -290 -1198  -911   159  -27  -236
 [4,]   832   465   -84  -466 -232   610
 [5,]  -169   209   982   697  628   915
 [6,] -1429  -959  -803  -848 -479 -1873
 [7,]   799   112   278  -125  313    NA
 [8,]   246   280  -769   612  -56    NA
 [9,]  -453   288  1848  -414 -269    NA
[10,]  -264 -1153 -1409    83 -438    NA
[11,]  1241   479  -418   676   50    NA
[12,]    -1   111   906 -1561   28    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/rcomp/tmp/1ej5a1291992881.ps",horizontal=F,onefile=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/rcomp/tmp/2ej5a1291992881.ps",horizontal=F,onefile=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/rcomp/tmp/3pa4d1291992881.ps",horizontal=F,onefile=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/rcomp/tmp/4pa4d1291992881.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] [,10] [,11] [,12]
[1,] 3224.0 3552 3940.0 3913 3681 4309.0 3830.0 3828 3809  3818  4109  3430
[2,] 3458.0 3970 4411.0 4205 4104 4801.0 3953.0 4143 4087  4026  4109  3830
[3,] 4439.5 4324 5167.5 4737 4745 5340.5 4231.5 4578 4440  4755  4146  4625
[4,] 4736.0 4839 5260.0 4929 5634 5592.0 4619.0 4731 5011  5299  4248  4785
[5,] 5731.0 5040 6102.0 5024 5761 6549.0 4676.0 4962 5208  5657  4248  5732

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3615.149 3763.467 4619.868 4269.997 3758.101 4830.279 3801.908 4162.521
[2,] 5263.851 4884.533 5715.132 5204.003 5731.899 5850.721 4661.092 4993.479
         [,9]    [,10]    [,11]    [,12]
[1,] 3787.104 3855.501 4047.783 3950.199
[2,] 5092.896 5654.499 4244.217 5299.801

$out
[1] 4491 3380

$group
[1] 11 11

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

Warning message:
In bxp(list(stats = c(3224, 3458, 4439.5, 4736, 5731, 3552, 3970,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5pa4d1291992881.ps",horizontal=F,onefile=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,] -691.0 -428.0 -1198 -466.0 -169.0 -1873.0 -125  -56 -453 -1409  -418    -1
[2,] -517.0  388.0  -911 -232.0  209.0 -1429.0  112  -56 -414 -1153    50    -1
[3,]  194.5  843.5  -263  190.5  662.5  -903.5  278  246 -269  -438   479    28
[4,]  328.0 1062.0   -27  610.0  915.0  -803.0  313  280  288  -264   676   111
[5,]  512.0 1290.0   159  832.0  982.0  -479.0  313  612  288    83  1241   111

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]       [,6]     [,7]
[1,] -350.5523  408.7482 -833.2086 -352.6172  207.1072 -1307.2902 135.9739
[2,]  739.5523 1278.2518  307.2086  733.6172 1117.8928  -499.7098 420.0261
           [,8]      [,9]      [,10]     [,11]     [,12]
[1,]   8.583246 -765.0314 -1066.1652  36.66998 -51.13892
[2,] 483.416754  227.0314   190.1652 921.33002 107.13892

$out
[1]   799  -769  1848   906 -1561

$group
[1]  7  8  9 12 12

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

Warning message:
In bxp(list(stats = c(-691, -517, 194.5, 328, 512, -428, 388, 843.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/602ly1291992881.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]
[1,] 4143.0 4146.0 3809.0 3828.0 3224.0 3458
[2,] 4460.0 4678.0 4163.0 4065.0 3491.0 4323
[3,] 4945.5 5025.5 4274.0 4425.5 3824.0 5024
[4,] 5405.5 5473.5 4919.5 4760.5 4013.5 5447
[5,] 5761.0 6102.0 5657.0 4839.0 4309.0 6549

$n
[1] 12 12 12 12 12  7

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 4514.251 4662.667 3928.955 4108.278 3585.684 4352.765
[2,] 5376.749 5388.333 4619.045 4742.722 4062.316 5695.235

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(4143, 4460, 4945.5, 5405.5, 5761, 4146, 4678,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/702ly1291992881.ps",horizontal=F,onefile=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,] 4074.800 4146.00 4178.500
[2,] 4339.750 4381.75 4289.938
[3,] 4495.700 4601.50 4493.000
[4,] 4744.667 4750.00 4750.250
[5,] 5322.000 5167.50 5232.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4311.015 4433.539 4283.048
[2,] 4680.385 4769.461 4702.952

$out
[1] 5340.5

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(4074.8, 4339.75, 4495.7, 4744.66666666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ej5a1291992881.ps tmp/1ej5a1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ej5a1291992881.ps tmp/2ej5a1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/3pa4d1291992881.ps tmp/3pa4d1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/4pa4d1291992881.ps tmp/4pa4d1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pa4d1291992881.ps tmp/5pa4d1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/602ly1291992881.ps tmp/602ly1291992881.png",intern=TRUE))
character(0)
> try(system("convert tmp/702ly1291992881.ps tmp/702ly1291992881.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.650   0.820   2.443