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(512238,519164,517009,509933,509127,500857,506971,569323,579714,577992,565464,547344,554788,562325,560854,555332,543599,536662,542722,593530,610763,612613,611324,594167,595454,590865,589379,584428,573100,567456,569028,620735,628884,628232,612117,595404,597141,593408,590072,579799,574205,572775,572942,619567,625809,619916,587625,565742,557274,560576,548854,531673,525919,511038,498662,555362,564591,541657,527070,509846,514258)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> ylab = ''
> xlab = ''
> main = ''
> #'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] 61
> (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 5 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]
 [1,] 512238 554788 595454 597141 557274 514258
 [2,] 519164 562325 590865 593408 560576     NA
 [3,] 517009 560854 589379 590072 548854     NA
 [4,] 509933 555332 584428 579799 531673     NA
 [5,] 509127 543599 573100 574205 525919     NA
 [6,] 500857 536662 567456 572775 511038     NA
 [7,] 506971 542722 569028 572942 498662     NA
 [8,] 569323 593530 620735 619567 555362     NA
 [9,] 579714 610763 628884 625809 564591     NA
[10,] 577992 612613 628232 619916 541657     NA
[11,] 565464 611324 612117 587625 527070     NA
[12,] 547344 594167 595404 565742 509846     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5] [,6]
 [1,]   6926   7537  -4589  -3733   3302   NA
 [2,]  -2155  -1471  -1486  -3336 -11722   NA
 [3,]  -7076  -5522  -4951 -10273 -17181   NA
 [4,]   -806 -11733 -11328  -5594  -5754   NA
 [5,]  -8270  -6937  -5644  -1430 -14881   NA
 [6,]   6114   6060   1572    167 -12376   NA
 [7,]  62352  50808  51707  46625  56700   NA
 [8,]  10391  17233   8149   6242   9229   NA
 [9,]  -1722   1850   -652  -5893 -22934   NA
[10,] -12528  -1289 -16115 -32291 -14587   NA
[11,] -18120 -17157 -16713 -21883 -17224   NA
[12,]   7444   1287   1737  -8468   4412   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/1nstr1226318618.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/279h61226318618.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/3mw031226318618.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/44wf01226318618.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,] 512238 519164 517009 509933 509127 500857 498662 555362 564591 541657
[2,] 514258 560576 548854 531673 525919 511038 506971 569323 579714 577992
[3,] 556031 562325 560854 555332 543599 536662 542722 593530 610763 612613
[4,] 595454 590865 589379 579799 573100 567456 569028 619567 625809 619916
[5,] 597141 593408 590072 584428 574205 572775 572942 620735 628884 628232
      [,11]  [,12]
[1,] 527070 509846
[2,] 565464 547344
[3,] 587625 565742
[4,] 611324 594167
[5,] 612117 595404

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

$conf
       [,1]     [,2]     [,3]     [,4]   [,5]     [,6]     [,7]     [,8]
[1,] 503657 540922.9 532219.1 521326.3 510261 496797.2 498872.7 558027.7
[2,] 608405 583727.1 589488.9 589337.7 576937 576526.8 586571.3 629032.3
         [,9]    [,10]    [,11]  [,12]
[1,] 578192.4 582989.6 555220.4 532657
[2,] 643333.6 642236.4 620029.6 598827

$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(512238, 514258, 556031, 595454, 597141, 519164,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5krml1226318618.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,] -4589 -3336 -17181 -11733 -8270  167 46625  6242 -5893 -16115 -18120  1287
[2,] -3733 -3336 -10273 -11328 -8270  167 50808  8149 -5893 -16115 -18120  1287
[3,]  3302 -2155  -7076  -5754 -6937 1572 51707  9229 -1722 -14587 -17224  1737
[4,]  6926 -1486  -5522  -5594 -5644 6060 56700 10391  -652 -12528 -17157  4412
[5,]  7537 -1471  -4951   -806 -5644 6114 62352 10391  1850 -12528 -16713  7444

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

$conf
          [,1]       [,2]       [,3]     [,4]      [,5]      [,6]     [,7]
[1,] -4229.623 -3462.2053 -10433.045 -9805.63 -8792.525 -2591.979 47543.73
[2,] 10833.623  -847.7947  -3718.955 -1702.37 -5081.475  5735.979 55870.27
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,]  7644.808 -5425.277 -17121.57 -17904.45 -471.1171
[2,] 10813.192  1981.277 -12052.43 -16543.55 3945.1171

$out
 [1] -11722  -1430 -14881 -12376  17233 -22934  -1289 -32291 -21883  -8468

$group
 [1]  2  5  5  6  8  9 10 10 11 12

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

Warning message:
In bxp(list(stats = c(-4589, -3733, 3302, 6926, 7537, -3336, -3336,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6xocb1226318618.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]
[1,] 500857.0 536662.0 567456.0 565742.0 498662.0 514258
[2,] 509530.0 549193.5 578764.0 573573.5 518478.5 514258
[3,] 518086.5 561589.5 593134.5 588848.5 536665.0 514258
[4,] 567393.5 602465.0 616426.0 608354.0 556318.0 514258
[5,] 579714.0 612613.0 628884.0 625809.0 564591.0 514258

$n
[1] 12 12 12 12 12  1

$conf
         [,1]   [,2]     [,3]     [,4]     [,5]   [,6]
[1,] 491694.6 537292 575956.6 572984.9 519406.2 514258
[2,] 544478.4 585887 610312.4 604712.1 553923.8 514258

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(500857, 509530, 518086.5, 567393.5, 579714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7wkrs1226318618.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,] 537757.6 536662.0 537999.5
[2,] 548711.5 549465.5 552329.6
[3,] 561867.1 561589.5 569936.0
[4,] 586211.7 590577.5 591419.5
[5,] 601952.2 612613.0 602761.5

$n
[1] 12 12 12

$conf
         [,1]   [,2]     [,3]
[1,] 544763.0 542838 552106.8
[2,] 578971.2 580341 587765.2

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(537757.6, 548711.5, 561867.1, 586211.7, 601952.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1nstr1226318618.ps tmp/1nstr1226318618.png")
> system("convert tmp/279h61226318618.ps tmp/279h61226318618.png")
> system("convert tmp/3mw031226318618.ps tmp/3mw031226318618.png")
> system("convert tmp/44wf01226318618.ps tmp/44wf01226318618.png")
> system("convert tmp/5krml1226318618.ps tmp/5krml1226318618.png")
> system("convert tmp/6xocb1226318618.ps tmp/6xocb1226318618.png")
> system("convert tmp/7wkrs1226318618.ps tmp/7wkrs1226318618.png")
> 
> 
> proc.time()
   user  system elapsed 
  3.376   2.369   3.689