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(178600
+ ,181600
+ ,178500
+ ,176700
+ ,183500
+ ,175900
+ ,179800
+ ,186500
+ ,182700
+ ,182800
+ ,178600
+ ,183300
+ ,182900
+ ,191400
+ ,189300
+ ,192200
+ ,187900
+ ,193900
+ ,189100
+ ,193100
+ ,194800
+ ,200200
+ ,211500
+ ,202100
+ ,200300
+ ,199200
+ ,204900
+ ,207300
+ ,200000
+ ,197700
+ ,202200
+ ,200200
+ ,208300
+ ,215100
+ ,210700
+ ,208100
+ ,209000
+ ,211000
+ ,210200
+ ,205500
+ ,211400
+ ,211700
+ ,209300
+ ,207500
+ ,203300
+ ,207100
+ ,206900
+ ,228700
+ ,226900
+ ,265000
+ ,227100
+ ,228100
+ ,226500
+ ,225200
+ ,217800
+ ,221300
+ ,215300
+ ,231300
+ ,227100
+ ,237800
+ ,230200
+ ,233400
+ ,231100
+ ,237200
+ ,243700
+ ,239700
+ ,248400
+ ,241000
+ ,254500
+ ,242800
+ ,268300
+ ,253900
+ ,262100
+ ,264100
+ ,261000
+ ,269300
+ ,260400
+ ,263200
+ ,279200
+ ,272200
+ ,269200
+ ,289600
+ ,283200
+ ,284300
+ ,283000
+ ,289100
+ ,289600
+ ,289100
+ ,287400
+ ,279600
+ ,289300
+ ,295000
+ ,299600
+ ,293600
+ ,294400
+ ,290200
+ ,301000
+ ,307900
+ ,298800
+ ,310300
+ ,293900
+ ,305000
+ ,311300
+ ,317300
+ ,296200
+ ,306800
+ ,291800
+ ,301900
+ ,314600
+ ,321500
+ ,329400
+ ,311700
+ ,309700
+ ,306500
+ ,307100
+ ,301300
+ ,292200
+ ,310100
+ ,316800
+ ,284400
+ ,284600
+ ,301200
+ ,287600
+ ,314300
+ ,298200
+ ,299400
+ ,301900
+ ,265500
+ ,287100
+ ,274000
+ ,290100
+ ,263100
+ ,245200
+ ,258600
+ ,259800
+ ,269800
+ ,274600
+ ,274800
+ ,271100
+ ,257800
+ ,290300
+ ,262200
+ ,270000
+ ,290600)
> 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] 144
> (np <- floor(n / par1))
[1] 12
> 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] 12 12 12 12 12 12 12 12 12 12 12 12
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 178600 182900 200300 209000 226900 230200 262100 283000 301000 314600
 [2,] 181600 191400 199200 211000 265000 233400 264100 289100 307900 321500
 [3,] 178500 189300 204900 210200 227100 231100 261000 289600 298800 329400
 [4,] 176700 192200 207300 205500 228100 237200 269300 289100 310300 311700
 [5,] 183500 187900 200000 211400 226500 243700 260400 287400 293900 309700
 [6,] 175900 193900 197700 211700 225200 239700 263200 279600 305000 306500
 [7,] 179800 189100 202200 209300 217800 248400 279200 289300 311300 307100
 [8,] 186500 193100 200200 207500 221300 241000 272200 295000 317300 301300
 [9,] 182700 194800 208300 203300 215300 254500 269200 299600 296200 292200
[10,] 182800 200200 215100 207100 231300 242800 289600 293600 306800 310100
[11,] 178600 211500 210700 206900 227100 268300 283200 294400 291800 316800
[12,] 183300 202100 208100 228700 237800 253900 284300 290200 301900 284400
       [,11]  [,12] [,13]
 [1,] 284600 245200    NA
 [2,] 301200 258600    NA
 [3,] 287600 259800    NA
 [4,] 314300 269800    NA
 [5,] 298200 274600    NA
 [6,] 299400 274800    NA
 [7,] 301900 271100    NA
 [8,] 265500 257800    NA
 [9,] 287100 290300    NA
[10,] 274000 262200    NA
[11,] 290100 270000    NA
[12,] 263100 290600    NA
> darr
       [,1]  [,2]  [,3]  [,4]   [,5]   [,6]  [,7]  [,8]   [,9]  [,10]  [,11]
 [1,]  3000  8500 -1100  2000  38100   3200  2000  6100   6900   6900  16600
 [2,] -3100 -2100  5700  -800 -37900  -2300 -3100   500  -9100   7900 -13600
 [3,] -1800  2900  2400 -4700   1000   6100  8300  -500  11500 -17700  26700
 [4,]  6800 -4300 -7300  5900  -1600   6500 -8900 -1700 -16400  -2000 -16100
 [5,] -7600  6000 -2300   300  -1300  -4000  2800 -7800  11100  -3200   1200
 [6,]  3900 -4800  4500 -2400  -7400   8700 16000  9700   6300    600   2500
 [7,]  6700  4000 -2000 -1800   3500  -7400 -7000  5700   6000  -5800 -36400
 [8,] -3800  1700  8100 -4200  -6000  13500 -3000  4600 -21100  -9100  21600
 [9,]   100  5400  6800  3800  16000 -11700 20400 -6000  10600  17900 -13100
[10,] -4200 11300 -4400  -200  -4200  25500 -6400   800 -15000   6700  16100
[11,]  4700 -9400 -2600 21800  10700 -14400  1100 -4200  10100 -32400 -27000
[12,]  -400 -1800   900 -1800  -7600   8200 -1300 10800  12700    200 -17900
       [,12] [,13]
 [1,]  13400    NA
 [2,]   1200    NA
 [3,]  10000    NA
 [4,]   4800    NA
 [5,]    200    NA
 [6,]  -3700    NA
 [7,] -13300    NA
 [8,]  32500    NA
 [9,] -28100    NA
[10,]   7800    NA
[11,]  20600    NA
[12,]     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/1qped1272035994.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/2qped1272035994.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/3qped1272035994.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/40yvy1272035994.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,] 178600 181600 178500 176700 183500 175900 179800 186500 182700 182800
[2,] 204650 205100 207550 206400 205700 204700 205750 203850 205800 211100
[3,] 237700 261350 245450 253250 252050 251450 259750 249400 261850 252500
[4,] 283800 295150 288600 299700 290650 289500 295600 283600 291250 291600
[5,] 314600 321500 329400 314300 309700 306500 311300 317300 299600 310100
      [,11]  [,12]
[1,] 178600 183300
[2,] 211100 218400
[3,] 269150 258500
[4,] 290950 287300
[5,] 316800 301900

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 201599.2 220277.6 208482.6 210695.2 213303.7 212772.2 218768.8 213025.5
[2,] 273800.8 302422.4 282417.4 295804.8 290796.3 290127.8 300731.2 285774.5
         [,9]    [,10]    [,11]    [,12]
[1,] 222875.7 215783.4 232729.9 227074.2
[2,] 300824.3 289216.6 305570.1 289925.8

$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(178600, 204650, 237700, 283800, 314600, 181600,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/50yvy1272035994.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]
[1,] -1100 -13600 -4700 -16400 -7800 -7400 -13300 -21100 -28100 -15000 -32400
[2,]  2500  -6100 -1150  -8100 -3600 -3050  -7200  -5100  -8850  -4300 -11900
[3,]  6500  -2200  2650  -1850  -550  3200  -1900   -650   4600    300   -750
[4,] 10950    850  9150   5350  2000  7500   4850  10800  13300   9550  10400
[5,] 16600   7900 11500   6800  6000 16000   6700  32500  20400  25500  21800
     [,12]
[1,] -7600
[2,] -1800
[3,]  -400
[4,]  4550
[5,] 12700

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 11

$conf
          [,1]       [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  2645.898 -5369.9417 -2047.899 -7984.635 -3104.198 -1611.926 -7396.086
[2,] 10354.102   969.9417  7347.899  4284.635  2004.198  8011.926  3596.086
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -7902.097 -5502.764 -6017.078 -10921.18 -3425.063
[2,]  6602.097 14702.764  6617.078   9421.18  2625.063

$out
[1]  38100 -37900 -17700  26700  11100 -36400 -17900

$group
[1]  1  2  3  3  5  7 12

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

Warning message:
In bxp(list(stats = c(-1100, 2500, 6500, 10950, 16600, -13600, -6100,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/60yvy1272035994.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]  [,10]
[1,] 175900 182900 197700 203300 215300 230200 260400 283000 291800 292200
[2,] 178550 189200 200100 207000 223250 235300 262650 288250 297500 303900
[3,] 180700 192650 203550 209150 227000 241900 269250 289450 303450 309900
[4,] 183050 197500 208200 211200 229700 251150 281200 294000 309100 315700
[5,] 186500 202100 215100 211700 237800 268300 289600 299600 317300 329400
      [,11]  [,12] [,13]
[1,] 263100 245200    NA
[2,] 279300 259200    NA
[3,] 288850 269900    NA
[4,] 300300 274700    NA
[5,] 314300 290600    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 178647.5 188864.3 199855.5 207234.4 224058.1 234670.7 260789.2 286827.4
[2,] 182752.5 196435.7 207244.5 211065.6 229941.9 249129.3 277710.8 292072.6
         [,9]    [,10]    [,11]    [,12] [,13]
[1,] 298159.2 304517.9 279271.8 262830.3    NA
[2,] 308740.8 315282.1 298428.2 276969.7    NA

$out
[1] 211500 228700 265000 279600 284400

$group
[1]  2  4  5  8 10

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

Warning message:
In bxp(list(stats = c(175900, 178550, 180700, 183050, 186500, 182900,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7t7ci1272035994.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,] 243200.0 237700 245112.5
[2,] 247495.8 250425 247431.2
[3,] 250000.0 252875 249450.0
[4,] 251650.0 260550 250768.8
[5,] 254116.7 269150 254700.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 248105.3 248256.9 247927.7
[2,] 251894.7 257493.1 250972.3

$out
[1] 241787.5

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(243200, 247495.833333333, 250000, 251650,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1qped1272035994.ps tmp/1qped1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/2qped1272035994.ps tmp/2qped1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/3qped1272035994.ps tmp/3qped1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/40yvy1272035994.ps tmp/40yvy1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/50yvy1272035994.ps tmp/50yvy1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/60yvy1272035994.ps tmp/60yvy1272035994.png",intern=TRUE))
character(0)
> try(system("convert tmp/7t7ci1272035994.ps tmp/7t7ci1272035994.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.212   0.927   2.348