R version 2.12.1 (2010-12-16)
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(335378,331533,318620,320838,326738,339196,336544,316150,309160,314901,313468,309914,306532,299149,288776,289441,293886,305957,305486,285179,280244,290616,292133,300414,306122,304198,300074,307836,313761,329573,327225,306307,316532,323953,327103,337606,344178,342892,344693,355072,372423,378519,379421,366121,367866,367980,370804,376126,377795,376740,372402,380575,387267,387820,389827,376911,377367,378670,385705,387495,386586,388884,384359,392058,384502,381903,382842,371551,369167,370861,378205,377632,376974,370490,364431,372720,375787,375021,374556,355178,352487,355996,353467,356427,354846)
> 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] 85
> (np <- floor(n / par1))
[1] 7
> 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] 8 7 7 7 7 7 7 7 7 7 7 7
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]
 [1,] 335378 306532 306122 344178 377795 386586 376974 354846
 [2,] 331533 299149 304198 342892 376740 388884 370490     NA
 [3,] 318620 288776 300074 344693 372402 384359 364431     NA
 [4,] 320838 289441 307836 355072 380575 392058 372720     NA
 [5,] 326738 293886 313761 372423 387267 384502 375787     NA
 [6,] 339196 305957 329573 378519 387820 381903 375021     NA
 [7,] 336544 305486 327225 379421 389827 382842 374556     NA
 [8,] 316150 285179 306307 366121 376911 371551 355178     NA
 [9,] 309160 280244 316532 367866 377367 369167 352487     NA
[10,] 314901 290616 323953 367980 378670 370861 355996     NA
[11,] 313468 292133 327103 370804 385705 378205 353467     NA
[12,] 309914 300414 337606 376126 387495 377632 356427     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7] [,8]
 [1,]  -3845  -7383  -1924  -1286  -1055   2298  -6484   NA
 [2,] -12913 -10373  -4124   1801  -4338  -4525  -6059   NA
 [3,]   2218    665   7762  10379   8173   7699   8289   NA
 [4,]   5900   4445   5925  17351   6692  -7556   3067   NA
 [5,]  12458  12071  15812   6096    553  -2599   -766   NA
 [6,]  -2652   -471  -2348    902   2007    939   -465   NA
 [7,] -20394 -20307 -20918 -13300 -12916 -11291 -19378   NA
 [8,]  -6990  -4935  10225   1745    456  -2384  -2691   NA
 [9,]   5741  10372   7421    114   1303   1694   3509   NA
[10,]  -1433   1517   3150   2824   7035   7344  -2529   NA
[11,]  -3554   8281  10503   5322   1790   -573   2960   NA
[12,]  -3382   5708   6572   1669   -909   -658  -1581   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/13vz91301929158.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/2ww7p1301929158.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/3znha1301929158.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/4a40q1301929158.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]
[1,] 306122.0 299149.0 288776.0 289441.0 293886.0 305957.0 305486.0 285179.0
[2,] 320955.0 317865.5 309347.0 314337.0 320249.5 334384.5 331884.5 311228.5
[3,] 349512.0 342892.0 344693.0 355072.0 372423.0 375021.0 374556.0 355178.0
[4,] 377384.5 373615.0 368416.5 376647.5 380144.5 380211.0 381131.5 368836.0
[5,] 386586.0 388884.0 384359.0 392058.0 387267.0 387820.0 389827.0 376911.0
         [,9]    [,10]    [,11]  [,12]
[1,] 280244.0 290616.0 292133.0 300414
[2,] 312846.0 319427.0 320285.5 323760
[3,] 352487.0 355996.0 353467.0 356427
[4,] 368516.5 369420.5 374504.5 376879
[5,] 377367.0 378670.0 385705.0 387495

$n
 [1] 8 7 7 7 7 7 7 7 7 7 7 7

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 317989.7 309599.3 309417.6 317861.2 336654.7 347654.2 345146.5 320775.7
[2,] 381034.3 376184.7 379968.4 392282.8 408191.3 402387.8 403965.5 389580.3
         [,9]    [,10]    [,11]    [,12]
[1,] 319241.5 326140.7 321088.3 324705.2
[2,] 385732.5 385851.3 385845.7 388148.8

$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(306122, 320955, 349512, 377384.5, 386586,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5ymrv1301929158.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]
[1,] -7383.0 -12913   665.0 3067.0 -2599.0 -2652.0 -20918.0 -6990.0   114.0
[2,] -5164.5  -8216  4958.5 3756.0  -106.5 -1409.5 -20350.5 -3813.0  1498.5
[3,] -1924.0  -4525  7762.0 5900.0  6096.0  -465.0 -19378.0 -2384.0  3509.0
[4,] -1170.5  -4231  8231.0 6308.5 12264.5   920.5 -13108.0  1100.5  6581.0
[5,]  2298.0  -4124 10379.0 6692.0 15812.0  2007.0 -11291.0  1745.0 10372.0
       [,10]   [,11]   [,12]
[1,] -2529.0 -3554.0 -3382.0
[2,]    42.0   608.5 -1245.0
[3,]  2824.0  2960.0  -658.0
[4,]  5092.5  6801.5  3688.5
[5,]  7344.0 10503.0  6572.0

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

$conf
           [,1]      [,2]     [,3]     [,4]      [,5]       [,6]      [,7]
[1,] -4309.1524 -6904.778 5807.716 4375.688 -1291.762 -1856.4384 -23703.10
[2,]   461.1524 -2145.222 9716.284 7424.312 13483.762   926.4384 -15052.90
           [,8]     [,9]     [,10]     [,11]     [,12]
[1,] -5318.2629  473.813 -192.0771 -738.3597 -3604.207
[2,]   550.2629 6544.187 5840.0771 6658.3597  2288.207

$out
[1]  1801 17351 -7556 10225

$group
[1] 2 4 4 8

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

Warning message:
In bxp(list(stats = c(-7383, -5164.5, -1924, -1170.5, 2298, -12913,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/67uzr1301929158.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]     [,7]   [,8]
[1,] 309160.0 280244.0 300074.0 342892.0 372402.0 369167.0 352487.0 354846
[2,] 314184.5 289108.5 306214.5 349882.5 377139.0 374591.5 355587.0 354846
[3,] 319729.0 293009.5 315146.5 367923.0 379622.5 382372.5 367460.5 354846
[4,] 333455.5 302950.0 327164.0 374274.5 387381.0 385544.0 374788.5 354846
[5,] 339196.0 306532.0 337606.0 379421.0 389827.0 392058.0 376974.0 354846

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]   [,6]     [,7]   [,8]
[1,] 310939.4 286696.3 305591.3 356797.6 374951.1 377377 358702.6 354846
[2,] 328518.6 299322.7 324701.7 379048.4 384293.9 387368 376218.4 354846

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(309160, 314184.5, 319729, 333455.5, 339196,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/743b51301929158.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,] 338974.7 342892.0 338881.8
[2,] 341455.3 350999.5 342552.5
[3,] 345673.2 355125.0 346567.6
[4,] 349998.4 364425.0 351496.2
[5,] 356855.6 375021.0 357297.8

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 341776.6 349001.5 342488.3
[2,] 349569.8 361248.5 350646.9

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(338974.714285714, 341455.285714286, 345673.214285714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/13vz91301929158.ps tmp/13vz91301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ww7p1301929158.ps tmp/2ww7p1301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/3znha1301929158.ps tmp/3znha1301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/4a40q1301929158.ps tmp/4a40q1301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ymrv1301929158.ps tmp/5ymrv1301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/67uzr1301929158.ps tmp/67uzr1301929158.png",intern=TRUE))
character(0)
> try(system("convert tmp/743b51301929158.ps tmp/743b51301929158.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.148   0.484   2.612