R version 2.15.2 (2012-10-26) -- "Trick or Treat"
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(235.1
+ ,280.7
+ ,264.6
+ ,240.7
+ ,201.4
+ ,240.8
+ ,241.1
+ ,223.8
+ ,206.1
+ ,174.7
+ ,203.3
+ ,220.5
+ ,299.5
+ ,347.4
+ ,338.3
+ ,327.7
+ ,351.6
+ ,396.6
+ ,438.8
+ ,395.6
+ ,363.5
+ ,378.8
+ ,357
+ ,369
+ ,464.8
+ ,479.1
+ ,431.3
+ ,366.5
+ ,326.3
+ ,355.1
+ ,331.6
+ ,261.3
+ ,249
+ ,205.5
+ ,235.6
+ ,240.9
+ ,264.9
+ ,253.8
+ ,232.3
+ ,193.8
+ ,177
+ ,213.2
+ ,207.2
+ ,180.6
+ ,188.6
+ ,175.4
+ ,199
+ ,179.6
+ ,225.8
+ ,234
+ ,200.2
+ ,183.6
+ ,178.2
+ ,203.2
+ ,208.5
+ ,191.8
+ ,172.8
+ ,148
+ ,159.4
+ ,154.5
+ ,213.2
+ ,196.4
+ ,182.8
+ ,176.4
+ ,153.6
+ ,173.2
+ ,171
+ ,151.2
+ ,161.9
+ ,157.2
+ ,201.7
+ ,236.4
+ ,356.1
+ ,398.3
+ ,403.7
+ ,384.6
+ ,365.8
+ ,368.1
+ ,367.9
+ ,347
+ ,343.3
+ ,292.9
+ ,311.5
+ ,300.9
+ ,366.9
+ ,356.9
+ ,329.7
+ ,316.2
+ ,269
+ ,289.3
+ ,266.2
+ ,253.6
+ ,233.8
+ ,228.4
+ ,253.6
+ ,260.1
+ ,306.6
+ ,309.2
+ ,309.5
+ ,271
+ ,279.9
+ ,317.9
+ ,298.4
+ ,246.7
+ ,227.3
+ ,209.1
+ ,259.9
+ ,266
+ ,320.6
+ ,308.5
+ ,282.2
+ ,262.7
+ ,263.5
+ ,313.1
+ ,284.3
+ ,252.6
+ ,250.3
+ ,246.5
+ ,312.7
+ ,333.2
+ ,446.4
+ ,511.6
+ ,515.5
+ ,506.4
+ ,483.2
+ ,522.3
+ ,509.8
+ ,460.7
+ ,405.8
+ ,375
+ ,378.5
+ ,406.8
+ ,467.8
+ ,469.8
+ ,429.8
+ ,355.8
+ ,332.7
+ ,378
+ ,360.5
+ ,334.7
+ ,319.5
+ ,323.1
+ ,363.6
+ ,352.1
+ ,411.9
+ ,388.6
+ ,416.4
+ ,360.7
+ ,338
+ ,417.2
+ ,388.4
+ ,371.1
+ ,331.5
+ ,353.7
+ ,396.7
+ ,447
+ ,533.5
+ ,565.4
+ ,542.3
+ ,488.7
+ ,467.1
+ ,531.3
+ ,496.1
+ ,444
+ ,403.4
+ ,386.3
+ ,394.1
+ ,404.1
+ ,462.1
+ ,448.1
+ ,432.3
+ ,386.3
+ ,395.2
+ ,421.9
+ ,382.9
+ ,384.2
+ ,345.5
+ ,323.4
+ ,372.6
+ ,376
+ ,462.7
+ ,487
+ ,444.2
+ ,399.3
+ ,394.9
+ ,455.4
+ ,414
+ ,375.5
+ ,347
+ ,339.4
+ ,385.8
+ ,378.8
+ ,451.8
+ ,446.1
+ ,422.5
+ ,383.1
+ ,352.8
+ ,445.3
+ ,367.5
+ ,355.1
+ ,326.2
+ ,319.8
+ ,331.8
+ ,340.9
+ ,394.1
+ ,417.2
+ ,369.9
+ ,349.2
+ ,321.4
+ ,405.7
+ ,342.9
+ ,316.5
+ ,284.2
+ ,270.9
+ ,288.8
+ ,278.8
+ ,324.4
+ ,310.9
+ ,299
+ ,273
+ ,279.3
+ ,359.2
+ ,305
+ ,282.1
+ ,250.3
+ ,246.5
+ ,257.9
+ ,266.5
+ ,315.9
+ ,318.4
+ ,295.4
+ ,266.4
+ ,245.8
+ ,362.8
+ ,324.9
+ ,294.2
+ ,289.5
+ ,295.2
+ ,290.3
+ ,272
+ ,307.4
+ ,328.7
+ ,292.9
+ ,249.1
+ ,230.4
+ ,361.5
+ ,321.7
+ ,277.2
+ ,260.7
+ ,251
+ ,257.6
+ ,241.8
+ ,287.5
+ ,292.3
+ ,274.7
+ ,254.2
+ ,230
+ ,339
+ ,318.2
+ ,287
+ ,295.8
+ ,284
+ ,271
+ ,262.7
+ ,340.6
+ ,379.4
+ ,373.3
+ ,355.2
+ ,338.4
+ ,466.9
+ ,451
+ ,422
+ ,429.2
+ ,425.9
+ ,460.7
+ ,463.6
+ ,541.4
+ ,544.2
+ ,517.5
+ ,469.4
+ ,439.4
+ ,549
+ ,533
+ ,506.1
+ ,484
+ ,457
+ ,481.5
+ ,469.5
+ ,544.7
+ ,541.2
+ ,521.5
+ ,469.7
+ ,434.4
+ ,542.6
+ ,517.3
+ ,485.7
+ ,465.8
+ ,447
+ ,426.6
+ ,411.6
+ ,467.5
+ ,484.5
+ ,451.2
+ ,417.4
+ ,379.9
+ ,484.7
+ ,455
+ ,420.8
+ ,416.5
+ ,376.3
+ ,405.6
+ ,405.8
+ ,500.8
+ ,514
+ ,475.5
+ ,430.1
+ ,414.4
+ ,538
+ ,526
+ ,488.5
+ ,520.2
+ ,504.4
+ ,568.5
+ ,610.6
+ ,818
+ ,830.9
+ ,835.9
+ ,782
+ ,762.3
+ ,856.9
+ ,820.9
+ ,769.6
+ ,752.2
+ ,724.4
+ ,723.1
+ ,719.5
+ ,817.4
+ ,803.3
+ ,752.5
+ ,689
+ ,630.4
+ ,765.5
+ ,757.7
+ ,732.2
+ ,702.6
+ ,683.3
+ ,709.5
+ ,702.2
+ ,784.8
+ ,810.9
+ ,755.6
+ ,656.8
+ ,615.1
+ ,745.3
+ ,694.1
+ ,675.7
+ ,643.7
+ ,622.1
+ ,634.6
+ ,588
+ ,689.7
+ ,673.9
+ ,647.9
+ ,568.8
+ ,545.7
+ ,632.6
+ ,643.8
+ ,593.1
+ ,579.7
+ ,546
+ ,562.9
+ ,572.5)
> par1 = '12'
> par1 <- '12'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Mean Plot (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanplot.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 372
> (np <- floor(n / par1))
[1] 31
> 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] 31 31 31 31 31 31 31 31 31 31 31 31
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 235.1 299.5 464.8 264.9 225.8 213.2 356.1 366.9 306.6 320.6 446.4 467.8
 [2,] 280.7 347.4 479.1 253.8 234.0 196.4 398.3 356.9 309.2 308.5 511.6 469.8
 [3,] 264.6 338.3 431.3 232.3 200.2 182.8 403.7 329.7 309.5 282.2 515.5 429.8
 [4,] 240.7 327.7 366.5 193.8 183.6 176.4 384.6 316.2 271.0 262.7 506.4 355.8
 [5,] 201.4 351.6 326.3 177.0 178.2 153.6 365.8 269.0 279.9 263.5 483.2 332.7
 [6,] 240.8 396.6 355.1 213.2 203.2 173.2 368.1 289.3 317.9 313.1 522.3 378.0
 [7,] 241.1 438.8 331.6 207.2 208.5 171.0 367.9 266.2 298.4 284.3 509.8 360.5
 [8,] 223.8 395.6 261.3 180.6 191.8 151.2 347.0 253.6 246.7 252.6 460.7 334.7
 [9,] 206.1 363.5 249.0 188.6 172.8 161.9 343.3 233.8 227.3 250.3 405.8 319.5
[10,] 174.7 378.8 205.5 175.4 148.0 157.2 292.9 228.4 209.1 246.5 375.0 323.1
[11,] 203.3 357.0 235.6 199.0 159.4 201.7 311.5 253.6 259.9 312.7 378.5 363.6
[12,] 220.5 369.0 240.9 179.6 154.5 236.4 300.9 260.1 266.0 333.2 406.8 352.1
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] 411.9 533.5 462.1 462.7 451.8 394.1 324.4 315.9 307.4 287.5 340.6 541.4
 [2,] 388.6 565.4 448.1 487.0 446.1 417.2 310.9 318.4 328.7 292.3 379.4 544.2
 [3,] 416.4 542.3 432.3 444.2 422.5 369.9 299.0 295.4 292.9 274.7 373.3 517.5
 [4,] 360.7 488.7 386.3 399.3 383.1 349.2 273.0 266.4 249.1 254.2 355.2 469.4
 [5,] 338.0 467.1 395.2 394.9 352.8 321.4 279.3 245.8 230.4 230.0 338.4 439.4
 [6,] 417.2 531.3 421.9 455.4 445.3 405.7 359.2 362.8 361.5 339.0 466.9 549.0
 [7,] 388.4 496.1 382.9 414.0 367.5 342.9 305.0 324.9 321.7 318.2 451.0 533.0
 [8,] 371.1 444.0 384.2 375.5 355.1 316.5 282.1 294.2 277.2 287.0 422.0 506.1
 [9,] 331.5 403.4 345.5 347.0 326.2 284.2 250.3 289.5 260.7 295.8 429.2 484.0
[10,] 353.7 386.3 323.4 339.4 319.8 270.9 246.5 295.2 251.0 284.0 425.9 457.0
[11,] 396.7 394.1 372.6 385.8 331.8 288.8 257.9 290.3 257.6 271.0 460.7 481.5
[12,] 447.0 404.1 376.0 378.8 340.9 278.8 266.5 272.0 241.8 262.7 463.6 469.5
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32]
 [1,] 544.7 467.5 500.8 818.0 817.4 784.8 689.7    NA
 [2,] 541.2 484.5 514.0 830.9 803.3 810.9 673.9    NA
 [3,] 521.5 451.2 475.5 835.9 752.5 755.6 647.9    NA
 [4,] 469.7 417.4 430.1 782.0 689.0 656.8 568.8    NA
 [5,] 434.4 379.9 414.4 762.3 630.4 615.1 545.7    NA
 [6,] 542.6 484.7 538.0 856.9 765.5 745.3 632.6    NA
 [7,] 517.3 455.0 526.0 820.9 757.7 694.1 643.8    NA
 [8,] 485.7 420.8 488.5 769.6 732.2 675.7 593.1    NA
 [9,] 465.8 416.5 520.2 752.2 702.6 643.7 579.7    NA
[10,] 447.0 376.3 504.4 724.4 683.3 622.1 546.0    NA
[11,] 426.6 405.6 568.5 723.1 709.5 634.6 562.9    NA
[12,] 411.6 405.8 610.6 719.5 702.2 588.0 572.5    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]  45.6  47.9  14.3 -11.1   8.2 -16.8  42.2 -10.0   2.6 -12.1  65.2   2.0
 [2,] -16.1  -9.1 -47.8 -21.5 -33.8 -13.6   5.4 -27.2   0.3 -26.3   3.9 -40.0
 [3,] -23.9 -10.6 -64.8 -38.5 -16.6  -6.4 -19.1 -13.5 -38.5 -19.5  -9.1 -74.0
 [4,] -39.3  23.9 -40.2 -16.8  -5.4 -22.8 -18.8 -47.2   8.9   0.8 -23.2 -23.1
 [5,]  39.4  45.0  28.8  36.2  25.0  19.6   2.3  20.3  38.0  49.6  39.1  45.3
 [6,]   0.3  42.2 -23.5  -6.0   5.3  -2.2  -0.2 -23.1 -19.5 -28.8 -12.5 -17.5
 [7,] -17.3 -43.2 -70.3 -26.6 -16.7 -19.8 -20.9 -12.6 -51.7 -31.7 -49.1 -25.8
 [8,] -17.7 -32.1 -12.3   8.0 -19.0  10.7  -3.7 -19.8 -19.4  -2.3 -54.9 -15.2
 [9,] -31.4  15.3 -43.5 -13.2 -24.8  -4.7 -50.4  -5.4 -18.2  -3.8 -30.8   3.6
[10,]  28.6 -21.8  30.1  23.6  11.4  44.5  18.6  25.2  50.8  66.2   3.5  40.5
[11,]  17.2  12.0   5.3 -19.4  -4.9  34.7 -10.6   6.5   6.1  20.5  28.3 -11.5
[12,]  79.0  95.8  24.0  46.2  58.7 119.7  66.0  46.5  54.6 113.2  61.0  59.8
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] -23.3  31.9 -14.0  24.3  -5.7  23.1 -13.5   2.5  21.3   4.8  38.8   2.8
 [2,]  27.8 -23.1 -15.8 -42.8 -23.6 -47.3 -11.9 -23.0 -35.8 -17.6  -6.1 -26.7
 [3,] -55.7 -53.6 -46.0 -44.9 -39.4 -20.7 -26.0 -29.0 -43.8 -20.5 -18.1 -48.1
 [4,] -22.7 -21.6   8.9  -4.4 -30.3 -27.8   6.3 -20.6 -18.7 -24.2 -16.8 -30.0
 [5,]  79.2  64.2  26.7  60.5  92.5  84.3  79.9 117.0 131.1 109.0 128.5 109.6
 [6,] -28.8 -35.2 -39.0 -41.4 -77.8 -62.8 -54.2 -37.9 -39.8 -20.8 -15.9 -16.0
 [7,] -17.3 -52.1   1.3 -38.5 -12.4 -26.4 -22.9 -30.7 -44.5 -31.2 -29.0 -26.9
 [8,] -39.6 -40.6 -38.7 -28.5 -28.9 -32.3 -31.8  -4.7 -16.5   8.8   7.2 -22.1
 [9,]  22.2 -17.1 -22.1  -7.6  -6.4 -13.3  -3.8   5.7  -9.7 -11.8  -3.3 -27.0
[10,]  43.0   7.8  49.2  46.4  12.0  17.9  11.4  -4.9   6.6 -13.0  34.8  24.5
[11,]  50.3  10.0   3.4  -7.0   9.1 -10.0   8.6 -18.3 -15.8  -8.3   2.9 -12.0
[12,]  86.5  58.0  86.7  73.0  53.2  45.6  49.4  35.4  45.7  77.9  77.8  75.2
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32]
 [1,]  -3.5  17.0  13.2  12.9 -14.1  26.1 -15.8    NA
 [2,] -19.7 -33.3 -38.5   5.0 -50.8 -55.3 -26.0    NA
 [3,] -51.8 -33.8 -45.4 -53.9 -63.5 -98.8 -79.1    NA
 [4,] -35.3 -37.5 -15.7 -19.7 -58.6 -41.7 -23.1    NA
 [5,] 108.2 104.8 123.6  94.6 135.1 130.2  86.9    NA
 [6,] -25.3 -29.7 -12.0 -36.0  -7.8 -51.2  11.2    NA
 [7,] -31.6 -34.2 -37.5 -51.3 -25.5 -18.4 -50.7    NA
 [8,] -19.9  -4.3  31.7 -17.4 -29.6 -32.0 -13.4    NA
 [9,] -18.8 -40.2 -15.8 -27.8 -19.3 -21.6 -33.7    NA
[10,] -20.4  29.3  64.1  -1.3  26.2  12.5  16.9    NA
[11,] -15.0   0.2  42.1  -3.6  -7.3 -46.6   9.6    NA
[12,]  55.9  95.0 207.4  97.9  82.6 101.7    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/fisher/rcomp/tmp/1k9fc1355234913.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/fisher/rcomp/tmp/2k7bp1355234913.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/fisher/rcomp/tmp/3a53s1355234913.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/fisher/rcomp/tmp/4mhsr1355234913.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,] 213.20 196.40 182.8 176.40 153.60 173.20 171.00 151.20 161.90 148.0 159.4
[2,] 311.65 314.65 297.2 268.70 266.25 347.05 311.60 269.25 250.30 246.5 258.9
[3,] 411.90 417.20 416.4 360.70 338.40 405.70 367.90 355.10 331.50 323.1 357.0
[4,] 484.30 512.80 495.5 449.75 424.40 526.80 502.95 452.35 422.85 406.1 416.1
[5,] 689.70 803.30 755.6 689.00 630.40 765.50 757.70 675.70 643.70 622.1 634.6
      [,12]
[1,] 154.50
[2,] 264.35
[3,] 352.10
[4,] 429.30
[5,] 610.60

$n
 [1] 31 31 31 31 31 31 31 31 31 31 31 31

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 362.906 360.9697 360.1271 309.3223 293.5208 354.6912 313.5994 303.1405
[2,] 460.894 473.4303 472.6729 412.0777 383.2792 456.7088 422.2006 407.0595
         [,9]    [,10]    [,11]    [,12]
[1,] 282.5344 277.8093 312.3904 305.2911
[2,] 380.4656 368.3907 401.6096 398.9089

$out
 [1] 818.0 817.4 784.8 830.9 810.9 835.9 782.0 762.3 856.9 820.9 769.6 732.2
[13] 752.2 702.6 724.4 683.3 723.1 709.5 719.5 702.2

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

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

> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/5uzqa1355234913.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]
[1,] -23.30 -55.30 -98.8 -47.20   2.30 -62.80 -70.30 -54.9 -50.40 -21.80 -19.4
[2,] -10.55 -34.80 -52.7 -30.15  38.55 -36.95 -40.85 -30.7 -25.90   9.60 -10.3
[3,]   4.80 -23.10 -38.5 -22.70  79.20 -23.10 -29.00 -19.0 -15.80  23.60   2.9
[4,]  23.70 -12.75 -20.0 -16.25 108.60  -9.90 -20.35  -4.5  -5.05  37.65   9.8
[5,]  65.20   5.40  -6.4   0.80 135.10  11.20   1.30  31.7  22.20  66.20  34.7
     [,12]
[1,]  24.0
[2,]  53.2
[3,]  69.5
[4,]  86.7
[5,] 119.7

$n
 [1] 31 31 31 31 31 31 31 31 31 31 31 30

$conf
          [,1]      [,2]      [,3]      [,4]     [,5]      [,6]      [,7]
[1,] -4.919341 -29.35727 -47.77949 -26.64449 59.32146 -30.77615 -34.81742
[2,] 14.519341 -16.84273 -29.22051 -18.75551 99.07854 -15.42385 -23.18258
          [,8]       [,9]    [,10]     [,11]    [,12]
[1,] -26.43494 -21.716737 15.64007 -2.803905 59.83635
[2,] -11.56506  -9.883263 31.55993  8.603905 79.16365

$out
 [1]  27.8  23.9   8.9   8.9   6.3 -58.6  42.2 -77.8  50.3  42.1 -46.6 207.4

$group
 [1]  2  4  4  4  4  4  6  6 11 11 11 12

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

> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/6kklk1355234913.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,] 174.70 299.50 205.50 175.40 148.00 151.20 292.90 228.40 209.10 246.50
[2,] 204.70 342.85 244.95 180.10 166.10 159.55 327.40 253.60 253.30 257.65
[3,] 229.45 360.25 328.95 196.40 187.70 174.80 360.95 267.60 275.45 283.25
[4,] 240.95 387.20 398.90 222.75 205.85 199.05 376.35 322.95 307.90 312.90
[5,] 280.70 438.80 479.10 264.90 234.00 236.40 403.70 366.90 317.90 333.20
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 375.00 319.50 331.50 386.30 323.40 339.40 319.80 270.90 246.50 245.80
[2,] 406.30 333.70 357.20 403.75 374.30 377.15 336.35 286.50 262.20 280.75
[3,] 471.95 358.15 388.50 477.90 385.25 397.10 361.30 332.15 280.70 294.70
[4,] 510.70 403.90 414.15 532.40 427.10 449.80 433.90 382.00 307.95 317.15
[5,] 522.30 469.80 447.00 565.40 462.10 487.00 451.80 417.20 359.20 362.80
      [,21]  [,22]  [,23]  [,24]  [,25]  [,26] [,27] [,28]  [,29]  [,30]  [,31]
[1,] 230.40 230.00 338.40 439.40 411.60 376.30 414.4 719.5 630.40 588.00 545.70
[2,] 250.05 266.85 364.25 469.45 440.70 405.70 482.0 738.3 695.60 628.35 565.85
[3,] 268.95 285.50 423.95 495.05 477.70 419.10 509.2 775.8 720.85 666.25 586.40
[4,] 314.55 294.05 455.85 537.20 531.35 461.25 532.0 825.9 761.60 750.45 645.85
[5,] 361.50 318.20 466.90 549.00 544.70 484.70 568.5 856.9 817.40 810.90 689.70
     [,32]
[1,]    NA
[2,]    NA
[3,]    NA
[4,]    NA
[5,]    NA

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 212.9161 340.0217 258.7324 176.947 169.5698 156.7838 338.6236 235.969
[2,] 245.9839 380.4783 399.1676 215.853 205.8302 192.8162 383.2764 299.231
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 250.5466 258.0501 424.3325 326.1313 362.5247 419.2219 361.1676 363.9638
[2,] 300.3534 308.4499 519.5675 390.1687 414.4753 536.5781 409.3324 430.2362
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 316.8068 288.5918 259.8331 278.0977 239.5311 273.0939 382.1706 464.1488
[2,] 405.7932 375.7082 301.5669 311.3023 298.3689 297.9061 465.7294 525.9512
        [,25]    [,26]    [,27]    [,28]   [,29]    [,30]    [,31] [,32]
[1,] 436.3539 393.7633 486.3947 735.8451 690.747 610.5594 549.9115    NA
[2,] 519.0461 444.4367 532.0053 815.7549 750.953 721.9406 622.8885    NA

$out
[1] 339.0 610.6

$group
[1] 22 27

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30"
[31] "31" NA  

Warning message:
In bxp(list(stats = c(174.7, 204.7, 229.45, 240.95, 280.7, 299.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/7tux41355234913.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,] 347.4581 323.10 326.3000
[2,] 366.2145 345.25 341.4125
[3,] 380.8710 358.85 360.0125
[4,] 426.8435 408.80 402.6250
[5,] 442.9258 417.20 436.9250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 353.2177 329.8644 332.0931
[2,] 408.5243 387.8356 387.9319

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(347.458064516129, 366.214516129032, 380.870967741935,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1k9fc1355234913.ps tmp/1k9fc1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/2k7bp1355234913.ps tmp/2k7bp1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/3a53s1355234913.ps tmp/3a53s1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/4mhsr1355234913.ps tmp/4mhsr1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/5uzqa1355234913.ps tmp/5uzqa1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/6kklk1355234913.ps tmp/6kklk1355234913.png",intern=TRUE))
character(0)
> try(system("convert tmp/7tux41355234913.ps tmp/7tux41355234913.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.598   0.985   4.574