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(192.1
+ ,201.8
+ ,186.9
+ ,218
+ ,214.4
+ ,227.5
+ ,204.1
+ ,225.8
+ ,223.7
+ ,244.7
+ ,243.9
+ ,257.3
+ ,234.5
+ ,251.4
+ ,243.8
+ ,247.4
+ ,245.3
+ ,262.5
+ ,270
+ ,259.9
+ ,262.2
+ ,244.9
+ ,249.3
+ ,268.2
+ ,231.2
+ ,264.3
+ ,252.7
+ ,275.5
+ ,261.5
+ ,275.5
+ ,272.3
+ ,268.6
+ ,270.4
+ ,267.7
+ ,275
+ ,272.6
+ ,248.6
+ ,279.4
+ ,270.5
+ ,292.8
+ ,297.8
+ ,296.8
+ ,290.9
+ ,282.8
+ ,312.8
+ ,303.2
+ ,301.4
+ ,289.8
+ ,279.6
+ ,302.2
+ ,299.1
+ ,319.7
+ ,310.9
+ ,315.2
+ ,338.5
+ ,315.6
+ ,321.2
+ ,318.5
+ ,342.7
+ ,261.4
+ ,287
+ ,331.5
+ ,326.9
+ ,338.6
+ ,337
+ ,358.4
+ ,344.5
+ ,345.7
+ ,344.1
+ ,317.4
+ ,354.5
+ ,345.2
+ ,314.1
+ ,352.5
+ ,361.2
+ ,365.9
+ ,332.5
+ ,364
+ ,359.1
+ ,345.6
+ ,366.9
+ ,370.2
+ ,359.9
+ ,366.6
+ ,336.3
+ ,368.5
+ ,374.2
+ ,384.3
+ ,358.9
+ ,407.7
+ ,433.3
+ ,404.7
+ ,392.7
+ ,409.7
+ ,416.5
+ ,414.3
+ ,404.3
+ ,421.4
+ ,372.6
+ ,404.7
+ ,420.2
+ ,438.4
+ ,449.1
+ ,445.8
+ ,413.8
+ ,420.5
+ ,442.3
+ ,438.9
+ ,394.5
+ ,416.8
+ ,402.9
+ ,424.5
+ ,432.3
+ ,484.1
+ ,492.7
+ ,496.3
+ ,471.9
+ ,491.2
+ ,512.9
+ ,482.4
+ ,407.9
+ ,448.5
+ ,431.1
+ ,498.8
+ ,497.1
+ ,517.1
+ ,487.7
+ ,512.5
+ ,550.1
+ ,532.5
+ ,524.1
+ ,515.7
+ ,461
+ ,529.3
+ ,467.4
+ ,559.8
+ ,536.5
+ ,531.9
+ ,546.5
+ ,547.4
+ ,536.1
+ ,482.8
+ ,551
+ ,532.9
+ ,484.1
+ ,554.8
+ ,537
+ ,558
+ ,511.4
+ ,502.9
+ ,558.6
+ ,545.1
+ ,574.3
+ ,542.2
+ ,600
+ ,588.6
+ ,524.4
+ ,618.5
+ ,580.9
+ ,557.2
+ ,571.2
+ ,597.5
+ ,601.7
+ ,558.9
+ ,600.9
+ ,601
+ ,615.7
+ ,578.1
+ ,495.9
+ ,526.8
+ ,522.1
+ ,605.1
+ ,574.4
+ ,609.7
+ ,580.7
+ ,565.1
+ ,590.7
+ ,571.5
+ ,601.3
+ ,567.3
+ ,467.9
+ ,588.9
+ ,579.4
+ ,502.6
+ ,568.7
+ ,616
+ ,586.2
+ ,575.5
+ ,599.9
+ ,568.2
+ ,516
+ ,493.4
+ ,496.8
+ ,529.9
+ ,491.7
+ ,543.2
+ ,490.8
+ ,554.7
+ ,625.7
+ ,605
+ ,645.2
+ ,645.2
+ ,611.8
+ ,600.3
+ ,549.8
+ ,635.5
+ ,617.7
+ ,643.5
+ ,485.7
+ ,689.5
+ ,692
+ ,677.3
+ ,704.7
+ ,668.6
+ ,717.8
+ ,689.8
+ ,640.4
+ ,675.2
+ ,528.1
+ ,538
+ ,527.2
+ ,655.6
+ ,650.6
+ ,623.7
+ ,748.4
+ ,727.4
+ ,750.5
+ ,678.9
+ ,659.5
+ ,691.9
+ ,639.8
+ ,663.8
+ ,572.9
+ ,592.5
+ ,734.8
+ ,696.1
+ ,589.2
+ ,662.9
+ ,661.2
+ ,672.1
+ ,583.7
+ ,705.5
+ ,631
+ ,733.3
+ ,674.9
+ ,695.5
+ ,634.1
+ ,630.6
+ ,635.2
+ ,554.1
+ ,623.9
+ ,679.3
+ ,565.6
+ ,564.1
+ ,637.2
+ ,650.8
+ ,602.7
+ ,587.5
+ ,619.2
+ ,616.5
+ ,637.9
+ ,557.9
+ ,594
+ ,668.7
+ ,603.3
+ ,674.5
+ ,573.4
+ ,706
+ ,693.7
+ ,627.5
+ ,550.7
+ ,592.3
+ ,660.2
+ ,597.3
+ ,641
+ ,663.6
+ ,595.9
+ ,638.4
+ ,665.4
+ ,671.4
+ ,637
+ ,685.7
+ ,705.8
+ ,704.8
+ ,734.4
+ ,674.2
+ ,748.6
+ ,763.4
+ ,658
+ ,627.5
+ ,528.9
+ ,488.3
+ ,575.5
+ ,735.6
+ ,685.3
+ ,613.6
+ ,629.5
+ ,634.7
+ ,652.6
+ ,728.3
+ ,634.3
+ ,690.7
+ ,676.3
+ ,675.4
+ ,595.6
+ ,712.4
+ ,735.8
+ ,544.4
+ ,567
+ ,510
+ ,564
+ ,630.7
+ ,496.7
+ ,660.9
+ ,601.2
+ ,655.2
+ ,591.6
+ ,606.1
+ ,560.7
+ ,368.3
+ ,371.6
+ ,413.9
+ ,413.9
+ ,389
+ ,399.2
+ ,429.8
+ ,395.6
+ ,472
+ ,486
+ ,525
+ ,396
+ ,511
+ ,525
+ ,492
+ ,517
+ ,525
+ ,474
+ ,539
+ ,468
+ ,543
+ ,532
+ ,565
+ ,535
+ ,534
+ ,546
+ ,494
+ ,552
+ ,511
+ ,451
+ ,537
+ ,494
+ ,549
+ ,544
+ ,598
+ ,583
+ ,582
+ ,589
+ ,578
+ ,561
+ ,592
+ ,504
+ ,545
+ ,547
+ ,585
+ ,562
+ ,520
+ ,581
+ ,590
+ ,562
+ ,548
+ ,567
+ ,542
+ ,473
+ ,531
+ ,462
+ ,479
+ ,533
+ ,552
+ ,547
+ ,562
+ ,524
+ ,479
+ ,445
+ ,406
+ ,475
+ ,589
+ ,495
+ ,484
+ ,536
+ ,555
+ ,565
+ ,564
+ ,573
+ ,569
+ ,588
+ ,546
+ ,508
+ ,560
+ ,558
+ ,516
+ ,549
+ ,595
+ ,586
+ ,597
+ ,592
+ ,538
+ ,590
+ ,576
+ ,451
+ ,538
+ ,555
+ ,532
+ ,530
+ ,553
+ ,626
+ ,601
+ ,573
+ ,569
+ ,562
+ ,468
+ ,483
+ ,460
+ ,411
+ ,458
+ ,455
+ ,600
+ ,605
+ ,545
+ ,549
+ ,415
+ ,568
+ ,577
+ ,517
+ ,558
+ ,518
+ ,489
+ ,502
+ ,569
+ ,540
+ ,550
+ ,557
+ ,542
+ ,542
+ ,582
+ ,525
+ ,584
+ ,562
+ ,639
+ ,613
+ ,604
+ ,613
+ ,625
+ ,654
+ ,638)
> 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] 454
> (np <- floor(n / par1))
[1] 37
> 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] 38 38 38 38 38 38 38 38 38 38 37 37
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 192.1 234.5 231.2 248.6 279.6 287.0 314.1 336.3 404.3 394.5 407.9 461.0
 [2,] 201.8 251.4 264.3 279.4 302.2 331.5 352.5 368.5 421.4 416.8 448.5 529.3
 [3,] 186.9 243.8 252.7 270.5 299.1 326.9 361.2 374.2 372.6 402.9 431.1 467.4
 [4,] 218.0 247.4 275.5 292.8 319.7 338.6 365.9 384.3 404.7 424.5 498.8 559.8
 [5,] 214.4 245.3 261.5 297.8 310.9 337.0 332.5 358.9 420.2 432.3 497.1 536.5
 [6,] 227.5 262.5 275.5 296.8 315.2 358.4 364.0 407.7 438.4 484.1 517.1 531.9
 [7,] 204.1 270.0 272.3 290.9 338.5 344.5 359.1 433.3 449.1 492.7 487.7 546.5
 [8,] 225.8 259.9 268.6 282.8 315.6 345.7 345.6 404.7 445.8 496.3 512.5 547.4
 [9,] 223.7 262.2 270.4 312.8 321.2 344.1 366.9 392.7 413.8 471.9 550.1 536.1
[10,] 244.7 244.9 267.7 303.2 318.5 317.4 370.2 409.7 420.5 491.2 532.5 482.8
[11,] 243.9 249.3 275.0 301.4 342.7 354.5 359.9 416.5 442.3 512.9 524.1 551.0
[12,] 257.3 268.2 272.6 289.8 261.4 345.2 366.6 414.3 438.9 482.4 515.7 532.9
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] 484.1 524.4 495.9 467.9 496.8 549.8 640.4 659.5 583.7 565.6 603.3 595.9
 [2,] 554.8 618.5 526.8 588.9 529.9 635.5 675.2 691.9 705.5 564.1 674.5 638.4
 [3,] 537.0 580.9 522.1 579.4 491.7 617.7 528.1 639.8 631.0 637.2 573.4 665.4
 [4,] 558.0 557.2 605.1 502.6 543.2 643.5 538.0 663.8 733.3 650.8 706.0 671.4
 [5,] 511.4 571.2 574.4 568.7 490.8 485.7 527.2 572.9 674.9 602.7 693.7 637.0
 [6,] 502.9 597.5 609.7 616.0 554.7 689.5 655.6 592.5 695.5 587.5 627.5 685.7
 [7,] 558.6 601.7 580.7 586.2 625.7 692.0 650.6 734.8 634.1 619.2 550.7 705.8
 [8,] 545.1 558.9 565.1 575.5 605.0 677.3 623.7 696.1 630.6 616.5 592.3 704.8
 [9,] 574.3 600.9 590.7 599.9 645.2 704.7 748.4 589.2 635.2 637.9 660.2 734.4
[10,] 542.2 601.0 571.5 568.2 645.2 668.6 727.4 662.9 554.1 557.9 597.3 674.2
[11,] 600.0 615.7 601.3 516.0 611.8 717.8 750.5 661.2 623.9 594.0 641.0 748.6
[12,] 588.6 578.1 567.3 493.4 600.3 689.8 678.9 672.1 679.3 668.7 663.6 763.4
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
 [1,] 658.0 634.3 496.7 399.2   474   451   504   473   475   508   451   483
 [2,] 627.5 690.7 660.9 429.8   539   537   545   531   589   560   538   460
 [3,] 528.9 676.3 601.2 395.6   468   494   547   462   495   558   555   411
 [4,] 488.3 675.4 655.2 472.0   543   549   585   479   484   516   532   458
 [5,] 575.5 595.6 591.6 486.0   532   544   562   533   536   549   530   455
 [6,] 735.6 712.4 606.1 525.0   565   598   520   552   555   595   553   600
 [7,] 685.3 735.8 560.7 396.0   535   583   581   547   565   586   626   605
 [8,] 613.6 544.4 368.3 511.0   534   582   590   562   564   597   601   545
 [9,] 629.5 567.0 371.6 525.0   546   589   562   524   573   592   573   549
[10,] 634.7 510.0 413.9 492.0   494   578   548   479   569   538   569   415
[11,] 652.6 564.0 413.9 517.0   552   561   567   445   588   590   562   568
[12,] 728.3 630.7 389.0 525.0   511   592   542   406   546   576   468   577
      [,37] [,38]
 [1,]   517   525
 [2,]   558   584
 [3,]   518   562
 [4,]   489   639
 [5,]   502   613
 [6,]   569   604
 [7,]   540   613
 [8,]   550   625
 [9,]   557   654
[10,]   542   638
[11,]   542    NA
[12,]   582    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]   9.7  16.9  33.1  30.8  22.6  44.5  38.4  32.2  17.1  22.3  40.6  68.3
 [2,] -14.9  -7.6 -11.6  -8.9  -3.1  -4.6   8.7   5.7 -48.8 -13.9 -17.4 -61.9
 [3,]  31.1   3.6  22.8  22.3  20.6  11.7   4.7  10.1  32.1  21.6  67.7  92.4
 [4,]  -3.6  -2.1 -14.0   5.0  -8.8  -1.6 -33.4 -25.4  15.5   7.8  -1.7 -23.3
 [5,]  13.1  17.2  14.0  -1.0   4.3  21.4  31.5  48.8  18.2  51.8  20.0  -4.6
 [6,] -23.4   7.5  -3.2  -5.9  23.3 -13.9  -4.9  25.6  10.7   8.6 -29.4  14.6
 [7,]  21.7 -10.1  -3.7  -8.1 -22.9   1.2 -13.5 -28.6  -3.3   3.6  24.8   0.9
 [8,]  -2.1   2.3   1.8  30.0   5.6  -1.6  21.3 -12.0 -32.0 -24.4  37.6 -11.3
 [9,]  21.0 -17.3  -2.7  -9.6  -2.7 -26.7   3.3  17.0   6.7  19.3 -17.6 -53.3
[10,]  -0.8   4.4   7.3  -1.8  24.2  37.1 -10.3   6.8  21.8  21.7  -8.4  68.2
[11,]  13.4  18.9  -2.4 -11.6 -81.3  -9.3   6.7  -2.2  -3.4 -30.5  -8.4 -18.1
[12,] -22.8 -37.0 -24.0 -10.2  25.6 -31.1 -30.3 -10.0 -44.4 -74.5 -54.7 -48.8
      [,13] [,14] [,15] [,16] [,17]  [,18]  [,19]  [,20]  [,21] [,22]  [,23]
 [1,]  70.7  94.1  30.9 121.0  33.1   85.7   34.8   32.4  121.8  -1.5   71.2
 [2,] -17.8 -37.6  -4.7  -9.5 -38.2  -17.8 -147.1  -52.1  -74.5  73.1 -101.1
 [3,]  21.0 -23.7  83.0 -76.8  51.5   25.8    9.9   24.0  102.3  13.6  132.6
 [4,] -46.6  14.0 -30.7  66.1 -52.4 -157.8  -10.8  -90.9  -58.4 -48.1  -12.3
 [5,]  -8.5  26.3  35.3  47.3  63.9  203.8  128.4   19.6   20.6 -15.2  -66.2
 [6,]  55.7   4.2 -29.0 -29.8  71.0    2.5   -5.0  142.3  -61.4  31.7  -76.8
 [7,] -13.5 -42.8 -15.6 -10.7 -20.7  -14.7  -26.9  -38.7   -3.5  -2.7   41.6
 [8,]  29.2  42.0  25.6  24.4  40.2   27.4  124.7 -106.9    4.6  21.4   67.9
 [9,] -32.1   0.1 -19.2 -31.7   0.0  -36.1  -21.0   73.7  -81.1 -80.0  -62.9
[10,]  57.8  14.7  29.8 -52.2 -33.4   49.2   23.1   -1.7   69.8  36.1   43.7
[11,] -11.4 -37.6 -34.0 -22.6 -11.5  -28.0  -71.6   10.9   55.4  74.7   22.6
[12,] -64.2 -82.2 -99.4   3.4 -50.5  -49.4  -19.4  -88.4 -113.7 -65.4  -67.7
       [,24] [,25]  [,26]  [,27]  [,28] [,29] [,30] [,31] [,32] [,33] [,34]
 [1,]   42.5 -30.5   56.4  164.2   30.6    65    86    41    58   114    52
 [2,]   27.0 -98.6  -14.4  -59.7  -34.2   -71   -43     2   -69   -94    -2
 [3,]    6.0 -40.6   -0.9   54.0   76.4    75    55    38    17   -11   -42
 [4,]  -34.4  87.2  -79.8  -63.6   14.0   -11    -5   -23    54    52    33
 [5,]   48.7 160.1  116.8   14.5   39.0    33    54   -42    19    19    46
 [6,]   20.1 -50.3   23.4  -45.4 -129.0   -30   -15    61    -5    10    -9
 [7,]   -1.0 -71.7 -191.4 -192.4  115.0    -1    -1     9    15    -1    11
 [8,]   29.6  15.9   22.6    3.3   14.0    12     7   -28   -38     9    -5
 [9,]  -60.2   5.2  -57.0   42.3  -33.0   -52   -11   -14   -45    -4   -54
[10,]   74.4  17.9   54.0    0.0   25.0    58   -17    19   -34    19    52
[11,]   14.8  75.7   66.7  -24.9    8.0   -41    31   -25   -39   -42   -14
[12,] -105.4 -94.0 -134.0   10.2  -51.0   -60   -88   -69    69   -38  -125
      [,35] [,36] [,37] [,38]
 [1,]    87   -23    41    59
 [2,]    17   -49   -40   -22
 [3,]   -23    47   -29    77
 [4,]    -2    -3    13   -26
 [5,]    23   145    67    -9
 [6,]    73     5   -29     9
 [7,]   -25   -60    10    12
 [8,]   -28     4     7    29
 [9,]    -4  -134   -15   -16
[10,]    -7   153     0    NA
[11,]   -94     9    40    NA
[12,]    15   -60   -57    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/1y5on1272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/28e581272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/38e581272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/41n4t1272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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] [,12]
[1,] 231.2 201.8 186.9 218.0 245.3 315.2 204.10 225.80 223.7 244.7 243.9 257.3
[2,] 399.2 421.4 395.6 424.5 432.3 484.1 449.10 445.80 413.8 415.0 442.3 414.3
[3,] 479.0 538.5 506.5 524.0 531.0 560.0 562.85 554.45 564.5 540.0 561.0 542.0
[4,] 525.0 618.5 573.4 605.1 572.9 606.1 619.20 601.00 600.9 578.0 601.3 600.3
[5,] 659.5 705.5 676.3 733.3 693.7 735.6 735.80 704.80 748.4 727.4 750.5 763.4

$n
 [1] 38 38 38 38 38 38 38 38 38 38 37 37

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 446.7562 487.9813 460.9281 477.7104 494.9628 528.7302 519.2517 514.6707
[2,] 511.2438 589.0187 552.0719 570.2896 567.0372 591.2698 606.4483 594.2293
         [,9]    [,10]    [,11]    [,12]
[1,] 516.5444 498.2215 519.6997 493.6864
[2,] 612.4556 581.7785 602.3003 590.3136

$out
[1] 192.1 214.4 227.5 262.5 275.5 296.8

$group
[1] 1 5 6 6 6 6

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

Warning message:
In bxp(list(stats = c(231.2, 399.2, 479, 525, 659.5, 201.8, 421.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/51n4t1272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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.0 -101.1 -42.00 -90.9 -15.2 -76.80 -60.0 -38.0 -81.10 -52.2 -81.3
[2,]  30.8  -52.1   4.70 -33.4  14.0 -29.00 -22.9  -2.1 -45.00  -0.8 -28.0
[3,]  41.0  -17.8  21.95  -9.8  22.2  -0.35  -3.6   8.0 -16.65  19.0  -9.3
[4,]  70.7   -4.7  54.00   7.8  48.8  20.10   3.6  27.4   0.00  43.7  13.4
[5,] 121.8   27.0 102.30  66.1  67.0  73.00  41.6  67.9  42.30  74.4  74.7
      [,12]
[1,] -134.0
[2,]  -74.5
[3,]  -51.0
[4,]  -24.0
[5,]   25.6

$n
 [1] 38 38 38 38 38 38 38 38 38 37 37 37

$conf
         [,1]       [,2]      [,3]        [,4]     [,5]      [,6]       [,7]
[1,] 30.77324 -29.949087  9.313925 -20.3599656 13.28042 -12.93481 -10.392211
[2,] 51.22676  -5.650913 34.586075   0.7599656 31.11958  12.23481   3.192211
           [,8]       [,9]     [,10]      [,11]     [,12]
[1,]  0.4388596 -28.183943  7.441107 -20.053667 -64.11740
[2,] 15.5611404  -5.116057 30.558893   1.453667 -37.88260

$out
 [1]  -30.5  164.2 -147.1   73.1  -76.8  132.6 -157.8   87.2  203.8  128.4
[11]  -66.2  160.1  116.8  -42.0  145.0  142.3 -129.0  -71.7 -191.4 -192.4
[21]  115.0  124.7 -106.9   73.7 -134.0  153.0   75.7  -94.0   69.0

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

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

Warning message:
In bxp(list(stats = c(-23, 30.8, 41, 70.7, 121.8, -101.1, -52.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/61n4t1272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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,] 186.90 234.50 252.7 270.50 279.60 317.40 332.50 336.30 372.60 394.50
[2,] 202.95 245.10 262.9 281.10 300.65 329.20 349.05 371.35 409.25 420.65
[3,] 220.85 250.35 269.5 291.85 315.40 341.35 360.55 398.70 420.95 477.15
[4,] 235.70 262.35 273.8 299.60 320.45 345.45 366.25 412.00 440.60 491.95
[5,] 257.30 270.00 275.5 312.80 342.70 358.40 370.20 433.30 449.10 512.90
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 407.90 461.00 484.10 524.40 495.90 467.90 490.80 549.80 527.20 572.90
[2,] 468.10 506.05 524.20 565.05 545.95 509.30 513.35 626.60 580.85 616.15
[3,] 505.65 534.50 549.95 589.20 572.95 572.10 577.50 672.95 653.10 662.05
[4,] 520.60 546.95 566.45 601.35 596.00 587.55 618.75 690.90 703.15 682.00
[5,] 550.10 559.80 600.00 618.50 609.70 616.00 645.20 717.80 750.50 734.80
      [,21]  [,22]  [,23]  [,24]  [,25] [,26]  [,27] [,28] [,29] [,30] [,31]
[1,] 554.10 557.90 550.70 595.90 488.30 510.0 368.30 395.6 468.0 494.0 504.0
[2,] 627.25 576.55 594.80 651.90 594.55 565.5 401.45 414.5 502.5 540.5 543.5
[3,] 634.65 609.60 634.25 679.95 632.10 632.5 528.70 489.0 534.5 569.5 555.0
[4,] 687.40 637.55 669.05 720.10 671.65 683.5 603.65 521.0 544.5 586.0 574.0
[5,] 733.30 668.70 706.00 763.40 735.60 735.8 660.90 525.0 565.0 598.0 590.0
     [,32] [,33] [,34] [,35] [,36] [,37] [,38]
[1,] 406.0 475.0 508.0   530 411.0 489.0   525
[2,] 467.5 515.5 543.5   531 456.5 517.5   584
[3,] 501.5 559.5 568.0   554 514.0 542.0   613
[4,] 540.0 571.0 591.0   571 572.5 557.5   638
[5,] 562.0 589.0 597.0   626 605.0 582.0   654

$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 12 12 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 205.9125 242.4822 264.5284 283.412 306.3691 333.9383 352.7050 380.1593
[2,] 235.7875 258.2178 274.4716 300.288 324.4309 348.7617 368.3950 417.2407
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 406.6511 444.6296 481.7044 515.8452 530.6795 572.6433 550.1219 536.4096
[2,] 435.2489 509.6704 529.5956 553.1548 569.2205 605.7567 595.7781 607.7904
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 529.4264 643.6223 597.3181 632.0154 607.2152 581.7775 600.3841 648.8435
[2,] 625.5736 702.2777 708.8819 692.0846 662.0848 637.4225 668.1159 711.0565
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30]    [,31]    [,32]
[1,] 596.9342 578.6794 436.4752 440.4246 515.3435 548.7471 541.0887 468.4323
[2,] 667.2658 686.3206 620.9248 537.5754 553.6565 590.2529 568.9113 534.5677
        [,33]    [,34]    [,35]    [,36]    [,37]    [,38]
[1,] 534.1861 546.3349 535.7557 461.0916 523.7557 586.0194
[2,] 584.8139 589.6651 572.2443 566.9084 560.2443 639.9806

$out
[1] 231.2 248.6 261.4 287.0 314.1 485.7 451.0 451.0 468.0

$group
[1]  3  4  5  6  7 18 30 35 35

$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" "32" "33" "34" "35" "36" "37" NA  

Warning message:
In bxp(list(stats = c(186.9, 202.95, 220.85, 235.7, 257.3, 234.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7uf4d1272311882.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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,] 480.6579 479.0 483.9875
[2,] 499.3961 527.5 500.8000
[3,] 513.8211 541.0 515.4813
[4,] 525.0069 560.5 525.5062
[5,] 531.2526 564.5 547.1875

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 502.1398 525.9485 504.2126
[2,] 525.5023 556.0515 526.7499

$out
[1] 460.7263 462.6625

$group
[1] 1 3

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

Warning message:
In bxp(list(stats = c(480.657894736842, 499.396052631579, 513.821052631579,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1y5on1272311882.ps tmp/1y5on1272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/28e581272311882.ps tmp/28e581272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/38e581272311882.ps tmp/38e581272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/41n4t1272311882.ps tmp/41n4t1272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/51n4t1272311882.ps tmp/51n4t1272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/61n4t1272311882.ps tmp/61n4t1272311882.png",intern=TRUE))
character(0)
> try(system("convert tmp/7uf4d1272311882.ps tmp/7uf4d1272311882.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.347   0.999   1.794