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(274412,272433,268361,268586,264768,269974,304744,309365,308347,298427,289231,291975,294912,293488,290555,284736,281818,287854,316263,325412,326011,328282,317480,317539,313737,312276,309391,302950,300316,304035,333476,337698,335932,323931,313927,314485,313218,309664,302963,298989,298423,301631,329765,335083,327616,309119,295916,291413,291542,284678,276475,272566,264981,263290,296806,303598,286994,276427,266424,267153,268381,262522,255542,253158,243803,250741,280445,285257,270976,261076,255603,260376,263903,264291,263276,262572,256167,264221,293860)
> 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] 79
> (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] 7 7 7 7 7 7 7 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 274412 294912 313737 313218 291542 268381 263903
 [2,] 272433 293488 312276 309664 284678 262522 264291
 [3,] 268361 290555 309391 302963 276475 255542 263276
 [4,] 268586 284736 302950 298989 272566 253158 262572
 [5,] 264768 281818 300316 298423 264981 243803 256167
 [6,] 269974 287854 304035 301631 263290 250741 264221
 [7,] 304744 316263 333476 329765 296806 280445 293860
 [8,] 309365 325412 337698 335083 303598 285257     NA
 [9,] 308347 326011 335932 327616 286994 270976     NA
[10,] 298427 328282 323931 309119 276427 261076     NA
[11,] 289231 317480 313927 295916 266424 255603     NA
[12,] 291975 317539 314485 291413 267153 260376     NA
> darr
       [,1]   [,2]   [,3]   [,4]   [,5]   [,6]  [,7]
 [1,] -1979  -1424  -1461  -3554  -6864  -5859   388
 [2,] -4072  -2933  -2885  -6701  -8203  -6980 -1015
 [3,]   225  -5819  -6441  -3974  -3909  -2384  -704
 [4,] -3818  -2918  -2634   -566  -7585  -9355 -6405
 [5,]  5206   6036   3719   3208  -1691   6938  8054
 [6,] 34770  28409  29441  28134  33516  29704 29639
 [7,]  4621   9149   4222   5318   6792   4812    NA
 [8,] -1018    599  -1766  -7467 -16604 -14281    NA
 [9,] -9920   2271 -12001 -18497 -10567  -9900    NA
[10,] -9196 -10802 -10004 -13203 -10003  -5473    NA
[11,]  2744     59    558  -4503    729   4773    NA
[12,]  2937  -3802  -1267    129   1228   3527    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/1pl691257968391.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/2e6rd1257968391.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/3faba1257968391.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/48gh61257968391.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,] 263903.0 262522 255542.0 253158.0 243803.0 250741.0 280445 285257.0 270976
[2,] 271396.5 268362 265818.5 265579.0 260467.5 263755.5 295333 303598.0 286994
[3,] 291542.0 284678 276475.0 272566.0 264981.0 269974.0 304744 317388.5 317179
[4,] 304065.0 301576 296759.0 291862.5 290120.5 294742.5 323014 335083.0 327616
[5,] 313737.0 312276 309391.0 302950.0 300316.0 304035.0 333476 337698.0 335932
      [,10]    [,11]  [,12]
[1,] 261076 255603.0 260376
[2,] 276427 266424.0 267153
[3,] 303773 292573.5 291694
[4,] 323931 313927.0 314485
[5,] 328282 317480.0 317539

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 272032.9 264843.1 257997.8 256869.9 247272.7 251469.1 288213.4 297079.7
[2,] 311051.1 304512.9 294952.2 288262.1 282689.3 288478.9 321274.6 337697.3
         [,9]    [,10]    [,11]    [,12]
[1,] 290976.5 273131.4 261932.5 261163.3
[2,] 343381.5 334414.6 323214.5 322224.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(263903, 271396.5, 291542, 304065, 313737,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5yf371257968391.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,] -6864.0 -8203.0 -6441.0 -9355 3208.0 28134 4222 -16604.0 -12001.0 -13203.0
[2,] -4706.5 -6840.5 -4896.5 -6995 3463.5 28925 4621 -14281.0 -12001.0 -10802.0
[3,] -1979.0 -4072.0 -3909.0 -3818 5206.0 29639 5065  -4616.5 -10243.5 -10003.5
[4,] -1442.5 -2909.0 -1544.0 -2776 6487.0 31610 6792  -1018.0  -9900.0  -9196.0
[5,]   388.0 -1015.0   225.0  -566 8054.0 34770 9149    599.0  -9900.0  -9196.0
      [,11]   [,12]
[1,]   59.0 -3802.0
[2,]   59.0 -1267.0
[3,]  643.5   678.5
[4,] 2744.0  2937.0
[5,] 4773.0  3527.0

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

$conf
            [,1]      [,2]      [,3]      [,4]     [,5]     [,6]     [,7]
[1,] -3928.20814 -6419.828 -5911.059 -6337.519 3400.415 28035.56 3664.635
[2,]   -29.79186 -1724.172 -1906.941 -1298.481 7011.585 31242.44 6465.365
           [,8]       [,9]      [,10]     [,11]     [,12]
[1,] -13171.563 -11598.713 -11039.422 -1088.412 -2033.216
[2,]   3938.563  -8888.287  -8967.578  2375.412  3390.216

$out
[1]  -1691   2271 -18497  -5473  -4503

$group
[1]  5  9  9 10 11

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

Warning message:
In bxp(list(stats = c(-6864, -4706.5, -1979, -1442.5, 388, -8203,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/63sf91257968391.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,] 264768.0 281818.0 300316.0 291413 263290.0 243803.0 262572
[2,] 269280.0 289204.5 306713.0 298706 266788.5 254350.0 262924
[3,] 281821.5 305587.5 313832.0 306041 276451.0 260726.0 263903
[4,] 301585.5 321475.5 328703.5 320417 289268.0 269678.5 264256
[5,] 309365.0 328282.0 337698.0 335083 303598.0 285257.0 264291

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

$conf
         [,1]     [,2]   [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 267086.7 290868.5 303802 296138.5 266197.9 253734.6 263107.6
[2,] 296556.3 320306.5 323862 315943.5 286704.1 267717.4 264698.4

$out
[1] 256167 293860

$group
[1] 7 7

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

Warning message:
In bxp(list(stats = c(264768, 269280, 281821.5, 301585.5, 309365,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/725xv1257968391.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,] 272896.6 264981.0 275294.0
[2,] 279294.3 274520.5 280268.9
[3,] 289175.0 291618.0 289252.9
[4,] 303726.0 304258.5 305125.5
[5,] 316068.8 317388.5 318852.5

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 278031.5 278054.3 277915.6
[2,] 300318.5 305181.7 300590.1

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(272896.571428571, 279294.285714286, 289174.964285714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/1pl691257968391.ps tmp/1pl691257968391.png")
> system("convert tmp/2e6rd1257968391.ps tmp/2e6rd1257968391.png")
> system("convert tmp/3faba1257968391.ps tmp/3faba1257968391.png")
> system("convert tmp/48gh61257968391.ps tmp/48gh61257968391.png")
> system("convert tmp/5yf371257968391.ps tmp/5yf371257968391.png")
> system("convert tmp/63sf91257968391.ps tmp/63sf91257968391.png")
> system("convert tmp/725xv1257968391.ps tmp/725xv1257968391.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.206   0.949   1.570