R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(553.12
+ ,568.15
+ ,552.75
+ ,510.65
+ ,524.93
+ ,532.95
+ ,540.84
+ ,540.22
+ ,553.75
+ ,529.69
+ ,525.93
+ ,527.31
+ ,527.31
+ ,512.03
+ ,502.76
+ ,496.62
+ ,492.23
+ ,495.12
+ ,469.93
+ ,492.36
+ ,497.87
+ ,480.21
+ ,462.29
+ ,456.03
+ ,456.22
+ ,460.41
+ ,466.59
+ ,441.37
+ ,455.31
+ ,426.96
+ ,419.7
+ ,419.7
+ ,416.44
+ ,404.04
+ ,388.63
+ ,397.65
+ ,390.38
+ ,378.1
+ ,384.87
+ ,419.19
+ ,427.96
+ ,413.81
+ ,408.04
+ ,410.3
+ ,405.66
+ ,400.9
+ ,387
+ ,388.25
+ ,390
+ ,416.44
+ ,436.11
+ ,428.79
+ ,424.16
+ ,409.12
+ ,392.01
+ ,388.37
+ ,373.97
+ ,358.93
+ ,371.96
+ ,353.92
+ ,364.95
+ ,340.52
+ ,353.67
+ ,361.19
+ ,364.7
+ ,359.43
+ ,371.21
+ ,385.24
+ ,389.63
+ ,433.23
+ ,407.79
+ ,400.9
+ ,385.87
+ ,406.67
+ ,406.04
+ ,418.19
+ ,429.22
+ ,420.95
+ ,402.78
+ ,391.01
+ ,416.94
+ ,397.14
+ ,406.67
+ ,419.44
+ ,422.2
+ ,435.98
+ ,470.81
+ ,504.51
+ ,497.25
+ ,508.52
+ ,522.43
+ ,551.62
+ ,537.59
+ ,559.76
+ ,558.26
+ ,563.27
+ ,558.01
+ ,563.27
+ ,564.15
+ ,582.81
+ ,592.96
+ ,602.73
+ ,581.06
+ ,595.47
+ ,605.36
+ ,615.39
+ ,602.48
+ ,565.65
+ ,565.65
+ ,566.9
+ ,558.63
+ ,547.48
+ ,543.72
+ ,517.42
+ ,526.31
+ ,512.4
+ ,496.12
+ ,503.63
+ ,501.13
+ ,499.88
+ ,501.13
+ ,494.86
+ ,488.6
+ ,482.34
+ ,447.26
+ ,440.99
+ ,418.44
+ ,418.44
+ ,412.18
+ ,394.64
+ ,334.5
+ ,328.24
+ ,319.47
+ ,323.23
+ ,328.24
+ ,365.82
+ ,371.88
+ ,340.6
+ ,337.48
+ ,337.22
+ ,313.96
+ ,308.39
+ ,307.57
+ ,295.4
+ ,297.04
+ ,341.89
+ ,341.18
+ ,352.12
+ ,367.36
+ ,364.88
+ ,363.09
+ ,351.25
+ ,349.29
+ ,335.47
+ ,320.1
+ ,310.7
+ ,312.39
+ ,309.68
+ ,309.67
+ ,328.92
+ ,337.01
+ ,327.79
+ ,324.38
+ ,313.93
+ ,310.83
+ ,316.62
+ ,325.5
+ ,320.03
+ ,320.1
+ ,338.13
+ ,379.25
+ ,376.82
+ ,398.37
+ ,394.46
+ ,411.95
+ ,425.91
+ ,444.48
+ ,433.45
+ ,446.07
+ ,426.04
+ ,447.06
+ ,467.94
+ ,516.73
+ ,520.27
+ ,516.49
+ ,519.41
+ ,537.66
+ ,547.44
+ ,507.04
+ ,495.99
+ ,436.58
+ ,453.11
+ ,456.77
+ ,450.38
+ ,439.18
+ ,416.56
+ ,440.06
+ ,447.61
+ ,420.32
+ ,417.67
+ ,404.38
+ ,416.41
+ ,419.48
+ ,417.72
+ ,408.62
+ ,442.94
+ ,425.82
+ ,451.19
+ ,467.49
+ ,478.76
+ ,478.56
+ ,427.63
+ ,448.81
+ ,435.41
+ ,434.67
+ ,413.62
+ ,399.02
+ ,406.64
+ ,384.83
+ ,379.81
+ ,355.7
+ ,348.24
+ ,308.83
+ ,296.93
+ ,280.11
+ ,286.85
+ ,294.93
+ ,294.77
+ ,299.37
+ ,287.1
+ ,297.46
+ ,298.88
+ ,288.74
+ ,288.32
+ ,286.32
+ ,254.34
+ ,247.09
+ ,247.29
+ ,255.49
+ ,267.26
+ ,276.44
+ ,260.07
+ ,267.1
+ ,273.81
+ ,290.37
+ ,293.98
+ ,302.36
+ ,289.92
+ ,283.27
+ ,279.87
+ ,267.66
+ ,286.88
+ ,309.69
+ ,323.95
+ ,315.36
+ ,327.52
+ ,325.69
+ ,326.92
+ ,328.26
+ ,348.94
+ ,340.53
+ ,330.29
+ ,335.91
+ ,376.13
+ ,444.11
+ ,516.35
+ ,529.16
+ ,525.07
+ ,519.78
+ ,548.84
+ ,539.68
+ ,534.99
+ ,584.34
+ ,664.34
+ ,691.81
+ ,689.34
+ ,725.81
+ ,734.89
+ ,681.58
+ ,685.72
+ ,633.01
+ ,680.28
+ ,684.95
+ ,653.47
+ ,647.28
+ ,602.73
+ ,589.76
+ ,588.41
+ ,613.51
+ ,611.93
+ ,587.69
+ ,554.63
+ ,533.09
+ ,560.59
+ ,553.05
+ ,528.97
+ ,500.93
+ ,508.86
+ ,537.86
+ ,547.3
+ ,556.94
+ ,549.33
+ ,545.18
+ ,543.69
+ ,543.49
+ ,553.32
+ ,563.57
+ ,531.87
+ ,517.85
+ ,500.83
+ ,481.51
+ ,479.73
+ ,496.6
+ ,520.76
+ ,528.71
+ ,515.46
+ ,522
+ ,515.63
+ ,522.87
+ ,534.29
+ ,521.83
+ ,524.66
+ ,640.02
+ ,644.21
+ ,715.51
+ ,706.96
+ ,724.75
+ ,742.55
+ ,788.25
+ ,787.58
+ ,761.64
+ ,732.83
+ ,765.24
+ ,807.45
+ ,828.19
+ ,817.45
+ ,840.86
+ ,843.77
+ ,835.3
+ ,823.97
+ ,781.49
+ ,778.46
+ ,793.63
+ ,717.57
+ ,656.35
+ ,510.28
+ ,450.98
+ ,477.83
+ ,464.59
+ ,449.72
+ ,460.47
+ ,516.6
+ ,599.31
+ ,609.15
+ ,609.3
+ ,655.49
+ ,690.91
+ ,743.59
+ ,756.8
+ ,780.4
+ ,842.74
+ ,818.72
+ ,847.84
+ ,817.97
+ ,763.94
+ ,762.81
+ ,777.29
+ ,786.94
+ ,766.28)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 369
> (np <- floor(n / par1))
[1] 30
> 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 30 30 30
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 553.12 527.31 456.22 390.38 390.00 364.95 385.87 422.20 558.01 565.65
 [2,] 568.15 512.03 460.41 378.10 416.44 340.52 406.67 435.98 563.27 566.90
 [3,] 552.75 502.76 466.59 384.87 436.11 353.67 406.04 470.81 564.15 558.63
 [4,] 510.65 496.62 441.37 419.19 428.79 361.19 418.19 504.51 582.81 547.48
 [5,] 524.93 492.23 455.31 427.96 424.16 364.70 429.22 497.25 592.96 543.72
 [6,] 532.95 495.12 426.96 413.81 409.12 359.43 420.95 508.52 602.73 517.42
 [7,] 540.84 469.93 419.70 408.04 392.01 371.21 402.78 522.43 581.06 526.31
 [8,] 540.22 492.36 419.70 410.30 388.37 385.24 391.01 551.62 595.47 512.40
 [9,] 553.75 497.87 416.44 405.66 373.97 389.63 416.94 537.59 605.36 496.12
[10,] 529.69 480.21 404.04 400.90 358.93 433.23 397.14 559.76 615.39 503.63
[11,] 525.93 462.29 388.63 387.00 371.96 407.79 406.67 558.26 602.48 501.13
[12,] 527.31 456.03 397.65 388.25 353.92 400.90 419.44 563.27 565.65 499.88
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 501.13 319.47 297.04 312.39 320.10 447.06 456.77 408.62 399.02 299.37
 [2,] 494.86 323.23 341.89 309.68 338.13 467.94 450.38 442.94 406.64 287.10
 [3,] 488.60 328.24 341.18 309.67 379.25 516.73 439.18 425.82 384.83 297.46
 [4,] 482.34 365.82 352.12 328.92 376.82 520.27 416.56 451.19 379.81 298.88
 [5,] 447.26 371.88 367.36 337.01 398.37 516.49 440.06 467.49 355.70 288.74
 [6,] 440.99 340.60 364.88 327.79 394.46 519.41 447.61 478.76 348.24 288.32
 [7,] 418.44 337.48 363.09 324.38 411.95 537.66 420.32 478.56 308.83 286.32
 [8,] 418.44 337.22 351.25 313.93 425.91 547.44 417.67 427.63 296.93 254.34
 [9,] 412.18 313.96 349.29 310.83 444.48 507.04 404.38 448.81 280.11 247.09
[10,] 394.64 308.39 335.47 316.62 433.45 495.99 416.41 435.41 286.85 247.29
[11,] 334.50 307.57 320.10 325.50 446.07 436.58 419.48 434.67 294.93 255.49
[12,] 328.24 295.40 310.70 320.03 426.04 453.11 417.72 413.62 294.77 267.26
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30]
 [1,] 276.44 309.69 444.11 725.81 613.51 556.94 496.60 715.51 840.86 464.59
 [2,] 260.07 323.95 516.35 734.89 611.93 549.33 520.76 706.96 843.77 449.72
 [3,] 267.10 315.36 529.16 681.58 587.69 545.18 528.71 724.75 835.30 460.47
 [4,] 273.81 327.52 525.07 685.72 554.63 543.69 515.46 742.55 823.97 516.60
 [5,] 290.37 325.69 519.78 633.01 533.09 543.49 522.00 788.25 781.49 599.31
 [6,] 293.98 326.92 548.84 680.28 560.59 553.32 515.63 787.58 778.46 609.15
 [7,] 302.36 328.26 539.68 684.95 553.05 563.57 522.87 761.64 793.63 609.30
 [8,] 289.92 348.94 534.99 653.47 528.97 531.87 534.29 732.83 717.57 655.49
 [9,] 283.27 340.53 584.34 647.28 500.93 517.85 521.83 765.24 656.35 690.91
[10,] 279.87 330.29 664.34 602.73 508.86 500.83 524.66 807.45 510.28 743.59
[11,] 267.66 335.91 691.81 589.76 537.86 481.51 640.02 828.19 450.98 756.80
[12,] 286.88 376.13 689.34 588.41 547.30 479.73 644.21 817.45 477.83 780.40
       [,31]
 [1,] 842.74
 [2,] 818.72
 [3,] 847.84
 [4,] 817.97
 [5,] 763.94
 [6,] 762.81
 [7,] 777.29
 [8,] 786.94
 [9,] 766.28
[10,]     NA
[11,]     NA
[12,]     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,]  15.03 -15.28   4.19 -12.28  26.44 -24.43  20.80  13.78   5.26   1.25
 [2,] -15.40  -9.27   6.18   6.77  19.67  13.15  -0.63  34.83   0.88  -8.27
 [3,] -42.10  -6.14 -25.22  34.32  -7.32   7.52  12.15  33.70  18.66 -11.15
 [4,]  14.28  -4.39  13.94   8.77  -4.63   3.51  11.03  -7.26  10.15  -3.76
 [5,]   8.02   2.89 -28.35 -14.15 -15.04  -5.27  -8.27  11.27   9.77 -26.30
 [6,]   7.89 -25.19  -7.26  -5.77 -17.11  11.78 -18.17  13.91 -21.67   8.89
 [7,]  -0.62  22.43   0.00   2.26  -3.64  14.03 -11.77  29.19  14.41 -13.91
 [8,]  13.53   5.51  -3.26  -4.64 -14.40   4.39  25.93 -14.03   9.89 -16.28
 [9,] -24.06 -17.66 -12.40  -4.76 -15.04  43.60 -19.80  22.17  10.03   7.51
[10,]  -3.76 -17.92 -15.41 -13.90  13.03 -25.44   9.53  -1.50 -12.91  -2.50
[11,]   1.38  -6.26   9.02   1.25 -18.04  -6.89  12.77   5.01 -36.83  -1.25
[12,]   0.00   0.19  -7.27   1.75  11.03 -15.03   2.76  -5.26   0.00   1.25
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,]  -6.27   3.76  44.85  -2.71  18.03  20.88  -6.39  34.32   7.62 -12.27
 [2,]  -6.26   5.01  -0.71  -0.01  41.12  48.79 -11.20 -17.12 -21.81  10.36
 [3,]  -6.26  37.58  10.94  19.25  -2.43   3.54 -22.62  25.37  -5.02   1.42
 [4,] -35.08   6.06  15.24   8.09  21.55  -3.78  23.50  16.30 -24.11 -10.14
 [5,]  -6.27 -31.28  -2.48  -9.22  -3.91   2.92   7.55  11.27  -7.46  -0.42
 [6,] -22.55  -3.12  -1.79  -3.41  17.49  18.25 -27.29  -0.20 -39.41  -2.00
 [7,]   0.00  -0.26 -11.84 -10.45  13.96   9.78  -2.65 -50.93 -11.90 -31.98
 [8,]  -6.26 -23.26  -1.96  -3.10  18.57 -40.40 -13.29  21.18 -16.82  -7.25
 [9,] -17.54  -5.57 -13.82   5.79 -11.03 -11.05  12.03 -13.40   6.74   0.20
[10,] -60.14  -0.82 -15.37   8.88  12.62 -59.41   3.07  -0.74   8.08   8.20
[11,]  -6.26 -12.17  -9.40  -5.47 -20.03  16.53  -1.76 -21.05  -0.16  11.77
[12,]  -8.77   1.64   1.69   0.07  21.02   3.66  -9.10 -14.60   4.60   9.18
       [,21]  [,22] [,23]  [,24]  [,25]  [,26]  [,27]  [,28]   [,29]  [,30]
 [1,] -16.37  14.26 72.24   9.08  -1.58  -7.61  24.16  -8.55    2.91 -14.87
 [2,]   7.03  -8.59 12.81 -53.31 -24.24  -4.15   7.95  17.79   -8.47  10.75
 [3,]   6.71  12.16 -4.09   4.14 -33.06  -1.49 -13.25  17.80  -11.33  56.13
 [4,]  16.56  -1.83 -5.29 -52.71 -21.54  -0.20   6.54  45.70  -42.48  82.71
 [5,]   3.61   1.23 29.06  47.27  27.50   9.83  -6.37  -0.67   -3.03   9.84
 [6,]   8.38   1.34 -9.16   4.67  -7.54  10.25   7.24 -25.94   15.17   0.15
 [7,] -12.44  20.68 -4.69 -31.48 -24.08 -31.70  11.42 -28.81  -76.06  46.19
 [8,]  -6.65  -8.41 49.35  -6.19 -28.04 -14.02 -12.46  32.41  -61.22  35.42
 [9,]  -3.40 -10.24 80.00 -44.55   7.93 -17.02   2.83  42.21 -146.07  52.68
[10,] -12.21   5.62 27.47 -12.97  29.00 -19.32 115.36  20.74  -59.30  13.21
[11,]  19.22  40.22 -2.47  -1.35   9.44  -1.78   4.19 -10.74   26.85  23.60
[12,]  22.81  67.98 36.47  25.10   9.64  16.87  71.30  23.41  -13.24  62.34
       [,31]
 [1,] -24.02
 [2,]  29.12
 [3,] -29.87
 [4,] -54.03
 [5,]  -1.13
 [6,]  14.48
 [7,]   9.65
 [8,] -20.66
 [9,]     NA
[10,]     NA
[11,]     NA
[12,]     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/rcomp/tmp/1ner11321397285.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/www/rcomp/tmp/2alij1321397285.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/www/rcomp/tmp/3rdxp1321397285.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/www/rcomp/tmp/4x1na1321397285.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]
[1,] 276.44 260.070 267.100 273.810 288.740 288.32 286.320 254.340 247.09
[2,] 375.41 359.995 382.040 378.315 385.125 379.67 381.610 386.805 381.80
[3,] 447.06 450.380 466.590 482.340 467.490 478.76 469.930 427.630 448.81
[4,] 555.03 556.300 548.965 534.380 538.290 551.08 546.945 543.830 545.67
[5,] 725.81 843.770 724.750 742.550 763.940 787.58 793.630 732.830 766.28
      [,10]   [,11]  [,12]
[1,] 247.29 255.490 267.26
[2,] 358.93 335.910 353.92
[3,] 434.43 435.625 422.74
[4,] 524.66 537.860 547.30
[5,] 743.59 828.190 817.45

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]   [,8]
[1,] 396.0881 394.6733 419.2206 438.0524 424.0254 430.1179 423.0118 383.07
[2,] 498.0319 506.0867 513.9594 526.6276 510.9546 527.4021 516.8482 472.19
         [,9]    [,10]    [,11]    [,12]
[1,] 402.3076 386.6223 377.3690 366.9562
[2,] 495.3124 482.2377 493.8810 478.5238

$out
 [1] 840.86 842.74 835.30 847.84 823.97 817.97 788.25 781.49 786.94 807.45

$group
 [1]  1  1  3  3  4  4  5  5  8 10

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

Warning message:
In bxp(list(stats = c(276.44, 375.41, 447.06, 555.03, 725.81, 260.07,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5pivm1321397285.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]
[1,] -24.43 -24.24 -42.100 -35.080 -28.350 -39.410 -31.980 -40.400 -44.550
[2,]  -8.08  -8.53  -9.235  -6.275  -6.915 -13.135 -13.175 -14.215 -15.040
[3,]   3.76   0.88   1.420   3.510  -0.670  -1.790  -2.650  -6.260  -5.165
[4,]  16.53  11.78  14.980  14.110   8.895   8.635  10.600   7.700   7.930
[5,]  44.85  41.12  37.580  23.500  29.060  18.250  46.190  35.420  42.210
      [,10]  [,11]   [,12]
[1,] -25.44 -21.05 -15.030
[2,] -15.37  -6.89   0.000
[3,]  -1.16  -1.30   2.255
[4,]   9.53   9.44  21.020
[5,]  29.00  26.85  36.470

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -3.223737 -4.883498 -5.451645 -2.274781 -5.156505 -7.967812 -9.396783
[2,] 10.743737  6.643498  8.291645  9.294781  3.816505  4.387812  4.096783
             [,8]       [,9]     [,10]     [,11]     [,12]
[1,] -12.47895931 -11.791092 -8.342834 -6.010670 -3.808581
[2,]  -0.04104069   1.461092  6.022834  3.410670  8.318581

$out
 [1]   72.24   48.79  -53.31   56.13  -52.71   45.70  -42.48   82.71  -54.03
[10]  -31.28   47.27  -50.93  -76.06   49.35  -61.22   43.60   80.00 -146.07
[19]   52.68  -60.14  -59.41  115.36  -59.30  -36.83   40.22   67.98   71.30
[28]   62.34

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

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

Warning message:
In bxp(list(stats = c(-24.43, -8.07999999999998, 3.75999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6omnv1321397285.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]
[1,] 510.650 456.030 388.630 378.100 353.920 340.520 385.870 422.200 558.010
[2,] 526.620 475.070 410.240 387.625 372.965 360.310 399.960 484.030 564.900
[3,] 536.585 493.740 423.330 403.280 391.005 368.080 406.670 515.475 587.885
[4,] 552.935 500.315 455.765 412.055 420.300 395.265 418.815 554.940 602.605
[5,] 568.150 527.310 466.590 427.960 436.110 433.230 429.220 563.270 615.390
       [,10]   [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18]
[1,] 496.120 328.240 295.400 297.040 309.670 320.100 436.580 404.380 408.620
[2,] 502.380 403.410 311.175 327.785 311.610 378.035 460.525 417.115 426.725
[3,] 521.865 429.715 325.735 345.590 318.325 405.160 511.765 419.900 439.175
[4,] 553.055 485.470 339.040 357.605 326.645 429.745 519.840 443.835 459.340
[5,] 566.900 501.130 371.880 367.360 337.010 446.070 547.440 456.770 478.760
       [,19]   [,20]   [,21]  [,22]   [,23]   [,24]   [,25]   [,26]   [,27]
[1,] 280.110 247.090 260.070 309.69 444.110 588.410 500.930 479.730 515.460
[2,] 294.850 254.915 270.735 324.82 522.425 617.870 531.030 509.340 518.195
[3,] 328.535 286.710 281.570 327.89 537.335 666.875 550.175 543.590 522.435
[4,] 382.320 293.100 290.145 338.22 624.340 685.335 574.140 551.325 531.500
[5,] 406.640 299.370 302.360 348.94 691.810 734.890 613.510 563.570 534.290
      [,28]   [,29]   [,30]  [,31]
[1,] 706.96 450.980 449.720 762.81
[2,] 728.79 583.315 490.595 766.28
[3,] 763.44 779.975 609.225 786.94
[4,] 797.85 829.635 717.250 818.72
[5,] 828.19 843.770 780.400 847.84

$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  9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 524.5826 482.2256 402.5657 392.1373 369.4152 352.1368 398.0701 483.1325
[2,] 548.5874 505.2544 444.0943 414.4227 412.5948 384.0232 415.2699 547.8175
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]   [,16]
[1,] 570.6875 498.7518 392.2869 313.0256 331.9889 311.4674 381.5747 484.711
[2,] 605.0825 544.9782 467.1431 338.4444 359.1911 325.1826 428.7453 538.819
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 407.7128 424.2991 288.6393 269.2936 272.7170 321.7782 490.8509 636.1038
[2,] 432.0872 454.0509 368.4307 304.1264 290.4230 334.0018 583.8191 697.6462
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30]    [,31]
[1,] 530.5122 524.4404 516.3665 731.9413 667.6268 505.8461 759.3216
[2,] 569.8378 562.7396 528.5035 794.9387 892.3232 712.6039 814.5584

$out
[1] 376.13 496.60 640.02 644.21

$group
[1] 22 27 27 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] NA  

Warning message:
In bxp(list(stats = c(510.65, 526.62, 536.585, 552.935, 568.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/77r2q1321397285.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,] 458.9177 422.7400 439.9000
[2,] 466.9145 435.0275 453.6175
[3,] 478.0716 449.5950 462.7212
[4,] 483.3623 468.7100 465.2687
[5,] 485.6655 482.3400 465.5025

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 470.5697 434.2322 457.4070
[2,] 485.5736 464.9578 468.0355

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(458.917666666667, 466.914483870968, 478.071612903226,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ner11321397285.ps tmp/1ner11321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/2alij1321397285.ps tmp/2alij1321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/3rdxp1321397285.ps tmp/3rdxp1321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/4x1na1321397285.ps tmp/4x1na1321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pivm1321397285.ps tmp/5pivm1321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/6omnv1321397285.ps tmp/6omnv1321397285.png",intern=TRUE))
character(0)
> try(system("convert tmp/77r2q1321397285.ps tmp/77r2q1321397285.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.650   0.250   2.888