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(414.89
+ ,444.50
+ ,481.29
+ ,491.09
+ ,419.70
+ ,432.88
+ ,438.01
+ ,412.84
+ ,402.91
+ ,416.20
+ ,411.80
+ ,393.60
+ ,381.70
+ ,388.34
+ ,370.89
+ ,385.96
+ ,394.26
+ ,381.37
+ ,377.40
+ ,377.70
+ ,347.71
+ ,347.70
+ ,340.20
+ ,341.00
+ ,342.00
+ ,319.54
+ ,302.79
+ ,299.10
+ ,313.50
+ ,326.80
+ ,316.00
+ ,316.50
+ ,317.20
+ ,330.40
+ ,323.35
+ ,325.85
+ ,321.50
+ ,321.90
+ ,347.48
+ ,338.89
+ ,345.70
+ ,340.44
+ ,342.40
+ ,342.70
+ ,348.34
+ ,376.66
+ ,417.73
+ ,423.51
+ ,397.56
+ ,390.92
+ ,408.26
+ ,401.12
+ ,408.91
+ ,438.35
+ ,460.23
+ ,449.59
+ ,450.52
+ ,461.15
+ ,460.38
+ ,465.35
+ ,467.57
+ ,486.24
+ ,476.58
+ ,442.07
+ ,443.61
+ ,451.55
+ ,451.07
+ ,451.33
+ ,437.63
+ ,431.28
+ ,413.46
+ ,406.78
+ ,420.17
+ ,419.05
+ ,404.01
+ ,387.51
+ ,390.15
+ ,384.40
+ ,371.05
+ ,367.60
+ ,375.04
+ ,365.14
+ ,361.75
+ ,366.88
+ ,394.26
+ ,409.39
+ ,410.11
+ ,416.81
+ ,393.06
+ ,374.24
+ ,369.05
+ ,352.33
+ ,362.53
+ ,394.73
+ ,389.32
+ ,380.74
+ ,381.73
+ ,376.95
+ ,383.64
+ ,363.83
+ ,363.34
+ ,358.38
+ ,356.95
+ ,366.72
+ ,367.69
+ ,356.31
+ ,348.74
+ ,358.69
+ ,360.17
+ ,361.73
+ ,354.45
+ ,353.91
+ ,344.34
+ ,338.62
+ ,337.24
+ ,340.81
+ ,352.72
+ ,343.06
+ ,345.43
+ ,344.38
+ ,335.02
+ ,334.82
+ ,329.01
+ ,329.31
+ ,330.08
+ ,342.15
+ ,367.18
+ ,371.89
+ ,392.19
+ ,378.84
+ ,355.28
+ ,364.18
+ ,373.83
+ ,383.30
+ ,386.88
+ ,381.91
+ ,384.13
+ ,377.27
+ ,381.43
+ ,385.64
+ ,385.49
+ ,380.36
+ ,391.58
+ ,389.77
+ ,384.39
+ ,379.29
+ ,378.55
+ ,376.64
+ ,382.12
+ ,391.03
+ ,385.22
+ ,387.56
+ ,386.23
+ ,383.67
+ ,383.06
+ ,383.14
+ ,385.31
+ ,387.44
+ ,399.45
+ ,404.76
+ ,396.21
+ ,392.85
+ ,391.93
+ ,385.27
+ ,383.47
+ ,387.35
+ ,383.14
+ ,381.07
+ ,377.85
+ ,369.00
+ ,355.11
+ ,346.58
+ ,351.81
+ ,344.47
+ ,343.84
+ ,340.76
+ ,324.10
+ ,324.01
+ ,322.82
+ ,324.87
+ ,306.04
+ ,288.74
+ ,289.10
+ ,297.49
+ ,295.94
+ ,308.29
+ ,299.10
+ ,292.32
+ ,292.87
+ ,284.11
+ ,288.98
+ ,295.93
+ ,294.12
+ ,291.68
+ ,287.08
+ ,287.33
+ ,285.96
+ ,282.62
+ ,276.44
+ ,261.31
+ ,256.08
+ ,256.69
+ ,264.74
+ ,310.72
+ ,293.18
+ ,283.07
+ ,284.32
+ ,299.86
+ ,286.39
+ ,279.69
+ ,275.19
+ ,285.73
+ ,281.59
+ ,274.47
+ ,273.68
+ ,270.00
+ ,266.01
+ ,271.45
+ ,265.49
+ ,261.87
+ ,263.03
+ ,260.48
+ ,272.36
+ ,270.23
+ ,267.53
+ ,272.39
+ ,283.42
+ ,283.06
+ ,276.16
+ ,275.85
+ ,281.51
+ ,295.50
+ ,294.06
+ ,302.68
+ ,314.49
+ ,321.18
+ ,313.29
+ ,310.26
+ ,319.14
+ ,316.56
+ ,319.07
+ ,331.92
+ ,356.86
+ ,358.97
+ ,340.55
+ ,328.18
+ ,355.68
+ ,356.35
+ ,351.02
+ ,359.77
+ ,378.95
+ ,378.92
+ ,389.91
+ ,406.95
+ ,413.79
+ ,404.88
+ ,406.67
+ ,403.26
+ ,383.78
+ ,392.37
+ ,398.09
+ ,400.51
+ ,405.28
+ ,420.46
+ ,439.38
+ ,442.08
+ ,424.03
+ ,423.35
+ ,433.85
+ ,429.23
+ ,421.87
+ ,430.66
+ ,424.48
+ ,437.93
+ ,456.05
+ ,469.90
+ ,476.67
+ ,510.10
+ ,549.86
+ ,555.00
+ ,557.09
+ ,610.65
+ ,675.39
+ ,596.15
+ ,633.71
+ ,632.59
+ ,598.19
+ ,585.78
+ ,627.83
+ ,629.79
+ ,631.17
+ ,664.75
+ ,654.90
+ ,679.37
+ ,667.31
+ ,655.66
+ ,665.38
+ ,665.41
+ ,712.65
+ ,754.60
+ ,806.25
+ ,803.20
+ ,889.60
+ ,922.30
+ ,968.43
+ ,909.71
+ ,888.66
+ ,889.49
+ ,939.77
+ ,839.03
+ ,829.93
+ ,806.62
+ ,760.86
+ ,816.09
+ ,858.69
+ ,943.00
+ ,924.27
+ ,890.20
+ ,928.65
+ ,945.67
+ ,934.23
+ ,949.38
+ ,996.59
+ ,1043.16
+ ,1127.04
+ ,1134.72
+ ,1117.96
+ ,1095.41
+ ,1113.34
+ ,1148.69
+ ,1205.43
+ ,1232.92
+ ,1192.97
+ ,1215.81
+ ,1270.98
+ ,1342.02
+ ,1369.89
+ ,1390.55
+ ,1356.40
+ ,1372.73
+ ,1424.00
+ ,1479.76
+ ,1512.60
+ ,1528.66
+ ,1572.21
+ ,1757.21
+ ,1770.95
+ ,1665.21
+ ,1738.11
+ ,1641.84
+ ,1652.21
+ ,1742.14
+ ,1673.77
+ ,1649.69
+ ,1591.19
+ ,1598.76
+ ,1589.90
+ ,1630.31
+ ,1744.81
+ ,1746.58
+ ,1721.64)
> 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] 361
> (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 30 30 30 30 30 30 30 30 30 30 30
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 414.89 381.70 342.00 321.50 397.56 467.57 420.17 394.26 381.73 360.17
 [2,] 444.50 388.34 319.54 321.90 390.92 486.24 419.05 409.39 376.95 361.73
 [3,] 481.29 370.89 302.79 347.48 408.26 476.58 404.01 410.11 383.64 354.45
 [4,] 491.09 385.96 299.10 338.89 401.12 442.07 387.51 416.81 363.83 353.91
 [5,] 419.70 394.26 313.50 345.70 408.91 443.61 390.15 393.06 363.34 344.34
 [6,] 432.88 381.37 326.80 340.44 438.35 451.55 384.40 374.24 358.38 338.62
 [7,] 438.01 377.40 316.00 342.40 460.23 451.07 371.05 369.05 356.95 337.24
 [8,] 412.84 377.70 316.50 342.70 449.59 451.33 367.60 352.33 366.72 340.81
 [9,] 402.91 347.71 317.20 348.34 450.52 437.63 375.04 362.53 367.69 352.72
[10,] 416.20 347.70 330.40 376.66 461.15 431.28 365.14 394.73 356.31 343.06
[11,] 411.80 340.20 323.35 417.73 460.38 413.46 361.75 389.32 348.74 345.43
[12,] 393.60 341.00 325.85 423.51 465.35 406.78 366.88 380.74 358.69 344.38
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 335.02 373.83 384.39 385.31 377.85 306.04 294.12 293.18 266.01 276.16
 [2,] 334.82 383.30 379.29 387.44 369.00 288.74 291.68 283.07 271.45 275.85
 [3,] 329.01 386.88 378.55 399.45 355.11 289.10 287.08 284.32 265.49 281.51
 [4,] 329.31 381.91 376.64 404.76 346.58 297.49 287.33 299.86 261.87 295.50
 [5,] 330.08 384.13 382.12 396.21 351.81 295.94 285.96 286.39 263.03 294.06
 [6,] 342.15 377.27 391.03 392.85 344.47 308.29 282.62 279.69 260.48 302.68
 [7,] 367.18 381.43 385.22 391.93 343.84 299.10 276.44 275.19 272.36 314.49
 [8,] 371.89 385.64 387.56 385.27 340.76 292.32 261.31 285.73 270.23 321.18
 [9,] 392.19 385.49 386.23 383.47 324.10 292.87 256.08 281.59 267.53 313.29
[10,] 378.84 380.36 383.67 387.35 324.01 284.11 256.69 274.47 272.39 310.26
[11,] 355.28 391.58 383.06 383.14 322.82 288.98 264.74 273.68 283.42 319.14
[12,] 364.18 389.77 383.14 381.07 324.87 295.93 310.72 270.00 283.06 316.56
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]   [,27]   [,28]   [,29]   [,30]
 [1,] 319.07 389.91 439.38 476.67 627.83 806.25  760.86 1127.04 1369.89 1738.11
 [2,] 331.92 406.95 442.08 510.10 629.79 803.20  816.09 1134.72 1390.55 1641.84
 [3,] 356.86 413.79 424.03 549.86 631.17 889.60  858.69 1117.96 1356.40 1652.21
 [4,] 358.97 404.88 423.35 555.00 664.75 922.30  943.00 1095.41 1372.73 1742.14
 [5,] 340.55 406.67 433.85 557.09 654.90 968.43  924.27 1113.34 1424.00 1673.77
 [6,] 328.18 403.26 429.23 610.65 679.37 909.71  890.20 1148.69 1479.76 1649.69
 [7,] 355.68 383.78 421.87 675.39 667.31 888.66  928.65 1205.43 1512.60 1591.19
 [8,] 356.35 392.37 430.66 596.15 655.66 889.49  945.67 1232.92 1528.66 1598.76
 [9,] 351.02 398.09 424.48 633.71 665.38 939.77  934.23 1192.97 1572.21 1589.90
[10,] 359.77 400.51 437.93 632.59 665.41 839.03  949.38 1215.81 1757.21 1630.31
[11,] 378.95 405.28 456.05 598.19 712.65 829.93  996.59 1270.98 1770.95 1744.81
[12,] 378.92 420.46 469.90 585.78 754.60 806.62 1043.16 1342.02 1665.21 1746.58
        [,31]
 [1,] 1721.64
 [2,]      NA
 [3,]      NA
 [4,]      NA
 [5,]      NA
 [6,]      NA
 [7,]      NA
 [8,]      NA
 [9,]      NA
[10,]      NA
[11,]      NA
[12,]      NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9] [,10]
 [1,]  29.61   6.64 -22.46   0.40  -6.64  18.67  -1.12  15.13  -4.78  1.56
 [2,]  36.79 -17.45 -16.75  25.58  17.34  -9.66 -15.04   0.72   6.69 -7.28
 [3,]   9.80  15.07  -3.69  -8.59  -7.14 -34.51 -16.50   6.70 -19.81 -0.54
 [4,] -71.39   8.30  14.40   6.81   7.79   1.54   2.64 -23.75  -0.49 -9.57
 [5,]  13.18 -12.89  13.30  -5.26  29.44   7.94  -5.75 -18.82  -4.96 -5.72
 [6,]   5.13  -3.97 -10.80   1.96  21.88  -0.48 -13.35  -5.19  -1.43 -1.38
 [7,] -25.17   0.30   0.50   0.30 -10.64   0.26  -3.45 -16.72   9.77  3.57
 [8,]  -9.93 -29.99   0.70   5.64   0.93 -13.70   7.44  10.20   0.97 11.91
 [9,]  13.29  -0.01  13.20  28.32  10.63  -6.35  -9.90  32.20 -11.38 -9.66
[10,]  -4.40  -7.50  -7.05  41.07  -0.77 -17.82  -3.39  -5.41  -7.57  2.37
[11,] -18.20   0.80   2.50   5.78   4.97  -6.68   5.13  -8.58   9.95 -1.05
[12,] -11.90   1.00  -4.35 -25.95   2.22  13.39  27.38   0.99   1.48 -9.36
       [,11] [,12] [,13] [,14]  [,15]  [,16]  [,17]  [,18] [,19] [,20]  [,21]
 [1,]  -0.20  9.47 -5.10  2.13  -8.85 -17.30  -2.44 -10.11  5.44 -0.31  12.85
 [2,]  -5.81  3.58 -0.74 12.01 -13.89   0.36  -4.60   1.25 -5.96  5.66  24.94
 [3,]   0.30 -4.97 -1.91  5.31  -8.53   8.39   0.25  15.54 -3.62 13.99   2.11
 [4,]   0.77  2.22  5.48 -8.55   5.23  -1.55  -1.37 -13.47  1.16 -1.44 -18.42
 [5,]  12.07 -6.86  8.91 -3.36  -7.34  12.35  -3.34  -6.70 -2.55  8.62 -12.37
 [6,]  25.03  4.16 -5.81 -0.92  -0.63  -9.19  -6.18  -4.50 11.88 11.81  27.50
 [7,]   4.71  4.21  2.34 -6.66  -3.08  -6.78 -15.13  10.54 -2.13  6.69   0.67
 [8,]  20.30 -0.15 -1.33 -1.80 -16.66   0.55  -5.23  -4.14 -2.70 -7.89  -5.33
 [9,] -13.35 -5.13 -2.56  3.88  -0.09  -8.76   0.61  -7.12  4.86 -3.03   8.75
[10,] -23.56 11.22 -0.61 -4.21  -1.19   4.87   8.05  -0.79 11.03  8.88  19.18
[11,]   8.90 -1.81  0.08 -2.07   2.05   6.95  45.98  -3.68 -0.36 -2.58  -0.03
[12,]   9.65 -5.38  2.17 -3.22 -18.83  -1.81 -17.54  -3.99 -6.90  2.51  10.99
       [,22]  [,23]  [,24]  [,25]   [,26]  [,27]  [,28]   [,29]  [,30] [,31]
 [1,]  17.04   2.70  33.43   1.96   -3.05  55.23   7.68   20.66 -96.27    NA
 [2,]   6.84 -18.05  39.76   1.38   86.40  42.60 -16.76  -34.15  10.37    NA
 [3,]  -8.91  -0.68   5.14  33.58   32.70  84.31 -22.55   16.33  89.93    NA
 [4,]   1.79  10.50   2.09  -9.85   46.13 -18.73  17.93   51.27 -68.37    NA
 [5,]  -3.41  -4.62  53.56  24.47  -58.72 -34.07  35.35   55.76 -24.08    NA
 [6,] -19.48  -7.36  64.74 -12.06  -21.05  38.45  56.74   32.84 -58.50    NA
 [7,]   8.59   8.79 -79.24 -11.65    0.83  17.02  27.49   16.06   7.57    NA
 [8,]   5.72  -6.18  37.56   9.72   50.28 -11.44 -39.95   43.55  -8.86    NA
 [9,]   2.42  13.45  -1.12   0.03 -100.74  15.15  22.84  185.00  40.41    NA
[10,]   4.77  18.12 -34.40  47.24   -9.10  47.21  55.17   13.74 114.50    NA
[11,]  15.18  13.85 -12.41  41.95  -23.31  46.57  71.04 -105.74   1.77    NA
[12,]  18.92   6.77  42.05  51.65  -45.76  83.88  27.87   72.90 -24.94    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/1p8wi1355659795.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/21rsa1355659795.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/3ipa51355659795.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/47opv1355659795.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,] 266.01 271.45 265.490 261.870 263.030 260.480 272.360 261.310 256.08
[2,] 338.51 331.92 347.480 338.890 340.550 338.620 342.400 340.810 347.71
[3,] 385.31 387.89 393.165 386.735 391.605 382.885 379.415 381.485 384.48
[4,] 472.12 486.24 481.290 491.090 443.610 451.550 460.230 451.330 450.52
[5,] 627.83 629.79 631.170 664.750 557.090 610.650 460.230 596.150 450.52
       [,10]  [,11]   [,12]
[1,] 256.690 264.74 270.000
[2,] 343.060 340.20 341.000
[3,] 382.015 386.23 382.105
[4,] 461.150 460.38 469.900
[5,] 632.590 598.19 585.780

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 347.3946 343.3737 354.5652 342.8303 361.8756 350.3084 345.4249 349.6036
[2,] 423.2254 432.4063 431.7648 430.6397 421.3344 415.4616 413.4051 413.3664
         [,9]    [,10]   [,11]    [,12]
[1,] 354.8227 347.9499 351.562 344.9216
[2,] 414.1373 416.0801 420.898 419.2884

$out
 [1]  806.25  760.86 1127.04 1369.89 1738.11 1721.64  803.20  816.09 1134.72
[10] 1390.55 1641.84  889.60  858.69 1117.96 1356.40 1652.21  922.30  943.00
[19] 1095.41 1372.73 1742.14  654.90  968.43  924.27 1113.34 1424.00 1673.77
[28]  679.37  909.71  890.20 1148.69 1479.76 1649.69  675.39  667.31  888.66
[37]  928.65 1205.43 1512.60 1591.19  655.66  889.49  945.67 1232.92 1528.66
[46] 1598.76  633.71  665.38  939.77  934.23 1192.97 1572.21 1589.90  665.41
[55]  839.03  949.38 1215.81 1757.21 1630.31  712.65  829.93  996.59 1270.98
[64] 1770.95 1744.81  754.60  806.62 1043.16 1342.02 1665.21 1746.58

$group
 [1]  1  1  1  1  1  1  2  2  2  2  2  3  3  3  3  3  4  4  4  4  4  5  5  5  5
[26]  5  5  6  6  6  6  6  6  7  7  7  7  7  7  7  8  8  8  8  8  8  9  9  9  9
[51]  9  9  9 10 10 10 10 10 10 11 11 11 11 11 11 12 12 12 12 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/5ycvp1355659795.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]
[1,] -22.46 -34.150 -34.510 -23.75 -34.070 -21.05 -25.170 -29.99 -13.35 -23.56
[2,]  -4.78  -9.660  -7.140  -9.57  -6.860  -7.36  -6.660  -7.89  -6.35  -5.41
[3,]   1.76   0.985   0.275   1.35  -3.385  -1.15   0.585  -0.74   0.32   0.88
[4,]  12.85  12.010  13.990   6.81  12.350  11.88   7.570   7.44  13.29  13.74
[5,]  33.43  42.600  33.580  17.93  35.350  38.45  27.490  20.30  40.41  41.07
       [,11]  [,12]
[1,] -18.200 -25.95
[2,]  -2.580  -6.90
[3,]   1.285   1.24
[4,]   8.900  13.39
[5,]  15.180  42.05

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -3.325677 -5.266085 -5.820312 -3.375093 -8.926455 -6.700109 -3.519888
[2,]  6.845677  7.236085  6.370312  6.075093  2.156455  4.400109  4.689888
          [,8]      [,9]     [,10]     [,11]  [,12]
[1,] -5.162202 -5.345496 -4.644147 -2.026604 -4.613
[2,]  3.682202  5.985496  6.404147  4.596604  7.093

$out
 [1]   55.23  -96.27   86.40   84.31   89.93  -71.39   46.13   51.27  -68.37
[10]   53.56  -58.72   55.76   64.74   56.74  -58.50  -79.24   37.56   50.28
[19]  -39.95   43.55 -100.74  185.00  -34.40   47.24   47.21   55.17  114.50
[28]   45.98   41.95  -23.31   46.57   71.04 -105.74   51.65  -45.76   83.88
[37]   72.90

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

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

Warning message:
In bxp(list(stats = c(-22.46, -4.78000000000003, 1.75999999999996,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/60f541355659795.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,] 393.600 340.200 299.100 321.500 390.920 406.780 361.75 352.330 348.740
[2,] 412.320 347.705 314.750 339.665 404.690 434.455 367.24 371.645 357.665
[3,] 417.950 377.550 318.370 344.200 443.970 447.340 379.72 391.190 363.585
[4,] 441.255 383.830 326.325 362.500 460.305 459.560 397.08 402.060 372.320
[5,] 481.290 394.260 342.000 376.660 465.350 486.240 420.17 416.810 383.640
       [,10]   [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]  [,18]
[1,] 337.240 329.010 373.830 376.640 381.070 322.820 284.110 256.080 270.00
[2,] 341.935 332.450 380.895 380.705 384.370 324.485 289.040 263.025 274.83
[3,] 344.905 348.715 383.715 383.405 387.395 344.155 294.400 284.290 282.33
[4,] 354.180 369.535 386.260 385.725 394.530 353.460 298.295 289.505 286.06
[5,] 361.730 392.190 391.580 391.030 404.760 377.850 308.290 310.720 299.86
       [,19]   [,20]   [,21]  [,22]   [,23]   [,24]   [,25]   [,26]    [,27]
[1,] 260.480 275.850 319.070 383.78 421.870 476.670 627.830 803.200  816.090
[2,] 264.260 287.785 336.235 395.23 424.255 552.430 643.035 818.275  874.445
[3,] 268.880 306.470 356.015 404.07 432.255 590.965 665.065 889.075  931.440
[4,] 272.375 315.525 359.370 406.81 440.730 621.620 673.340 916.005  947.525
[5,] 283.420 321.180 378.950 420.46 456.050 675.390 712.650 968.430 1043.160
        [,28]   [,29]    [,30]   [,31]
[1,] 1095.410 1356.40 1589.900 1721.64
[2,] 1122.500 1381.64 1614.535 1721.64
[3,] 1170.830 1496.18 1650.950 1721.64
[4,] 1224.365 1618.71 1740.125 1721.64
[5,] 1342.020 1770.95 1746.580 1721.64

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 404.7526 361.0731 313.0906 333.7848 418.6036 435.8894 366.1098 377.3175
[2,] 431.1474 394.0269 323.6494 354.6152 469.3364 458.7906 393.3302 405.0625
         [,9]  [,10]    [,11]   [,12]    [,13]   [,14]    [,15]    [,16]
[1,] 356.9008 339.32 331.8003 381.268 381.1153 382.761 330.9393 290.1787
[2,] 370.2692 350.49 365.6297 386.162 385.6947 392.029 357.3707 298.6213
        [,17]    [,18]    [,19]    [,20]   [,21]    [,22]    [,23]   [,24]
[1,] 272.2123 277.2079 265.1787 293.8176 345.463 398.7883 424.7406 559.407
[2,] 296.3677 287.4521 272.5813 319.1224 366.567 409.3517 439.7694 622.523
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30]   [,31]
[1,] 651.2427 844.4997 898.1077 1124.369 1388.051 1593.668 1721.64
[2,] 678.8873 933.6503 964.7723 1217.291 1604.309 1708.232 1721.64

$out
[1] 491.09 417.73 423.51 469.90 754.60 760.86

$group
[1]  1  4  4 23 25 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(393.6, 412.32, 417.95, 441.255, 481.29, 340.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/7w3ah1355659795.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,] 509.6813 379.4150 391.3338
[2,] 520.3412 382.0600 396.8356
[3,] 524.2265 384.8950 400.3762
[4,] 536.7518 387.3125 405.3094
[5,] 546.7777 393.1650 414.6675

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 516.7415 382.4993 396.5113
[2,] 531.7115 387.2907 404.2412

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(509.681333333333, 520.341166666667, 524.2265,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1p8wi1355659795.ps tmp/1p8wi1355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/21rsa1355659795.ps tmp/21rsa1355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ipa51355659795.ps tmp/3ipa51355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/47opv1355659795.ps tmp/47opv1355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ycvp1355659795.ps tmp/5ycvp1355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/60f541355659795.ps tmp/60f541355659795.png",intern=TRUE))
character(0)
> try(system("convert tmp/7w3ah1355659795.ps tmp/7w3ah1355659795.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.678   1.016   4.713