R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)

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(265.2
+ ,261.3
+ ,253.7
+ ,249.2
+ ,239.1
+ ,236.4
+ ,235.2
+ ,245.2
+ ,246.2
+ ,247.7
+ ,251.4
+ ,253.3
+ ,254.8
+ ,250
+ ,249.3
+ ,241.5
+ ,243.3
+ ,248
+ ,253
+ ,252.9
+ ,251.5
+ ,251.6
+ ,253.5
+ ,259.8
+ ,334.1
+ ,448
+ ,445.8
+ ,445
+ ,448.2
+ ,438.2
+ ,439.8
+ ,423.4
+ ,410.8
+ ,408.4
+ ,406.7
+ ,405.9
+ ,402.7
+ ,405.1
+ ,399.6
+ ,386.5
+ ,381.4
+ ,375.2
+ ,357.7
+ ,359
+ ,355
+ ,352.7
+ ,344.4
+ ,343.8
+ ,338
+ ,339
+ ,333.3
+ ,334.4
+ ,328.3
+ ,330.7
+ ,330
+ ,331.6
+ ,351.2
+ ,389.4
+ ,410.9
+ ,442.8
+ ,462.8
+ ,466.9
+ ,461.7
+ ,439.2
+ ,430.3
+ ,416.1
+ ,402.5
+ ,397.3
+ ,403.3
+ ,395.9
+ ,387.8
+ ,378.6
+ ,377.1
+ ,370.4
+ ,362
+ ,350.3
+ ,348.2
+ ,344.6
+ ,343.5
+ ,342.8
+ ,347.6
+ ,346.6
+ ,349.5
+ ,342.1
+ ,342
+ ,342.8
+ ,339.3
+ ,348.2
+ ,333.7
+ ,334.7
+ ,354
+ ,367.7
+ ,363.3
+ ,358.4
+ ,353.1
+ ,343.1
+ ,344.6
+ ,344.4
+ ,333.9
+ ,331.7
+ ,324.3
+ ,321.2
+ ,322.4
+ ,321.7
+ ,320.5
+ ,312.8
+ ,309.7
+ ,315.6
+ ,309.7
+ ,304.6
+ ,302.5
+ ,301.5
+ ,298.8
+ ,291.3
+ ,293.6
+ ,294.6
+ ,285.9
+ ,297.6
+ ,301.1
+ ,293.8
+ ,297.7
+ ,292.9
+ ,292.1
+ ,287.2
+ ,288.2
+ ,283.8
+ ,299.9
+ ,292.4
+ ,293.3
+ ,300.8
+ ,293.7
+ ,293.1
+ ,294.4
+ ,292.1
+ ,291.9
+ ,282.5
+ ,277.9
+ ,287.5
+ ,289.2
+ ,285.6
+ ,293.2
+ ,290.8
+ ,283.1
+ ,275
+ ,287.8
+ ,287.8
+ ,287.4
+ ,284
+ ,277.8
+ ,277.6
+ ,304.9
+ ,294
+ ,300.9
+ ,324
+ ,332.9
+ ,341.6
+ ,333.4
+ ,348.2
+ ,344.7
+ ,344.7
+ ,329.3
+ ,323.5
+ ,323.2
+ ,317.4
+ ,330.1
+ ,329.2
+ ,334.9
+ ,315.8
+ ,315.4
+ ,319.6
+ ,317.3
+ ,313.8
+ ,315.8
+ ,311.3)
> par1 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 174
> (np <- floor(n / par1))
[1] 14
> 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] 15 15 15 15 15 15 14 14 14 14 14 14
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 265.2 254.8 334.1 402.7 338.0 462.8 377.1 342.0 344.6 309.7 297.7 294.4
 [2,] 261.3 250.0 448.0 405.1 339.0 466.9 370.4 342.8 344.4 304.6 292.9 292.1
 [3,] 253.7 249.3 445.8 399.6 333.3 461.7 362.0 339.3 333.9 302.5 292.1 291.9
 [4,] 249.2 241.5 445.0 386.5 334.4 439.2 350.3 348.2 331.7 301.5 287.2 282.5
 [5,] 239.1 243.3 448.2 381.4 328.3 430.3 348.2 333.7 324.3 298.8 288.2 277.9
 [6,] 236.4 248.0 438.2 375.2 330.7 416.1 344.6 334.7 321.2 291.3 283.8 287.5
 [7,] 235.2 253.0 439.8 357.7 330.0 402.5 343.5 354.0 322.4 293.6 299.9 289.2
 [8,] 245.2 252.9 423.4 359.0 331.6 397.3 342.8 367.7 321.7 294.6 292.4 285.6
 [9,] 246.2 251.5 410.8 355.0 351.2 403.3 347.6 363.3 320.5 285.9 293.3 293.2
[10,] 247.7 251.6 408.4 352.7 389.4 395.9 346.6 358.4 312.8 297.6 300.8 290.8
[11,] 251.4 253.5 406.7 344.4 410.9 387.8 349.5 353.1 309.7 301.1 293.7 283.1
[12,] 253.3 259.8 405.9 343.8 442.8 378.6 342.1 343.1 315.6 293.8 293.1 275.0
      [,13] [,14] [,15]
 [1,] 287.8 333.4 315.4
 [2,] 287.8 348.2 319.6
 [3,] 287.4 344.7 317.3
 [4,] 284.0 344.7 313.8
 [5,] 277.8 329.3 315.8
 [6,] 277.6 323.5 311.3
 [7,] 304.9 323.2    NA
 [8,] 294.0 317.4    NA
 [9,] 300.9 330.1    NA
[10,] 324.0 329.2    NA
[11,] 332.9 334.9    NA
[12,] 341.6 315.8    NA
> darr
       [,1] [,2]  [,3]  [,4] [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]  -3.9 -4.8 113.9   2.4  1.0   4.1  -6.7   0.8  -0.2  -5.1  -4.8  -2.3
 [2,]  -7.6 -0.7  -2.2  -5.5 -5.7  -5.2  -8.4  -3.5 -10.5  -2.1  -0.8  -0.2
 [3,]  -4.5 -7.8  -0.8 -13.1  1.1 -22.5 -11.7   8.9  -2.2  -1.0  -4.9  -9.4
 [4,] -10.1  1.8   3.2  -5.1 -6.1  -8.9  -2.1 -14.5  -7.4  -2.7   1.0  -4.6
 [5,]  -2.7  4.7 -10.0  -6.2  2.4 -14.2  -3.6   1.0  -3.1  -7.5  -4.4   9.6
 [6,]  -1.2  5.0   1.6 -17.5 -0.7 -13.6  -1.1  19.3   1.2   2.3  16.1   1.7
 [7,]  10.0 -0.1 -16.4   1.3  1.6  -5.2  -0.7  13.7  -0.7   1.0  -7.5  -3.6
 [8,]   1.0 -1.4 -12.6  -4.0 19.6   6.0   4.8  -4.4  -1.2  -8.7   0.9   7.6
 [9,]   1.5  0.1  -2.4  -2.3 38.2  -7.4  -1.0  -4.9  -7.7  11.7   7.5  -2.4
[10,]   3.7  1.9  -1.7  -8.3 21.5  -8.1   2.9  -5.3  -3.1   3.5  -7.1  -7.7
[11,]   1.9  6.3  -0.8  -0.6 31.9  -9.2  -7.4 -10.0   5.9  -7.3  -0.6  -8.1
[12,]   1.5 74.3  -3.2  -5.8 20.0  -1.5  -0.1   1.5  -5.9   3.9   1.3  12.8
      [,13] [,14] [,15]
 [1,]   0.0  14.8   4.2
 [2,]  -0.4  -3.5  -2.3
 [3,]  -3.4   0.0  -3.5
 [4,]  -6.2 -15.4   2.0
 [5,]  -0.2  -5.8  -4.5
 [6,]  27.3  -0.3    NA
 [7,] -10.9  -5.8    NA
 [8,]   6.9  12.7    NA
 [9,]  23.1  -0.9    NA
[10,]   8.9   5.7    NA
[11,]   8.7 -19.1    NA
[12,]  -8.2  -0.4    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/wessaorg/rcomp/tmp/1niyv1350832091.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2wfqg1350832091.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3q9om1350832091.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4ozz01350832091.ps",horizontal=F,onefile=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] [,11]
[1,] 254.80 250.0 249.30 241.50 239.10 236.40 235.2 245.20 246.2 247.7 251.4
[2,] 296.05 292.5 292.00 285.60 283.05 285.65 293.6 292.40 293.2 297.6 293.7
[3,] 333.40 339.0 333.30 331.70 324.30 321.20 322.8 319.55 325.3 326.6 333.9
[4,] 343.30 359.3 353.35 349.25 340.95 339.65 354.0 359.00 355.0 358.4 353.1
[5,] 402.70 448.0 399.60 439.20 381.40 416.10 439.8 423.40 410.8 408.4 410.9
     [,12]
[1,] 253.3
[2,] 293.1
[3,] 328.7
[4,] 343.8
[5,] 405.9

$n
 [1] 15 15 15 15 15 15 14 14 14 14 14 14

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 314.1242 311.7487 308.272 305.7337 300.6794 299.1705 297.2947 291.4266
[2,] 352.6758 366.2513 358.328 357.6663 347.9206 343.2295 348.3053 347.6734
         [,9]    [,10]   [,11]    [,12]
[1,] 299.2035 300.9258 308.817 307.2908
[2,] 351.3965 352.2742 358.983 350.1092

$out
[1] 462.8 466.9 445.8 461.7 445.0 448.2 430.3 438.2 442.8

$group
[1]  1  2  3  3  4  5  5  6 12

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

Warning message:
In bxp(list(stats = c(254.8, 296.05, 333.4, 343.3, 402.7, 250, 292.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5neek1350832091.ps",horizontal=F,onefile=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,] -6.70 -10.50 -13.1 -15.40 -14.2 -1.2 -16.4 -12.60 -7.70  -8.3 -19.1  -8.2
[2,] -4.35  -5.60  -8.6  -8.15  -6.0 -1.1  -5.8  -4.00 -2.40  -7.1  -8.1  -3.2
[3,]  0.00  -3.50  -3.5  -5.10  -3.6  1.4  -0.7   0.95 -0.95   0.1  -0.7   0.6
[4,]  3.25  -1.45  -0.9  -0.55   0.4  5.0   1.3   6.90  7.50   3.7   5.9   3.9
[5,]  4.20  -0.20   8.9   3.20   9.6  5.0  10.0  19.60 11.70   8.9   8.7  12.8

$n
 [1] 15 15 15 15 15 14 14 14 14 14 14 14

$conf
          [,1]     [,2]       [,3]      [,4]       [,5]      [,6]      [,7]
[1,] -3.100452 -5.19301 -6.6412477 -8.200452 -6.2109072 -1.175864 -3.698137
[2,]  3.100452 -1.80699 -0.3587523 -1.999548 -0.9890928  3.975864  2.298137
          [,8]    [,9]     [,10]     [,11]     [,12]
[1,] -3.652773 -5.1305 -4.460546 -6.611819 -2.398137
[2,]  5.552773  3.2305  4.660546  5.211819  3.598137

$out
 [1] 113.9  14.8 -22.5 -17.5 -13.6  19.3  16.1  27.3  13.7  38.2  23.1  21.5
[13]  31.9  74.3  20.0

$group
 [1]  1  1  3  6  6  6  6  6  7  9  9 10 11 12 12

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

Warning message:
In bxp(list(stats = c(-6.70000000000005, -4.34999999999999, 0, 3.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/670al1350832091.ps",horizontal=F,onefile=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,] 235.20 243.30 405.90 343.80 328.30 378.60 342.10 333.70 309.70 285.9
[2,] 242.15 248.65 407.55 353.85 331.15 396.60 344.05 340.65 318.05 293.7
[3,] 248.45 251.55 430.80 367.10 336.20 409.70 347.90 345.65 322.05 298.2
[4,] 253.50 253.25 445.40 393.05 370.30 450.45 356.15 356.20 332.80 302.0
[5,] 265.20 259.80 448.20 405.10 410.90 466.90 370.40 367.70 344.60 309.7
      [,11]  [,12]  [,13]  [,14] [,15]
[1,] 283.80 275.00 277.60 315.80 311.3
[2,] 290.15 282.80 285.70 323.35 313.8
[3,] 293.00 288.35 290.90 329.70 315.6
[4,] 295.70 292.00 314.45 339.80 317.3
[5,] 300.80 294.40 341.60 348.20 319.6

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 243.2732 249.4519 413.5364 349.2206 318.3434 385.1387 342.3811 338.5575
[2,] 253.6268 253.6481 448.0636 384.9794 354.0566 434.2613 353.4189 352.7425
         [,9]    [,10]    [,11]    [,12]    [,13]   [,14]    [,15]
[1,] 315.3224 294.4143 290.4686 284.1538 277.7869 322.197 313.3424
[2,] 328.7776 301.9857 295.5314 292.5462 304.0131 337.203 317.8576

$out
[1] 241.5 334.1 442.8 377.1

$group
[1] 2 3 5 7

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

Warning message:
In bxp(list(stats = c(235.2, 242.15, 248.45, 253.5, 265.2, 243.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/79hj31350832091.ps",horizontal=F,onefile=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,] 321.3400 319.55 312.0000
[2,] 324.6140 323.55 317.9375
[3,] 328.9357 327.65 322.9750
[4,] 330.0626 333.35 323.8750
[5,] 338.2067 339.00 327.6875

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 326.4506 323.1802 320.2669
[2,] 331.4208 332.1198 325.6831

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(321.34, 324.614047619048, 328.935714285714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1niyv1350832091.ps tmp/1niyv1350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/2wfqg1350832091.ps tmp/2wfqg1350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/3q9om1350832091.ps tmp/3q9om1350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/4ozz01350832091.ps tmp/4ozz01350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/5neek1350832091.ps tmp/5neek1350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/670al1350832091.ps tmp/670al1350832091.png",intern=TRUE))
character(0)
> try(system("convert tmp/79hj31350832091.ps tmp/79hj31350832091.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.453   0.606   4.068