R version 2.10.1 (2009-12-14)
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(93.2
+ ,96
+ ,95.2
+ ,77.1
+ ,70.9
+ ,64.8
+ ,70.1
+ ,77.3
+ ,79.5
+ ,100.6
+ ,100.7
+ ,107.1
+ ,95.9
+ ,82.8
+ ,83.3
+ ,80
+ ,80.4
+ ,67.5
+ ,75.7
+ ,71.1
+ ,89.3
+ ,101.1
+ ,105.2
+ ,114.1
+ ,96.3
+ ,84.4
+ ,91.2
+ ,81.9
+ ,80.5
+ ,70.4
+ ,74.8
+ ,75.9
+ ,86.3
+ ,98.7
+ ,100.9
+ ,113.8
+ ,89.8
+ ,84.4
+ ,87.2
+ ,85.6
+ ,72
+ ,69.2
+ ,77.5
+ ,78.1
+ ,94.3
+ ,97.7
+ ,100.2
+ ,116.4
+ ,97.1
+ ,93
+ ,96
+ ,80.5
+ ,76.1
+ ,69.9
+ ,73.6
+ ,92.6
+ ,94.2
+ ,93.5
+ ,108.5
+ ,109.4
+ ,105.1
+ ,92.5
+ ,97.1
+ ,81.4
+ ,79.1
+ ,72.1
+ ,78.7
+ ,87.1
+ ,91.4
+ ,109.9
+ ,116.3
+ ,113
+ ,100
+ ,84.8
+ ,94.3
+ ,87.1
+ ,90.3
+ ,72.4
+ ,84.9
+ ,92.7
+ ,92.2
+ ,114.9
+ ,112.5
+ ,118.3
+ ,106
+ ,91.2
+ ,96.6
+ ,96.3
+ ,88.2
+ ,70.2
+ ,86.5
+ ,88.2
+ ,102.8
+ ,119.1
+ ,119.2
+ ,125.1
+ ,106.1
+ ,102.1
+ ,105.2
+ ,101
+ ,84.3
+ ,87.5
+ ,92.7
+ ,94.4
+ ,113
+ ,113.9
+ ,122.9
+ ,132.7
+ ,106.9
+ ,96.6
+ ,127.3
+ ,98.2
+ ,100.2
+ ,89.4
+ ,95.3
+ ,104.2
+ ,106.4
+ ,116.2
+ ,135.9
+ ,134
+ ,104.6
+ ,107.1
+ ,123.5
+ ,98.8
+ ,98.6
+ ,90.6
+ ,89.1
+ ,105.2
+ ,114
+ ,122.1
+ ,138
+ ,142.2
+ ,116.4
+ ,112.6
+ ,123.8
+ ,103.6
+ ,113.9
+ ,98.6
+ ,95
+ ,116
+ ,113.9
+ ,127.5
+ ,131.4
+ ,145.9
+ ,131.5
+ ,131
+ ,130.5
+ ,118.9
+ ,114.3
+ ,85.7
+ ,104.6
+ ,105.1
+ ,117.3
+ ,142.5
+ ,140
+ ,159.8
+ ,131.2
+ ,125.4
+ ,126.5
+ ,119.4
+ ,113.5
+ ,98.7
+ ,114.5
+ ,113.8
+ ,133.1
+ ,143.4
+ ,137.3
+ ,165.2
+ ,126.9
+ ,124
+ ,135.7
+ ,130
+ ,109.4
+ ,117.8
+ ,120.3
+ ,121
+ ,132.3
+ ,142.9
+ ,147.4
+ ,175.9
+ ,132.6
+ ,123.7
+ ,153.3
+ ,134
+ ,119.6
+ ,116.2
+ ,118.6
+ ,130.7
+ ,129.3
+ ,144.4
+ ,163.2
+ ,179.4
+ ,128.1
+ ,138.4
+ ,152.7
+ ,120
+ ,140.5
+ ,116.2
+ ,121.4
+ ,127.8
+ ,143.6
+ ,157.6
+ ,166.2
+ ,182.3
+ ,153.1
+ ,147.6
+ ,157.7
+ ,137.2
+ ,151.5
+ ,98.7
+ ,145.8
+ ,151.7
+ ,129.4
+ ,174.1
+ ,197
+ ,193.9
+ ,164.1
+ ,142.8
+ ,157.9
+ ,159.2
+ ,162.2
+ ,123.1
+ ,130
+ ,150.1
+ ,169.4
+ ,179.7
+ ,182.1
+ ,194.3
+ ,161.4
+ ,169.4
+ ,168.8
+ ,158.1
+ ,158.5
+ ,135.3
+ ,149.3
+ ,143.4
+ ,142.2
+ ,188.4
+ ,166.2
+ ,199.2
+ ,182.7
+ ,145.2
+ ,182.1
+ ,158.7
+ ,141.6
+ ,132.6
+ ,139.6
+ ,147
+ ,166.6
+ ,157
+ ,180.4
+ ,210.2
+ ,159.8
+ ,157.8
+ ,168.2
+ ,158.4
+ ,152
+ ,142.2
+ ,137.2
+ ,152.6
+ ,166.8
+ ,165.6
+ ,198.6
+ ,201.5
+ ,170.7
+ ,164.4
+ ,179.7
+ ,157
+ ,168
+ ,139.3
+ ,138.6
+ ,153.4
+ ,138.9
+ ,172.1
+ ,198.4
+ ,217.8
+ ,173.7
+ ,153.8
+ ,175.6
+ ,147.1
+ ,160.3
+ ,135.2
+ ,148.8
+ ,151
+ ,148.2
+ ,182.2
+ ,189.2
+ ,183.1
+ ,170
+ ,158.4
+ ,176.1
+ ,156.2
+ ,153.2
+ ,117.9
+ ,149.8
+ ,156.6
+ ,166.7
+ ,156.8
+ ,158.6
+ ,210.8
+ ,203.6
+ ,175.2
+ ,168.7
+ ,155.9
+ ,147.3
+ ,137
+ ,141.1
+ ,167.4
+ ,160.2
+ ,191.9
+ ,174.4
+ ,208.2
+ ,159.4
+ ,161.1
+ ,172.1
+ ,158.4
+ ,114.6
+ ,159.6
+ ,159.7
+ ,159.4
+ ,160.7
+ ,165.5
+ ,205
+ ,205.2
+ ,141.6
+ ,148.1
+ ,184.9
+ ,132.5
+ ,137.3
+ ,135.5
+ ,121.7
+ ,166.1
+ ,146.8
+ ,162.8
+ ,186.8
+ ,185.5
+ ,151.5
+ ,158.1
+ ,143
+ ,151.2
+ ,147.6
+ ,130.7
+ ,137.5
+ ,146.1
+ ,133.6
+ ,167.9
+ ,181.9
+ ,202
+ ,166.5
+ ,151.3
+ ,146.2
+ ,148.3
+ ,144.7
+ ,123.6
+ ,151.6
+ ,133.9
+ ,137.4
+ ,181.6
+ ,182
+ ,190
+ ,161.2
+ ,155.5
+ ,141.9
+ ,164.6
+ ,136.2
+ ,126.8
+ ,152.5
+ ,126.6
+ ,150.1
+ ,186.3
+ ,147.5
+ ,200.4
+ ,177.2
+ ,127.4
+ ,177.1
+ ,154.4
+ ,135.2
+ ,126.4
+ ,147.3
+ ,140.6
+ ,152.3
+ ,151.2
+ ,172.2
+ ,215.3
+ ,154.1
+ ,159.3
+ ,160.4
+ ,151.9
+ ,148.4
+ ,139.6
+ ,148.2
+ ,153.5
+ ,145.1
+ ,183.7
+ ,210.5
+ ,203.3
+ ,153.3
+ ,144.3
+ ,169.6
+ ,143.7
+ ,160.1
+ ,135.6
+ ,141.8
+ ,159.9
+ ,145.7
+ ,183.5
+ ,198.2
+ ,186.8
+ ,172
+ ,150.6
+ ,163.3
+ ,153.7
+ ,152.9
+ ,135.5
+ ,148.5
+ ,148.4
+ ,133.6
+ ,194.1
+ ,208.6
+ ,197.3
+ ,164.4
+ ,148.1
+ ,152
+ ,144.1
+ ,155
+ ,124.5
+ ,153
+ ,146
+ ,138
+ ,190
+ ,192
+ ,192
+ ,147
+ ,133
+ ,163
+ ,150
+ ,129
+ ,131
+ ,145
+ ,137
+ ,138
+ ,168
+ ,176
+ ,188
+ ,139
+ ,143
+ ,150
+ ,154
+ ,137
+ ,129
+ ,128
+ ,140
+ ,143
+ ,151
+ ,177
+ ,184
+ ,151
+ ,134
+ ,164
+ ,126
+ ,131
+ ,125
+ ,127
+ ,143
+ ,143
+ ,160
+ ,190
+ ,182
+ ,138
+ ,136
+ ,152
+ ,127
+ ,151
+ ,130
+ ,119
+ ,153)
> 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] 476
> (np <- floor(n / par1))
[1] 39
> 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] 40 40 40 40 40 40 40 40 39 39 39 39
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]  93.2  95.9  96.3  89.8  97.1 105.1 100.0 106.0 106.1 106.9 104.6 116.4
 [2,]  96.0  82.8  84.4  84.4  93.0  92.5  84.8  91.2 102.1  96.6 107.1 112.6
 [3,]  95.2  83.3  91.2  87.2  96.0  97.1  94.3  96.6 105.2 127.3 123.5 123.8
 [4,]  77.1  80.0  81.9  85.6  80.5  81.4  87.1  96.3 101.0  98.2  98.8 103.6
 [5,]  70.9  80.4  80.5  72.0  76.1  79.1  90.3  88.2  84.3 100.2  98.6 113.9
 [6,]  64.8  67.5  70.4  69.2  69.9  72.1  72.4  70.2  87.5  89.4  90.6  98.6
 [7,]  70.1  75.7  74.8  77.5  73.6  78.7  84.9  86.5  92.7  95.3  89.1  95.0
 [8,]  77.3  71.1  75.9  78.1  92.6  87.1  92.7  88.2  94.4 104.2 105.2 116.0
 [9,]  79.5  89.3  86.3  94.3  94.2  91.4  92.2 102.8 113.0 106.4 114.0 113.9
[10,] 100.6 101.1  98.7  97.7  93.5 109.9 114.9 119.1 113.9 116.2 122.1 127.5
[11,] 100.7 105.2 100.9 100.2 108.5 116.3 112.5 119.2 122.9 135.9 138.0 131.4
[12,] 107.1 114.1 113.8 116.4 109.4 113.0 118.3 125.1 132.7 134.0 142.2 145.9
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] 131.5 131.2 126.9 132.6 128.1 153.1 164.1 161.4 182.7 159.8 170.7 173.7
 [2,] 131.0 125.4 124.0 123.7 138.4 147.6 142.8 169.4 145.2 157.8 164.4 153.8
 [3,] 130.5 126.5 135.7 153.3 152.7 157.7 157.9 168.8 182.1 168.2 179.7 175.6
 [4,] 118.9 119.4 130.0 134.0 120.0 137.2 159.2 158.1 158.7 158.4 157.0 147.1
 [5,] 114.3 113.5 109.4 119.6 140.5 151.5 162.2 158.5 141.6 152.0 168.0 160.3
 [6,]  85.7  98.7 117.8 116.2 116.2  98.7 123.1 135.3 132.6 142.2 139.3 135.2
 [7,] 104.6 114.5 120.3 118.6 121.4 145.8 130.0 149.3 139.6 137.2 138.6 148.8
 [8,] 105.1 113.8 121.0 130.7 127.8 151.7 150.1 143.4 147.0 152.6 153.4 151.0
 [9,] 117.3 133.1 132.3 129.3 143.6 129.4 169.4 142.2 166.6 166.8 138.9 148.2
[10,] 142.5 143.4 142.9 144.4 157.6 174.1 179.7 188.4 157.0 165.6 172.1 182.2
[11,] 140.0 137.3 147.4 163.2 166.2 197.0 182.1 166.2 180.4 198.6 198.4 189.2
[12,] 159.8 165.2 175.9 179.4 182.3 193.9 194.3 199.2 210.2 201.5 217.8 183.1
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
 [1,] 170.0 203.6 159.4 141.6 151.5 166.5 161.2 177.2 154.1 153.3 172.0 164.4
 [2,] 158.4 175.2 161.1 148.1 158.1 151.3 155.5 127.4 159.3 144.3 150.6 148.1
 [3,] 176.1 168.7 172.1 184.9 143.0 146.2 141.9 177.1 160.4 169.6 163.3 152.0
 [4,] 156.2 155.9 158.4 132.5 151.2 148.3 164.6 154.4 151.9 143.7 153.7 144.1
 [5,] 153.2 147.3 114.6 137.3 147.6 144.7 136.2 135.2 148.4 160.1 152.9 155.0
 [6,] 117.9 137.0 159.6 135.5 130.7 123.6 126.8 126.4 139.6 135.6 135.5 124.5
 [7,] 149.8 141.1 159.7 121.7 137.5 151.6 152.5 147.3 148.2 141.8 148.5 153.0
 [8,] 156.6 167.4 159.4 166.1 146.1 133.9 126.6 140.6 153.5 159.9 148.4 146.0
 [9,] 166.7 160.2 160.7 146.8 133.6 137.4 150.1 152.3 145.1 145.7 133.6 138.0
[10,] 156.8 191.9 165.5 162.8 167.9 181.6 186.3 151.2 183.7 183.5 194.1 190.0
[11,] 158.6 174.4 205.0 186.8 181.9 182.0 147.5 172.2 210.5 198.2 208.6 192.0
[12,] 210.8 208.2 205.2 185.5 202.0 190.0 200.4 215.3 203.3 186.8 197.3 192.0
      [,37] [,38] [,39] [,40]
 [1,]   147   139   151   138
 [2,]   133   143   134   136
 [3,]   163   150   164   152
 [4,]   150   154   126   127
 [5,]   129   137   131   151
 [6,]   131   129   125   130
 [7,]   145   128   127   119
 [8,]   137   140   143   153
 [9,]   138   143   143    NA
[10,]   168   151   160    NA
[11,]   176   177   190    NA
[12,]   188   184   182    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]   2.8 -13.1 -11.9  -5.4  -4.1 -12.6 -15.2 -14.8  -4.0 -10.3   2.5  -3.8
 [2,]  -0.8   0.5   6.8   2.8   3.0   4.6   9.5   5.4   3.1  30.7  16.4  11.2
 [3,] -18.1  -3.3  -9.3  -1.6 -15.5 -15.7  -7.2  -0.3  -4.2 -29.1 -24.7 -20.2
 [4,]  -6.2   0.4  -1.4 -13.6  -4.4  -2.3   3.2  -8.1 -16.7   2.0  -0.2  10.3
 [5,]  -6.1 -12.9 -10.1  -2.8  -6.2  -7.0 -17.9 -18.0   3.2 -10.8  -8.0 -15.3
 [6,]   5.3   8.2   4.4   8.3   3.7   6.6  12.5  16.3   5.2   5.9  -1.5  -3.6
 [7,]   7.2  -4.6   1.1   0.6  19.0   8.4   7.8   1.7   1.7   8.9  16.1  21.0
 [8,]   2.2  18.2  10.4  16.2   1.6   4.3  -0.5  14.6  18.6   2.2   8.8  -2.1
 [9,]  21.1  11.8  12.4   3.4  -0.7  18.5  22.7  16.3   0.9   9.8   8.1  13.6
[10,]   0.1   4.1   2.2   2.5  15.0   6.4  -2.4   0.1   9.0  19.7  15.9   3.9
[11,]   6.4   8.9  12.9  16.2   0.9  -3.3   5.8   5.9   9.8  -1.9   4.2  14.5
[12,] -11.2 -17.8 -24.0 -19.3  -4.3 -13.0 -12.3 -19.0 -25.8 -29.4 -25.8 -14.4
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,]  -0.5  -5.8  -2.9  -8.9  10.3  -5.5 -21.3   8.0 -37.5  -2.0  -6.3 -19.9
 [2,]  -0.5   1.1  11.7  29.6  14.3  10.1  15.1  -0.6  36.9  10.4  15.3  21.8
 [3,] -11.6  -7.1  -5.7 -19.3 -32.7 -20.5   1.3 -10.7 -23.4  -9.8 -22.7 -28.5
 [4,]  -4.6  -5.9 -20.6 -14.4  20.5  14.3   3.0   0.4 -17.1  -6.4  11.0  13.2
 [5,] -28.6 -14.8   8.4  -3.4 -24.3 -52.8 -39.1 -23.2  -9.0  -9.8 -28.7 -25.1
 [6,]  18.9  15.8   2.5   2.4   5.2  47.1   6.9  14.0   7.0  -5.0  -0.7  13.6
 [7,]   0.5  -0.7   0.7  12.1   6.4   5.9  20.1  -5.9   7.4  15.4  14.8   2.2
 [8,]  12.2  19.3  11.3  -1.4  15.8 -22.3  19.3  -1.2  19.6  14.2 -14.5  -2.8
 [9,]  25.2  10.3  10.6  15.1  14.0  44.7  10.3  46.2  -9.6  -1.2  33.2  34.0
[10,]  -2.5  -6.1   4.5  18.8   8.6  22.9   2.4 -22.2  23.4  33.0  26.3   7.0
[11,]  19.8  27.9  28.5  16.2  16.1  -3.1  12.2  33.0  29.8   2.9  19.4  -6.1
[12,] -28.6 -38.3 -43.3 -51.3 -29.2 -29.8 -32.9 -16.5 -50.4 -30.8 -44.1 -13.1
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
 [1,] -11.6 -28.4   1.7   6.5   6.6 -15.2  -5.7 -49.8   5.2  -9.0 -21.4 -16.3
 [2,]  17.7  -6.5  11.0  36.8 -15.1  -5.1 -13.6  49.7   1.1  25.3  12.7   3.9
 [3,] -19.9 -12.8 -13.7 -52.4   8.2   2.1  22.7 -22.7  -8.5 -25.9  -9.6  -7.9
 [4,]  -3.0  -8.6 -43.8   4.8  -3.6  -3.6 -28.4 -19.2  -3.5  16.4  -0.8  10.9
 [5,] -35.3 -10.3  45.0  -1.8 -16.9 -21.1  -9.4  -8.8  -8.8 -24.5 -17.4 -30.5
 [6,]  31.9   4.1   0.1 -13.8   6.8  28.0  25.7  20.9   8.6   6.2  13.0  28.5
 [7,]   6.8  26.3  -0.3  44.4   8.6 -17.7 -25.9  -6.7   5.3  18.1  -0.1  -7.0
 [8,]  10.1  -7.2   1.3 -19.3 -12.5   3.5  23.5  11.7  -8.4 -14.2 -14.8  -8.0
 [9,]  -9.9  31.7   4.8  16.0  34.3  44.2  36.2  -1.1  38.6  37.8  60.5  52.0
[10,]   1.8 -17.5  39.5  24.0  14.0   0.4 -38.8  21.0  26.8  14.7  14.5   2.0
[11,]  52.2  33.8   0.2  -1.3  20.1   8.0  52.9  43.1  -7.2 -11.4 -11.3   0.0
[12,]  -7.2 -48.8 -63.6 -34.0 -35.5 -28.8 -23.2 -61.2 -50.0 -14.8 -32.9 -45.0
      [,37] [,38] [,39] [,40]
 [1,]   -14     4   -17    -2
 [2,]    30     7    30    16
 [3,]   -13     4   -38   -25
 [4,]   -21   -17     5    24
 [5,]     2    -8    -6   -21
 [6,]    14    -1     2   -11
 [7,]    -8    12    16    34
 [8,]     1     3     0    NA
 [9,]    30     8    17    NA
[10,]     8    26    30    NA
[11,]    12     7    -8    NA
[12,]   -49   -33   -44    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/1y27y1272274987.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/rcomp/tmp/2y27y1272274987.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/rcomp/tmp/3y27y1272274987.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/rcomp/tmp/4rb601272274987.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] [,10]  [,11]
[1,]  89.80  82.80  83.30  77.1  70.90  64.80  70.10  71.10  79.50  93.5 100.20
[2,] 106.50 104.60 123.65  99.9  99.40  88.45  93.85 104.65 113.45 120.6 133.65
[3,] 144.30 137.20 152.00 135.6 135.70 123.35 127.50 138.50 137.40 157.0 166.20
[4,] 162.75 152.55 168.45 154.2 151.25 133.90 146.55 151.35 146.25 176.9 188.00
[5,] 203.60 175.20 184.90 164.6 168.00 159.60 159.70 167.40 169.40 194.1 210.50
     [,12]
[1,] 107.1
[2,] 138.1
[3,] 184.0
[4,] 199.8
[5,] 217.8

$n
 [1] 40 40 40 40 40 40 40 40 39 39 39 39

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 130.2476 125.2211 140.8081 122.0348 122.7468 111.9957 114.3345 126.8334
[2,] 158.3524 149.1789 163.1919 149.1652 148.6532 134.7043 140.6655 150.1666
         [,9]    [,10]    [,11]    [,12]
[1,] 129.1015 142.7560 152.4493 168.3897
[2,] 145.6985 171.2440 179.9507 199.6103

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(89.8, 106.5, 144.3, 162.75, 203.6, 82.8, 104.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5rb601272274987.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]  [,10] [,11]
[1,] -28.40 -15.10 -38.0 -28.40 -39.10 -13.80 -17.7 -22.30 -9.90 -22.20 -11.4
[2,] -15.00   1.95 -22.7 -11.10 -22.15   2.45   0.2  -2.45  8.95   1.90   0.1
[3,]  -6.05  10.25 -12.9  -3.25 -10.55   6.70   6.6   2.20 16.00   8.00   8.9
[4,]  -1.25  17.05  -6.4   4.00  -6.60  14.00  15.1  13.20 33.60  20.35  19.6
[5,]  10.30  36.90   8.2  24.00   8.40  28.50  34.0  23.50 60.50  39.50  43.1
      [,12]
[1,] -63.60
[2,] -43.65
[3,] -29.40
[4,] -18.40
[5,]  -4.30

$n
 [1] 40 40 40 40 40 40 40 39 39 39 39 39

$conf
          [,1]      [,2]       [,3]      [,4]     [,5]    [,6]      [,7]
[1,] -9.485024  6.477719 -16.972065 -7.022281 -14.4347 3.81458  2.877683
[2,] -2.614976 14.022281  -8.827935  0.522281  -6.6653 9.58542 10.322317
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -1.759489  9.763489  3.332104  3.966452 -35.78831
[2,]  6.159489 22.236511 12.667896 13.833548 -23.01169

$out
 [1] -37.5 -49.8  49.7 -52.4  22.7 -43.8 -52.8  45.0  47.1  31.9  44.4 -25.9
[13] -38.8  52.2  52.9

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

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

> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/623n41272274987.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]  [,10]
[1,]  64.80  67.50  70.40  69.2  69.90  72.10  72.40  70.20  84.30  89.40
[2,]  74.00  77.85  78.20  77.8  78.30  80.25  86.00  88.20  93.55  97.40
[3,]  86.35  83.05  85.35  86.4  93.25  91.95  92.45  96.45 103.65 105.30
[4,]  98.30  98.50  97.50  96.0  96.55 107.50 106.25 112.55 113.45 121.75
[5,] 107.10 114.10 113.80 116.4 109.40 116.30 118.30 125.10 132.70 135.90
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17] [,18]  [,19]  [,20]
[1,]  89.10  95.00  85.70  98.70 109.40 116.20 116.20 129.4 123.10 135.30
[2,]  98.70 108.10 109.70 114.15 120.65 121.65 124.60 141.5 146.45 146.35
[3,] 106.15 114.95 124.70 125.95 128.45 131.65 139.45 151.6 160.70 159.95
[4,] 122.80 125.65 135.75 135.20 139.30 148.85 155.15 165.9 174.55 169.10
[5,] 142.20 145.90 159.80 165.20 147.40 179.40 182.30 197.0 194.30 199.20
      [,21] [,22]  [,23]  [,24]  [,25]  [,26] [,27]  [,28]  [,29]  [,30]  [,31]
[1,] 132.60 137.2 138.60 135.20 149.80 137.00 158.4 121.70 130.70 123.60 126.60
[2,] 143.40 152.3 146.35 148.50 154.70 151.60 159.4 136.40 140.25 141.05 139.05
[3,] 157.85 159.1 166.20 157.05 157.60 168.05 160.2 147.45 149.40 149.80 151.30
[4,] 181.25 167.5 175.90 178.90 168.35 183.55 168.8 175.50 163.00 174.05 162.90
[5,] 210.20 168.2 217.80 189.20 176.10 208.20 172.1 186.80 181.90 190.00 186.30
      [,32]  [,33]  [,34]  [,35]  [,36] [,37] [,38] [,39] [,40]
[1,] 126.40 139.60 135.60 133.60 124.50 129.0 128.0   125 119.0
[2,] 137.90 148.30 144.00 148.45 145.05 135.0 138.0   129 128.5
[3,] 151.75 153.80 156.60 153.30 152.50 146.0 143.0   143 137.0
[4,] 174.65 172.05 176.55 183.05 177.20 165.5 152.5   162 151.5
[5,] 215.30 203.30 198.20 208.60 192.00 188.0 154.0   190 153.0

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

$conf
        [,1]    [,2]     [,3]     [,4]      [,5]     [,6]      [,7]     [,8]
[1,] 75.2666 73.6314 76.54714 78.09886  84.92605  79.5211  83.21384  85.3438
[2,] 97.4334 92.4686 94.15286 94.70114 101.57395 104.3789 101.68616 107.5562
          [,9]    [,10]     [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,]  94.57348  94.1938  95.15783 106.9453 112.8184 116.3490 119.9436 119.2439
[2,] 112.72652 116.4062 117.14217 122.9547 136.5816 135.5510 136.9564 144.0561
        [,17]   [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 125.5159 140.471 147.8834 149.5736 140.5864 152.1672 152.7220 143.1844
[2,] 153.3841 162.729 173.5166 170.3264 175.1136 166.0328 179.6780 170.9156
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30]    [,31]    [,32]
[1,] 151.3741 153.4774 155.9126 129.6162 139.0236 134.7485 140.4219 134.9881
[2,] 163.8259 182.6226 164.4874 165.2838 159.7764 164.8515 162.1781 168.5119
        [,33]    [,34]    [,35]    [,36]    [,37]    [,38]    [,39]    [,40]
[1,] 142.9675 141.7537 137.5187 137.8362 132.0887 136.3865 127.9485 124.1519
[2,] 164.6325 171.4463 169.0813 167.1638 159.9113 149.6135 158.0515 149.8481

$out
 [1] 175.9  98.7 198.6 201.5 117.9 210.8 114.6 205.0 205.2 202.0 200.4 210.5
[13] 177.0 184.0

$group
 [1] 15 18 22 22 25 25 27 27 27 29 31 33 38 38

$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" "38" "39" NA  

Warning message:
In bxp(list(stats = c(64.8, 74, 86.35, 98.3, 107.1, 67.5, 77.85,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/723n41272274987.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,] 111.7825 123.35 111.0875
[2,] 126.4288 135.65 126.3688
[3,] 130.6685 137.95 129.3687
[4,] 146.2918 154.50 147.3875
[5,] 171.4205 166.20 168.9500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 121.6088 129.3524 119.7820
[2,] 139.7281 146.5476 138.9555

$out
[1] 184

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(111.7825, 126.42875, 130.668461538462, 146.291762820513,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1y27y1272274987.ps tmp/1y27y1272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/2y27y1272274987.ps tmp/2y27y1272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/3y27y1272274987.ps tmp/3y27y1272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/4rb601272274987.ps tmp/4rb601272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/5rb601272274987.ps tmp/5rb601272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/623n41272274987.ps tmp/623n41272274987.png",intern=TRUE))
character(0)
> try(system("convert tmp/723n41272274987.ps tmp/723n41272274987.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
   1.68    1.48    2.72