R version 2.12.1 (2010-12-16)
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(684.28
+ ,722.57
+ ,695.96
+ ,688.13
+ ,720.76
+ ,737.26
+ ,736.9
+ ,727.38
+ ,728.83
+ ,715.58
+ ,735.93
+ ,758.69
+ ,758.69
+ ,740.99
+ ,719.44
+ ,721.24
+ ,737.38
+ ,756.52
+ ,745.32
+ ,733.76
+ ,738.46
+ ,736.65
+ ,737.5
+ ,725.22
+ ,722.45
+ ,720.44
+ ,733.6
+ ,662.93
+ ,763.87
+ ,746.67
+ ,712.59
+ ,655.02
+ ,660.2
+ ,633.59
+ ,627.08
+ ,635.63
+ ,668.87
+ ,656.7
+ ,631.9
+ ,610.83
+ ,632.98
+ ,631.06
+ ,644.3
+ ,647.55
+ ,655.14
+ ,674.4
+ ,676.21
+ ,670.43
+ ,683.43
+ ,703.42
+ ,708.48
+ ,714.5
+ ,706.19
+ ,694.15
+ ,658.15
+ ,648.4
+ ,630.82
+ ,634.19
+ ,658.15
+ ,636.11
+ ,640.33
+ ,601.8
+ ,611.07
+ ,634.31
+ ,625.76
+ ,596.86
+ ,604.09
+ ,587.11
+ ,583.13
+ ,579.16
+ ,551.47
+ ,541.83
+ ,569.53
+ ,564.71
+ ,572.06
+ ,578.56
+ ,595.78
+ ,567.84
+ ,534.25
+ ,519.08
+ ,531.96
+ ,551.83
+ ,560.62
+ ,579.52
+ ,593.13
+ ,626.72
+ ,715.1
+ ,832.38
+ ,837.2
+ ,879.46
+ ,881.99
+ ,1044.9
+ ,924.73
+ ,987.95
+ ,966.15
+ ,1016.72
+ ,1107.27
+ ,1296.67
+ ,1379.75
+ ,1543.03
+ ,1570.12
+ ,1538.81
+ ,1484.63
+ ,1451.27
+ ,1414.91
+ ,1456.93
+ ,1319.19
+ ,1267.89
+ ,1349.77
+ ,1240.2
+ ,1189.51
+ ,1117.87
+ ,1080.06
+ ,1054.77
+ ,1019.47
+ ,1049.96
+ ,1060.79
+ ,1070.42
+ ,1075.24
+ ,1080.06
+ ,1074.04
+ ,1062
+ ,1064.4
+ ,1071.63
+ ,1057.18
+ ,898.24
+ ,895.83
+ ,951.22
+ ,936.77
+ ,901.85
+ ,888.61
+ ,870.55
+ ,887.41
+ ,596.02
+ ,586.39
+ ,596.02
+ ,721
+ ,777.51
+ ,723.48
+ ,680.64
+ ,613.45
+ ,558.18
+ ,641.49
+ ,652.19
+ ,619.37
+ ,655.83
+ ,667.93
+ ,667.7
+ ,663.07
+ ,633.89
+ ,595.28
+ ,568.94
+ ,572.72
+ ,535.26
+ ,508.04
+ ,512.94
+ ,495.22
+ ,469.37
+ ,469.37
+ ,429.69
+ ,468.13
+ ,470.06
+ ,464.93
+ ,450.74
+ ,423.51
+ ,454.84
+ ,497.77
+ ,465.45
+ ,542.31
+ ,606.03
+ ,609.58
+ ,645.79
+ ,719.63
+ ,779.41
+ ,773.5
+ ,806.82
+ ,876.1
+ ,824.64
+ ,881.7
+ ,878.41
+ ,904.18
+ ,892.34
+ ,887.13
+ ,867.85
+ ,839.28
+ ,826.06
+ ,751.11
+ ,789.25
+ ,732.98
+ ,622.07
+ ,600.95
+ ,590.53
+ ,584.39
+ ,525.31
+ ,573.83
+ ,597.67
+ ,743.54
+ ,701.36
+ ,671.43
+ ,751.65
+ ,738.33
+ ,681.61
+ ,616.97
+ ,632.94
+ ,677.73
+ ,730.96
+ ,719.66
+ ,764.21
+ ,805
+ ,829.35
+ ,826.26
+ ,765.93
+ ,801.91
+ ,769.16
+ ,739.89
+ ,688.07
+ ,636.11
+ ,631.72
+ ,625.92
+ ,627.82
+ ,606.13
+ ,595.3
+ ,583.14
+ ,500.19
+ ,462.89
+ ,417.47
+ ,472.27
+ ,474.81
+ ,489.07
+ ,493.14
+ ,626.64
+ ,680.43
+ ,620.3
+ ,676.74
+ ,690.03
+ ,631.04
+ ,623.26
+ ,619.83
+ ,631.74
+ ,648.77
+ ,724.21
+ ,727.09
+ ,767.31
+ ,801.42
+ ,817.72
+ ,764.33
+ ,746.93
+ ,717.29
+ ,695.9
+ ,688.38
+ ,663.53
+ ,688.39
+ ,716.14
+ ,733.28
+ ,688.23
+ ,760.63
+ ,716.77
+ ,683.81
+ ,630.79
+ ,617.91
+ ,524.19
+ ,441.98
+ ,466.09
+ ,501.44
+ ,599.55
+ ,621.99
+ ,607.57
+ ,614.56
+ ,619.09
+ ,603.14
+ ,569.12
+ ,575.72
+ ,642.41
+ ,748.18
+ ,768.39
+ ,763.16
+ ,800.63
+ ,778.49
+ ,733.7
+ ,740.17
+ ,678.33
+ ,697.09
+ ,678.54
+ ,670.87
+ ,674.52
+ ,664.39
+ ,646.52
+ ,661.25
+ ,729.36
+ ,721.52
+ ,708.75
+ ,706.12
+ ,676.93
+ ,708.15
+ ,717.64
+ ,714.35
+ ,703.6
+ ,718.96
+ ,736.03
+ ,717.96
+ ,730.77
+ ,734.66
+ ,728.58
+ ,729.18
+ ,717.13
+ ,684.6
+ ,692.06
+ ,668.84
+ ,647.4
+ ,622.57
+ ,593.69
+ ,583.24
+ ,639.17
+ ,709.73
+ ,698.75
+ ,697.7
+ ,732.25
+ ,714.49
+ ,704.85
+ ,717.56
+ ,704.91
+ ,690.74
+ ,790.34
+ ,790.69
+ ,893.66
+ ,875.62
+ ,914.02
+ ,940.37
+ ,989.08
+ ,987.73
+ ,936.34
+ ,914.11
+ ,943.08
+ ,1080.64
+ ,1102.71
+ ,1097.26
+ ,1157.37
+ ,1154.25
+ ,1136
+ ,1133.42
+ ,1062.44
+ ,1041.1
+ ,1047.97
+ ,941.71
+ ,896.79
+ ,734.91
+ ,646.17
+ ,666.39
+ ,625.45
+ ,586.26
+ ,617.15
+ ,686.86
+ ,761.34
+ ,741.73
+ ,763.62
+ ,822.57
+ ,867.58
+ ,944.85
+ ,953.95
+ ,970.08
+ ,1003.22
+ ,972.81
+ ,1003.86
+ ,991.83
+ ,919.13
+ ,919.52
+ ,916.17
+ ,936.18
+ ,960.65)
> 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,] 684.28 758.69 722.45 668.87 683.43 640.33 569.53  593.13 1107.27 1349.77
 [2,] 722.57 740.99 720.44 656.70 703.42 601.80 564.71  626.72 1296.67 1240.20
 [3,] 695.96 719.44 733.60 631.90 708.48 611.07 572.06  715.10 1379.75 1189.51
 [4,] 688.13 721.24 662.93 610.83 714.50 634.31 578.56  832.38 1543.03 1117.87
 [5,] 720.76 737.38 763.87 632.98 706.19 625.76 595.78  837.20 1570.12 1080.06
 [6,] 737.26 756.52 746.67 631.06 694.15 596.86 567.84  879.46 1538.81 1054.77
 [7,] 736.90 745.32 712.59 644.30 658.15 604.09 534.25  881.99 1484.63 1019.47
 [8,] 727.38 733.76 655.02 647.55 648.40 587.11 519.08 1044.90 1451.27 1049.96
 [9,] 728.83 738.46 660.20 655.14 630.82 583.13 531.96  924.73 1414.91 1060.79
[10,] 715.58 736.65 633.59 674.40 634.19 579.16 551.83  987.95 1456.93 1070.42
[11,] 735.93 737.50 627.08 676.21 658.15 551.47 560.62  966.15 1319.19 1075.24
[12,] 758.69 725.22 635.63 670.43 636.11 541.83 579.52 1016.72 1267.89 1080.06
        [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 1074.04 887.41 619.37 495.22 542.31 904.18 584.39 677.73 636.11 489.07
 [2,] 1062.00 596.02 655.83 469.37 606.03 892.34 525.31 730.96 631.72 493.14
 [3,] 1064.40 586.39 667.93 469.37 609.58 887.13 573.83 719.66 625.92 626.64
 [4,] 1071.63 596.02 667.70 429.69 645.79 867.85 597.67 764.21 627.82 680.43
 [5,] 1057.18 721.00 663.07 468.13 719.63 839.28 743.54 805.00 606.13 620.30
 [6,]  898.24 777.51 633.89 470.06 779.41 826.06 701.36 829.35 595.30 676.74
 [7,]  895.83 723.48 595.28 464.93 773.50 751.11 671.43 826.26 583.14 690.03
 [8,]  951.22 680.64 568.94 450.74 806.82 789.25 751.65 765.93 500.19 631.04
 [9,]  936.77 613.45 572.72 423.51 876.10 732.98 738.33 801.91 462.89 623.26
[10,]  901.85 558.18 535.26 454.84 824.64 622.07 681.61 769.16 417.47 619.83
[11,]  888.61 641.49 508.04 497.77 881.70 600.95 616.97 739.89 472.27 631.74
[12,]  870.55 652.19 512.94 465.45 878.41 590.53 632.94 688.07 474.81 648.77
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]   [,28]   [,29]  [,30]
 [1,] 724.21 716.14 599.55 800.63 729.36 730.77 639.17  893.66 1157.37 625.45
 [2,] 727.09 733.28 621.99 778.49 721.52 734.66 709.73  875.62 1154.25 586.26
 [3,] 767.31 688.23 607.57 733.70 708.75 728.58 698.75  914.02 1136.00 617.15
 [4,] 801.42 760.63 614.56 740.17 706.12 729.18 697.70  940.37 1133.42 686.86
 [5,] 817.72 716.77 619.09 678.33 676.93 717.13 732.25  989.08 1062.44 761.34
 [6,] 764.33 683.81 603.14 697.09 708.15 684.60 714.49  987.73 1041.10 741.73
 [7,] 746.93 630.79 569.12 678.54 717.64 692.06 704.85  936.34 1047.97 763.62
 [8,] 717.29 617.91 575.72 670.87 714.35 668.84 717.56  914.11  941.71 822.57
 [9,] 695.90 524.19 642.41 674.52 703.60 647.40 704.91  943.08  896.79 867.58
[10,] 688.38 441.98 748.18 664.39 718.96 622.57 690.74 1080.64  734.91 944.85
[11,] 663.53 466.09 768.39 646.52 736.03 593.69 790.34 1102.71  646.17 953.95
[12,] 688.39 501.44 763.16 661.25 717.96 583.24 790.69 1097.26  666.39 970.08
        [,31]
 [1,] 1003.22
 [2,]  972.81
 [3,] 1003.86
 [4,]  991.83
 [5,]  919.13
 [6,]  919.52
 [7,]  916.17
 [8,]  936.18
 [9,]  960.65
[10,]      NA
[11,]      NA
[12,]      NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]    [,8]    [,9]   [,10]
 [1,]  38.29 -17.70  -2.01 -12.17  19.99 -38.53  -4.82   33.59  189.40 -109.57
 [2,] -26.61 -21.55  13.16 -24.80   5.06   9.27   7.35   88.38   83.08  -50.69
 [3,]  -7.83   1.80 -70.67 -21.07   6.02  23.24   6.50  117.28  163.28  -71.64
 [4,]  32.63  16.14 100.94  22.15  -8.31  -8.55  17.22    4.82   27.09  -37.81
 [5,]  16.50  19.14 -17.20  -1.92 -12.04 -28.90 -27.94   42.26  -31.31  -25.29
 [6,]  -0.36 -11.20 -34.08  13.24 -36.00   7.23 -33.59    2.53  -54.18  -35.30
 [7,]  -9.52 -11.56 -57.57   3.25  -9.75 -16.98 -15.17  162.91  -33.36   30.49
 [8,]   1.45   4.70   5.18   7.59 -17.58  -3.98  12.88 -120.17  -36.36   10.83
 [9,] -13.25  -1.81 -26.61  19.26   3.37  -3.97  19.87   63.22   42.02    9.63
[10,]  20.35   0.85  -6.51   1.81  23.96 -27.69   8.79  -21.80 -137.74    4.82
[11,]  22.76 -12.28   8.55  -5.78 -22.04  -9.64  18.90   50.57  -51.30    4.82
[12,]   0.00  -2.77  33.24  13.00   4.22  27.70  13.61   90.55   81.88   -6.02
        [,11]   [,12]  [,13]  [,14]  [,15]   [,16]  [,17]  [,18]  [,19]  [,20]
 [1,]  -12.04 -291.39  36.46 -25.85  63.72  -11.84 -59.08  53.23  -4.39   4.07
 [2,]    2.40   -9.63  12.10   0.00   3.55   -5.21  48.52 -11.30  -5.80 133.50
 [3,]    7.23    9.63  -0.23 -39.68  36.21  -19.28  23.84  44.55   1.90  53.79
 [4,]  -14.45  124.98  -4.63  38.44  73.84  -28.57 145.87  40.79 -21.69 -60.13
 [5,] -158.94   56.51 -29.18   1.93  59.78  -13.22 -42.18  24.35 -10.83  56.44
 [6,]   -2.41  -54.03 -38.61  -5.13  -5.91  -74.95 -29.93  -3.09 -12.16  13.29
 [7,]   55.39  -42.84 -26.34 -14.19  33.32   38.14  80.22 -60.33 -82.95 -58.99
 [8,]  -14.45  -67.19   3.78 -27.23  69.28  -56.27 -13.32  35.98 -37.30  -7.78
 [9,]  -34.92  -55.27 -37.46  31.33 -51.46 -110.91 -56.72 -32.75 -45.42  -3.43
[10,]  -13.24   83.31 -27.22  42.93  57.06  -21.12 -64.64 -29.27  54.80  11.91
[11,]  -18.06   10.70   4.90 -32.32  -3.29  -10.42  15.97 -51.82   2.54  17.03
[12,]   16.86  -32.82 -17.72  76.86  25.77   -6.14  44.79 -51.96  14.26  75.44
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]   [,29]  [,30]
 [1,]   2.88  17.14  22.44 -22.14  -7.84   3.89  70.56 -18.04   -3.12 -39.19
 [2,]  40.22 -45.05 -14.42 -44.79 -12.77  -6.08 -10.98  38.40  -18.25  30.89
 [3,]  34.11  72.40   6.99   6.47  -2.63   0.60  -1.05  26.35   -2.58  69.71
 [4,]  16.30 -43.86   4.53 -61.84 -29.19 -12.05  34.55  48.71  -70.98  74.48
 [5,] -53.39 -32.96 -15.95  18.76  31.22 -32.53 -17.76  -1.35  -21.34 -19.61
 [6,] -17.40 -53.02 -34.02 -18.55   9.49   7.46  -9.64 -51.39    6.87  21.89
 [7,] -29.64 -12.88   6.60  -7.67  -3.29 -23.22  12.71 -22.23 -106.26  58.95
 [8,] -21.39 -93.72  66.69   3.65 -10.75 -21.44 -12.65  28.97  -44.92  45.01
 [9,]  -7.52 -82.21 105.77 -10.13  15.36 -24.83 -14.17 137.56 -161.88  77.27
[10,] -24.85  24.11  20.21 -17.87  17.07 -28.88  99.60  22.07  -88.74   9.10
[11,]  24.86  35.35  -5.23  14.73 -18.07 -10.45   0.35  -5.45   20.22  16.13
[12,]  27.75  98.11  37.47  68.11  12.81  55.93 102.97  60.11  -40.94  33.14
       [,31]
 [1,] -30.41
 [2,]  31.05
 [3,] -12.03
 [4,] -72.70
 [5,]   0.39
 [6,]  -3.35
 [7,]  20.01
 [8,]  24.47
 [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/1gtfe1321397510.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/2nlmy1321397510.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/3p2q51321397510.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/4x9jf1321397510.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,]  489.07 469.37 469.370  429.69  468.13  470.060  464.930  450.740  423.51
[2,]  622.41 614.01 621.535  640.05  670.00  680.275  651.225  639.295  627.04
[3,]  684.28 720.44 708.480  706.12  721.00  737.260  717.640  717.290  703.60
[4,]  844.02 759.74 750.505  816.90  827.46  827.705  799.880  814.695  871.84
[5,] 1157.37 972.81 914.020 1071.63 1062.44 1041.100 1019.470 1049.960 1060.79
       [,10]  [,11]   [,12]
[1,] 417.470 466.09  465.45
[2,] 619.830 600.95  590.53
[3,] 684.995 660.84  668.41
[4,] 769.160 790.34  790.69
[5,] 987.950 966.15 1080.06

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 621.3923 679.0853 671.8814 655.9341 676.3166 695.4228 675.4552 667.5156
[2,] 747.1677 761.7947 745.0786 756.3059 765.6834 779.0972 759.8248 767.0644
         [,9]    [,10]    [,11]    [,12]
[1,] 634.1315 641.9182 606.2072 610.6704
[2,] 773.0685 728.0718 715.4728 726.1496

$out
 [1] 1349.77 1296.67 1240.20 1062.00 1154.25 1379.75 1189.51 1064.40 1136.00
[10] 1003.86 1543.03 1117.87 1133.42 1570.12 1080.06 1538.81 1054.77 1484.63
[19] 1047.97 1451.27 1414.91 1456.93 1070.42 1080.64 1319.19 1075.24 1102.71
[28] 1267.89 1097.26

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

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

Warning message:
In bxp(list(stats = c(489.07, 622.41, 684.28, 844.02, 1157.37, 469.37,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5aexy1321397510.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,] -59.080 -50.690 -39.680 -72.700 -53.39 -74.950 -82.95 -67.190 -110.910
[2,] -20.090 -13.595  -2.605 -25.130 -28.42 -34.690 -27.99 -24.335  -37.460
[3,]  -4.390   0.000   6.470   4.820 -13.22 -11.200 -11.56  -7.780   -8.825
[4,]  21.215  22.025  30.230  36.495  17.63   1.085  16.36   9.210   19.260
[5,]  70.560  48.520  72.400 124.980  59.78  21.890  80.22  45.010   77.270
       [,10]   [,11]   [,12]
[1,] -88.740 -32.320 -51.960
[2,] -24.850 -10.450   0.000
[3,]   3.315   1.445  26.735
[4,]  22.070  16.130  60.110
[5,]  83.310  50.570 102.970

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

$conf
           [,1]      [,2]      [,3]      [,4]        [,5]       [,6]       [,7]
[1,] -16.111383 -10.10811 -2.847797 -12.66772 -26.2879022 -21.352100 -24.145482
[2,]   7.331383  10.10811 15.787797  22.30772  -0.1520978  -1.047900   1.025482
           [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -17.299279 -25.18686 -10.21988 -6.222459  9.395236
[2,]   1.739279   7.53686  16.84988  9.112459 44.074764

$out
 [1]  189.40 -109.57 -291.39   88.38   83.08  133.50  -70.67  117.28  163.28
[10]  -71.64  145.87 -158.94  162.91 -106.26 -120.17   69.28  -93.72   66.69
[19]  105.77  137.56 -161.88 -137.74   99.60  -51.30  -51.82

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

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

Warning message:
In bxp(list(stats = c(-59.08, -20.0900000000000, -4.38999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6f9sa1321397510.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,] 684.280 719.440 627.080 610.830 630.820 541.830 519.080  593.130 1107.27
[2,] 705.770 729.490 645.325 632.440 642.255 581.145 543.040  773.740 1307.93
[3,] 724.975 737.440 687.760 651.345 670.790 599.330 566.275  880.725 1433.09
[4,] 736.415 743.155 728.025 669.650 704.805 618.415 575.310  977.050 1511.72
[5,] 758.690 758.690 763.870 676.210 714.500 640.330 595.780 1044.900 1570.12
       [,10]    [,11]  [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18]
[1,] 1019.47  870.550 558.18 508.040 429.690 542.310 590.530 525.310 677.730
[2,] 1057.78  897.035 596.02 552.100 452.790 627.685 677.525 591.030 725.310
[3,] 1077.65  943.995 646.84 607.325 466.790 776.455 807.655 652.185 765.070
[4,] 1153.69 1063.200 722.24 659.450 469.715 850.370 877.490 719.845 803.455
[5,] 1240.20 1074.040 887.41 667.930 470.060 881.700 904.180 751.650 829.350
      [,19]   [,20]   [,21]   [,22]   [,23]   [,24]   [,25]   [,26]   [,27]
[1,] 417.47 619.830 663.530 441.980 569.120 646.520 703.600 583.240 690.740
[2,] 473.54 620.065 692.145 512.815 601.345 667.630 707.135 634.985 698.225
[3,] 589.22 628.840 725.650 657.300 616.825 678.435 715.995 688.330 707.320
[4,] 626.87 662.755 765.820 716.455 695.295 736.935 720.240 728.880 724.905
[5,] 636.11 690.030 817.720 760.630 768.390 800.630 736.030 734.660 732.250
        [,28]    [,29]   [,30]   [,31]
[1,]  875.620  646.170 586.260  916.17
[2,]  914.065  815.850 656.155  919.52
[3,]  941.725 1044.535 762.480  960.65
[4,] 1034.860 1134.710 906.215  991.83
[5,] 1102.710 1157.370 970.080 1003.86

$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,] 710.9976 731.2073 650.04 634.3733 642.2605 582.3309 551.5564 787.994
[2,] 738.9524 743.6727 725.48 668.3167 699.3195 616.3291 580.9936 973.456
        [,9]    [,10]    [,11]    [,12]   [,13]    [,14]    [,15]    [,16]
[1,] 1340.14 1033.905  868.206 589.2702 558.362 459.0704 674.8869 716.4496
[2,] 1526.04 1121.395 1019.784 704.4098 656.288 474.5096 878.0231 898.8604
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 593.4316 729.4275 519.2852 609.3688 692.0463 564.4184 573.9738 646.8245
[2,] 710.9384 800.7125 659.1548 648.3112 759.2537 750.1816 659.6762 710.0455
        [,25]    [,26]    [,27]    [,28]     [,29]   [,30]    [,31]
[1,] 710.0177 645.5039 695.1511 886.6296  899.1008 648.426 922.5667
[2,] 721.9723 731.1561 719.4889 996.8204 1189.9692 876.534 998.7333

$out
 [1] 1349.77  495.22  423.51  497.77  489.07  493.14  676.93  639.17  790.34
[10]  790.69

$group
 [1] 10 14 14 14 20 20 25 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(684.28, 705.77, 724.975, 736.415, 758.69,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7oclz1321397510.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,] 725.1463 660.8400 686.0200
[2,] 733.2918 684.6375 692.3112
[3,] 751.0495 707.3000 726.2738
[4,] 762.1865 719.0400 740.9725
[5,] 780.7603 737.2600 753.9900

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 737.8705 691.6088 704.0790
[2,] 764.2285 722.9912 748.4685

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(725.146333333333, 733.291838709677, 751.049516129032,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1gtfe1321397510.ps tmp/1gtfe1321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nlmy1321397510.ps tmp/2nlmy1321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/3p2q51321397510.ps tmp/3p2q51321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/4x9jf1321397510.ps tmp/4x9jf1321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/5aexy1321397510.ps tmp/5aexy1321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/6f9sa1321397510.ps tmp/6f9sa1321397510.png",intern=TRUE))
character(0)
> try(system("convert tmp/7oclz1321397510.ps tmp/7oclz1321397510.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.284   0.480   2.739