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(244752,244576,241572,240541,236089,236997,264579,270349,269645,267037,258113,262813,267413,267366,264777,258863,254844,254868,277267,285351,286602,283042,276687,277915,277128,277103,275037,270150,267140,264993,287259,291186,292300,288186,281477,282656,280190,280408,276836,275216,274352,271311,289802,290726,292300,278506,269826,265861,269034,264176,255198,253353,246057,235372,258556,260993,254663,250643,243422,247105,248541,245039,237080,237085,225554,226839,247934,248333,246969,245098,246263,255765,264319,268347,273046,273963,267430,271993,292710)
> 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,] 244752 267413 277128 280190 269034 248541 264319
 [2,] 244576 267366 277103 280408 264176 245039 268347
 [3,] 241572 264777 275037 276836 255198 237080 273046
 [4,] 240541 258863 270150 275216 253353 237085 273963
 [5,] 236089 254844 267140 274352 246057 225554 267430
 [6,] 236997 254868 264993 271311 235372 226839 271993
 [7,] 264579 277267 287259 289802 258556 247934 292710
 [8,] 270349 285351 291186 290726 260993 248333     NA
 [9,] 269645 286602 292300 292300 254663 246969     NA
[10,] 267037 283042 288186 278506 250643 245098     NA
[11,] 258113 276687 281477 269826 243422 246263     NA
[12,] 262813 277915 282656 265861 247105 255765     NA
> darr
       [,1]  [,2]  [,3]   [,4]   [,5]   [,6]  [,7]
 [1,]  -176   -47   -25    218  -4858  -3502  4028
 [2,] -3004 -2589 -2066  -3572  -8978  -7959  4699
 [3,] -1031 -5914 -4887  -1620  -1845      5   917
 [4,] -4452 -4019 -3010   -864  -7296 -11531 -6533
 [5,]   908    24 -2147  -3041 -10685   1285  4563
 [6,] 27582 22399 22266  18491  23184  21095 20717
 [7,]  5770  8084  3927    924   2437    399    NA
 [8,]  -704  1251  1114   1574  -6330  -1364    NA
 [9,] -2608 -3560 -4114 -13794  -4020  -1871    NA
[10,] -8924 -6355 -6709  -8680  -7221   1165    NA
[11,]  4700  1228  1179  -3965   3683   9502    NA
[12,]  4600  -787 -2466   3173   1436   8554    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/1i98y1257968130.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/2cznr1257968130.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/3ulha1257968130.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/4yxf91257968130.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,] 244752 244576.0 237080.0 237085.0 225554 226839.0 247934.0 248333 246969.0
[2,] 256430 254607.5 248385.0 246947.0 241073 236184.5 261567.5 260993 254663.0
[3,] 267413 267366.0 264777.0 258863.0 254844 254868.0 277267.0 277850 278123.5
[4,] 273081 272725.0 274041.5 272056.5 267285 268152.0 288530.5 290726 292300.0
[5,] 280190 280408.0 276836.0 275216.0 274352 271993.0 292710.0 291186 292300.0
        [,10]    [,11]  [,12]
[1,] 245098.0 243422.0 247105
[2,] 250643.0 246263.0 255765
[3,] 272771.5 263969.5 264337
[4,] 283042.0 276687.0 277915
[5,] 288186.0 281477.0 282656

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

$conf
         [,1]     [,2]     [,3]   [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 257469.3 256546.5 249455.4 243868 239190.6 235777.5 261165.1 258671.3
[2,] 277356.7 278185.5 280098.6 273858 270497.4 273958.5 293368.9 297028.7
         [,9]    [,10]    [,11]    [,12]
[1,] 253846.4 251873.1 244345.0 250049.5
[2,] 302400.6 293669.9 283594.0 278624.5

$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(244752, 256430, 267413, 273081, 280190, 244576,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5trbl1257968130.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,] -3502.0 -8978.0 -5914 -11531.0 -3041.0 18491.0  399 -1364 -4114 -8924
[2,] -1839.0 -5765.5 -3366  -6914.5 -2594.0 20906.0  924 -1364 -4114 -8680
[3,]   -47.0 -3004.0 -1620  -4452.0    24.0 22266.0 3182   205 -3790 -6965
[4,]    96.5 -2327.5  -513  -3514.5  1096.5 22791.5 5770  1251 -2608 -6355
[5,]   218.0 -2066.0   917   -864.0  4563.0 23184.0 8084  1574 -1871 -6355
       [,11]   [,12]
[1,] -3965.0 -2466.0
[2,]  1179.0  -787.0
[3,]  2455.5  2304.5
[4,]  4700.0  4600.0
[5,]  9502.0  8554.0

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

$conf
          [,1]       [,2]        [,3]      [,4]      [,5]     [,6]       [,7]
[1,] -1202.849 -5057.1181 -3323.76557 -6482.425 -2179.907 21140.01   56.17348
[2,]  1108.849  -950.8819    83.76557 -2421.575  2227.907 23391.99 6307.82652
          [,8]      [,9]   [,10]     [,11]     [,12]
[1,] -1481.759 -4761.419 -8464.7  184.3413 -1170.289
[2,]  1891.759 -2818.581 -5465.3 4726.6587  5779.289

$out
[1]  -4858   4028   4699 -10685  27582  -6330 -13794   1165

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

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

Warning message:
In bxp(list(stats = c(-3502, -1839, -47, 96.5, 218, -8978, -5765.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/66vxw1257968130.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,] 236089.0 254844.0 264993.0 265861.0 235372.0 225554.0 264319.0
[2,] 241056.5 261820.0 272593.5 272831.5 246581.0 237082.5 267888.5
[3,] 251432.5 272050.0 279302.5 277671.0 254008.0 245680.5 271993.0
[4,] 265808.0 280478.5 287722.5 285105.0 259774.5 248133.5 273504.5
[5,] 270349.0 286602.0 292300.0 292300.0 269034.0 255765.0 273963.0

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

$conf
         [,1]     [,2]     [,3]   [,4]     [,5]     [,6]     [,7]
[1,] 240143.2 263539.7 272402.1 272073 247990.4 240640.1 268639.2
[2,] 262721.8 280560.3 286202.9 283269 260025.6 250720.9 275346.8

$out
[1] 292710

$group
[1] 7

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

Warning message:
In bxp(list(stats = c(236089, 241056.5, 251432.5, 265808, 270349,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/73mxd1257968130.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,] 251767.6 254844.0 252168.2
[2,] 259479.8 261416.2 260357.5
[3,] 264170.9 266071.5 264210.9
[4,] 271249.2 275019.2 271483.4
[5,] 274489.7 278123.5 276357.1

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 258802.7 259867.1 259136.3
[2,] 269539.0 272275.9 269285.5

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(251767.571428571, 259479.785714286, 264170.857142857,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> system("convert tmp/1i98y1257968130.ps tmp/1i98y1257968130.png")
> system("convert tmp/2cznr1257968130.ps tmp/2cznr1257968130.png")
> system("convert tmp/3ulha1257968130.ps tmp/3ulha1257968130.png")
> system("convert tmp/4yxf91257968130.ps tmp/4yxf91257968130.png")
> system("convert tmp/5trbl1257968130.ps tmp/5trbl1257968130.png")
> system("convert tmp/66vxw1257968130.ps tmp/66vxw1257968130.png")
> system("convert tmp/73mxd1257968130.ps tmp/73mxd1257968130.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.306   1.021  12.692