R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(357704,281463,282445,319107,315278,328499,321151,328025,326280,313444,319639,324067,386918,293009,294822,338844,335407,345080,350608,351285,355147,332791,335615,343202,404868,317902,313552,361505,351436,373350,366310,361669,375078,345547,348117,356089,416856,328087,322747,373626,358275,391287,376371,371848,387261,353159,367855,376822,425283,342191,344062,373587,370144,399979,380431,385909,384798,352554,352479,338788,387964,313593,304056,334149,336155,354668,351418,354316,359483,330411,344726,347175)
> 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] 72
> (np <- floor(n / par1))
[1] 6
> 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 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6] [,7]
 [1,] 357704 386918 404868 416856 425283 387964   NA
 [2,] 281463 293009 317902 328087 342191 313593   NA
 [3,] 282445 294822 313552 322747 344062 304056   NA
 [4,] 319107 338844 361505 373626 373587 334149   NA
 [5,] 315278 335407 351436 358275 370144 336155   NA
 [6,] 328499 345080 373350 391287 399979 354668   NA
 [7,] 321151 350608 366310 376371 380431 351418   NA
 [8,] 328025 351285 361669 371848 385909 354316   NA
 [9,] 326280 355147 375078 387261 384798 359483   NA
[10,] 313444 332791 345547 353159 352554 330411   NA
[11,] 319639 335615 348117 367855 352479 344726   NA
[12,] 324067 343202 356089 376822 338788 347175   NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6] [,7]
 [1,] -76241 -93909 -86966 -88769 -83092 -74371   NA
 [2,]    982   1813  -4350  -5340   1871  -9537   NA
 [3,]  36662  44022  47953  50879  29525  30093   NA
 [4,]  -3829  -3437 -10069 -15351  -3443   2006   NA
 [5,]  13221   9673  21914  33012  29835  18513   NA
 [6,]  -7348   5528  -7040 -14916 -19548  -3250   NA
 [7,]   6874    677  -4641  -4523   5478   2898   NA
 [8,]  -1745   3862  13409  15413  -1111   5167   NA
 [9,] -12836 -22356 -29531 -34102 -32244 -29072   NA
[10,]   6195   2824   2570  14696    -75  14315   NA
[11,]   4428   7587   7972   8967 -13691   2449   NA
[12,]  62851  61666  60767  48461  49176     NA   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/1ej171272314970.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/2ej171272314970.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/3psia1272314970.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/4psia1272314970.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]
[1,] 357704 281463.0 282445 319107.0 315278.0 328499 321151 328025.0 326280.0
[2,] 386918 293009.0 294822 334149.0 335407.0 345080 350608 351285.0 355147.0
[3,] 396416 315747.5 308804 350174.5 343795.5 364009 358864 357992.5 367280.5
[4,] 416856 328087.0 322747 373587.0 358275.0 391287 376371 371848.0 384798.0
[5,] 425283 342191.0 344062 373626.0 370144.0 399979 380431 385909.0 387261.0
      [,10]    [,11]    [,12]
[1,] 313444 319639.0 324067.0
[2,] 330411 335615.0 338788.0
[3,] 339169 346421.5 345188.5
[4,] 352554 352479.0 356089.0
[5,] 353159 367855.0 376822.0

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

$conf
       [,1]     [,2]     [,3]     [,4]     [,5]   [,6]   [,7]     [,8]     [,9]
[1,] 377105 293121.1 290791.5 324735.7 329044.9 334204 342246 344728.7 348154.6
[2,] 415727 338373.9 326816.5 375613.3 358546.1 393814 375482 371256.3 386406.4
        [,10]    [,11]    [,12]
[1,] 324886.0 335543.7 334028.8
[2,] 353452.0 357299.3 356348.2

$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(357704, 386918, 396416, 416856, 425283, 281463,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5hk0v1272314970.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,] -93909 -9537 29525 -15351  9673.0 -19548 -4641.0 -1745.0 -34102.0   -75.0
[2,] -88769 -5340 30093 -10069 13221.0 -14916 -4523.0 -1111.0 -32244.0  2570.0
[3,] -85029 -1684 40342  -3636 20213.5  -7194  1787.5  4514.5 -29301.5  4509.5
[4,] -76241  1813 47953  -3437 29835.0  -3250  5478.0 13409.0 -22356.0 14315.0
[5,] -74371  1871 50879   2006 33012.0   5528  6874.0 15413.0 -12836.0 14696.0
      [,11] [,12]
[1,] 2449.0 48461
[2,] 2449.0 49176
[3,] 6007.5 60767
[4,] 7972.0 61666
[5,] 8967.0 62851

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

$conf
          [,1]      [,2]     [,3]       [,4]      [,5]        [,6]      [,7]
[1,] -93109.96 -6297.916 28821.72 -7913.8542  9496.933 -14718.9468 -4663.468
[2,] -76948.04  2929.916 51862.28   641.8542 30930.067    330.9468  8238.468
          [,8]      [,9]     [,10]    [,11]   [,12]
[1,] -4851.369 -35679.58 -3066.404 2444.987 51941.6
[2,] 13880.369 -22923.42 12085.404 9570.013 69592.4

$out
[1] -13691

$group
[1] 11

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

Warning message:
In bxp(list(stats = c(-93909, -88769, -85029, -76241, -74371, -9537,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6hk0v1272314970.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]
[1,] 313444.0 332791.0 313552 322747.0 338788.0 304056.0   NA
[2,] 314361.0 334099.0 346832 355717.0 348270.5 332280.0   NA
[3,] 320395.0 341023.0 358797 372737.0 371865.5 345950.5   NA
[4,] 327152.5 350946.5 369830 382041.5 385353.5 354492.0   NA
[5,] 328499.0 355147.0 375078 416856.0 425283.0 359483.0   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 314560.7 333338.7 348307.5 360730.2 354951.7 335819.5   NA
[2,] 326229.3 348707.3 369286.5 384743.8 388779.3 356081.5   NA

$out
[1] 357704 281463 282445 386918 293009 294822 404868 387964

$group
[1] 1 1 1 2 2 2 3 6

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

Warning message:
In bxp(list(stats = c(313444, 314361, 320395, 327152.5, 328499,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7hk0v1272314970.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,] 312707.5 315747.5 311847.9
[2,] 341216.8 341482.2 342772.4
[3,] 348913.4 348298.0 349910.3
[4,] 361758.2 361436.5 364736.5
[5,] 365477.2 367280.5 369299.5

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 339544.3 339196.7 339892.3
[2,] 358282.5 357399.3 359928.3

$out
[1] 396598.8 310280.7 396416.0 308804.0 400519.2 308789.4

$group
[1] 1 1 2 2 3 3

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

Warning message:
In bxp(list(stats = c(312707.5, 341216.75, 348913.416666667, 361758.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ej171272314970.ps tmp/1ej171272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ej171272314970.ps tmp/2ej171272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/3psia1272314970.ps tmp/3psia1272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/4psia1272314970.ps tmp/4psia1272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/5hk0v1272314970.ps tmp/5hk0v1272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/6hk0v1272314970.ps tmp/6hk0v1272314970.png",intern=TRUE))
character(0)
> try(system("convert tmp/7hk0v1272314970.ps tmp/7hk0v1272314970.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.216   0.931   3.192