R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.

  Natural language support but running in an English locale

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)
> 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] 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/html/freestat/rcomp/tmp/1qmpc1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/2qmpc1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/3jeox1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/4jeox1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/5jeox1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/6cnni1289315573.ps",horizontal=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/html/freestat/rcomp/tmp/7cnni1289315573.ps",horizontal=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/1qmpc1289315573.ps tmp/1qmpc1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/2qmpc1289315573.ps tmp/2qmpc1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jeox1289315573.ps tmp/3jeox1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jeox1289315573.ps tmp/4jeox1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jeox1289315573.ps tmp/5jeox1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/6cnni1289315573.ps tmp/6cnni1289315573.png",intern=TRUE))
character(0)
> try(system("convert tmp/7cnni1289315573.ps tmp/7cnni1289315573.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.020   1.472   2.288