R version 2.6.2 (2008-02-08)
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(412218,492982,501175,506949,491951,565441,575772,572507,548752,621486,615237,588119,552118,624447,615482,595130,526883,562254,545427,522084,483414,528797,532749,511380,472941,516118,502940,476118,432418,475525,453638,431417,390934,436414,418451,399528,367749,423433,420450,415906,392949,453203,455926,451879,434996,498811,505940,517395,508456,585132,587971,584027,557196,613433,600049,588993,559271,622580,616645,603243,557949,608882,582930,570492,542907,598067,568717,551773,514465,569055,528897,515229,481141,535612,498547,478587,445911,503412,469797,458365,436761,502205,481627,473698,457200,521671,513354,515369,505652,575676,555865,559504,540994,605635,600315,588224,569861,625950,601554,587760,573307,621764,570214,547034,511873,553870)
> 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] 106
> (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 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,] 412218 552118 472941 367749 508456 557949 481141 457200 569861
 [2,] 492982 624447 516118 423433 585132 608882 535612 521671 625950
 [3,] 501175 615482 502940 420450 587971 582930 498547 513354 601554
 [4,] 506949 595130 476118 415906 584027 570492 478587 515369 587760
 [5,] 491951 526883 432418 392949 557196 542907 445911 505652 573307
 [6,] 565441 562254 475525 453203 613433 598067 503412 575676 621764
 [7,] 575772 545427 453638 455926 600049 568717 469797 555865 570214
 [8,] 572507 522084 431417 451879 588993 551773 458365 559504 547034
 [9,] 548752 483414 390934 434996 559271 514465 436761 540994 511873
[10,] 621486 528797 436414 498811 622580 569055 502205 605635 553870
[11,] 615237 532749 418451 505940 616645 528897 481627 600315     NA
[12,] 588119 511380 399528 517395 603243 515229 473698 588224     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,]  80764  72329  43177  55684  76676  50933  54471  64471  56089
 [2,]   8193  -8965 -13178  -2983   2839 -25952 -37065  -8317 -24396
 [3,]   5774 -20352 -26822  -4544  -3944 -12438 -19960   2015 -13794
 [4,] -14998 -68247 -43700 -22957 -26831 -27585 -32676  -9717 -14453
 [5,]  73490  35371  43107  60254  56237  55160  57501  70024  48457
 [6,]  10331 -16827 -21887   2723 -13384 -29350 -33615 -19811 -51550
 [7,]  -3265 -23343 -22221  -4047 -11056 -16944 -11432   3639 -23180
 [8,] -23755 -38670 -40483 -16883 -29722 -37308 -21604 -18510 -35161
 [9,]  72734  45383  45480  63815  63309  54590  65444  64641  41997
[10,]  -6249   3952 -17963   7129  -5935 -40158 -20578  -5320     NA
[11,] -27118 -21369 -18923  11455 -13402 -13668  -7929 -12091     NA
[12,] -36001 -38439 -31779  -8939 -45294 -34088 -16498 -18363     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/144ad1209378902.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/2y12p1209378902.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/39p3u1209378902.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/45jl71209378902.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,] 367749 423433 420450 415906 392949 453203 453638 431417 390934 436414
[2,] 457200 516118 501175 478587 445911 503412 469797 458365 436761 502205
[3,] 481141 535612 513354 515369 505652 565441 555865 547034 511873 553870
[4,] 552118 608882 587971 584027 542907 598067 570214 559504 540994 605635
[5,] 569861 625950 615482 595130 573307 621764 600049 588993 559271 622580
        [,11]    [,12]
[1,] 418451.0 399528.0
[2,] 493783.5 492539.0
[3,] 530823.0 516312.0
[4,] 607776.0 588171.5
[5,] 616645.0 603243.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 431150.9 486756.3 467641.4 459837.3 454567.4 515589.4 502978.7 493767.5
[2,] 531131.1 584467.7 559066.6 570900.7 556736.6 615292.6 608751.3 600300.5
         [,9]    [,10]    [,11]    [,12]
[1,] 456977.0 499396.9 467145.2 462890.3
[2,] 566769.0 608343.1 594500.8 569733.7

$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(367749, 457200, 481141, 552118, 569861, 423433,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5awkn1209378903.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]
[1,] 43177 -37065 -26822 -43700 35371 -51550 -23343 -40483 41997 -40158.0
[2,] 54471 -24396 -19960 -32676 48457 -29350 -22221 -37308 45480 -19270.5
[3,] 56089  -8965 -12438 -26831 56237 -19811 -11432 -29722 63309  -6092.0
[4,] 72329  -2983  -3944 -14998 60254 -13384  -4047 -21604 64641   -684.0
[5,] 80764   8193   5774  -9717 73490  10331   3639 -16883 72734   7129.0
      [,11]    [,12]
[1,] -27118 -45294.0
[2,] -20146 -37220.0
[3,] -13535 -32933.5
[4,] -10010 -17430.5
[5,]  -7929  -8939.0

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

$conf
         [,1]       [,2]       [,3]      [,4]     [,5]      [,6]      [,7]
[1,] 46683.79 -20242.513 -20873.093 -36141.41 50023.91 -28219.76 -21003.64
[2,] 65494.21   2312.513  -4002.907 -17520.59 62450.09 -11402.24  -1860.36
          [,8]     [,9]      [,10]      [,11]    [,12]
[1,] -37992.77 53217.54 -16474.686 -19197.115 -43988.2
[2,] -21451.23 73400.46   4290.686  -7872.885 -21878.8

$out
[1] -68247  11455

$group
[1]  4 11

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

Warning message:
In bxp(list(stats = c(43177, 54471, 56089, 72329, 80764, -37065,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/68e0f1209378903.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,] 412218.0 483414.0 390934.0 367749.0 557196 514465 436761 457200.0 511873.0
[2,] 497078.5 524483.5 424934.0 418178.0 571649 535902 464081 514361.5 553870.0
[3,] 557096.5 539088.0 445026.0 443437.5 588482 563333 479864 548429.5 571760.5
[4,] 581945.5 578692.0 475821.5 477368.5 608338 576711 500376 581950.0 601554.0
[5,] 621486.0 624447.0 516118.0 517395.0 622580 608882 535612 605635.0 625950.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 518388.1 514363.1 421815.9 416440.3 571747.9 544719.7 463309.6 517601.9
[2,] 595804.9 563812.9 468236.1 470434.7 605216.1 581946.3 496418.4 579257.1
         [,9]
[1,] 547935.7
[2,] 595585.3

$out
[1] 508456

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(412218, 497078.5, 557096.5, 581945.5, 621486,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/78pw91209378903.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,] 486625.9 481141.0 488877.5
[2,] 508485.0 512613.5 506796.8
[3,] 529207.9 523567.5 537940.0
[4,] 542865.0 550452.0 551346.6
[5,] 552086.1 565441.0 562500.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 513527.0 506309.1 517620.5
[2,] 544888.9 540825.9 558259.5

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(486625.888888889, 508485, 529207.944444444,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/144ad1209378902.ps tmp/144ad1209378902.png")
> system("convert tmp/2y12p1209378902.ps tmp/2y12p1209378902.png")
> system("convert tmp/39p3u1209378902.ps tmp/39p3u1209378902.png")
> system("convert tmp/45jl71209378902.ps tmp/45jl71209378902.png")
> system("convert tmp/5awkn1209378903.ps tmp/5awkn1209378903.png")
> system("convert tmp/68e0f1209378903.ps tmp/68e0f1209378903.png")
> system("convert tmp/78pw91209378903.ps tmp/78pw91209378903.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.477   0.963   1.834