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(31.82
+ ,32.23
+ ,30.51
+ ,30.68
+ ,31.32
+ ,30.77
+ ,30.88
+ ,31.92
+ ,29.68
+ ,28.56
+ ,32.85
+ ,36.35
+ ,35.13
+ ,30.05
+ ,30.80
+ ,32.00
+ ,31.78
+ ,29.92
+ ,31.77
+ ,29.93
+ ,31.05
+ ,32.25
+ ,33.34
+ ,35.35
+ ,33.45
+ ,32.48
+ ,31.55
+ ,32.02
+ ,29.49
+ ,27.18
+ ,28.66
+ ,30.05
+ ,31.14
+ ,29.66
+ ,30.05
+ ,26.32
+ ,25.32
+ ,23.06
+ ,24.45
+ ,24.91
+ ,26.65
+ ,28.20
+ ,30.50
+ ,26.90
+ ,28.00
+ ,26.12
+ ,29.80
+ ,32.23
+ ,31.64
+ ,29.66
+ ,31.30
+ ,33.42
+ ,35.65
+ ,35.87
+ ,34.73
+ ,34.40
+ ,33.56
+ ,34.36
+ ,33.45
+ ,32.98
+ ,32.72
+ ,32.75
+ ,34.21
+ ,36.69
+ ,37.50
+ ,37.52
+ ,37.23
+ ,35.72
+ ,35.61
+ ,39.34
+ ,39.31
+ ,40.62
+ ,40.62
+ ,41.80
+ ,42.60
+ ,43.67
+ ,45.47
+ ,46.40
+ ,45.84
+ ,46.59
+ ,46.70
+ ,48.81
+ ,48.92
+ ,46.95
+ ,44.90
+ ,43.63
+ ,45.28
+ ,46.05
+ ,46.51
+ ,46.95
+ ,46.45
+ ,49.82
+ ,51.82
+ ,52.50
+ ,53.25
+ ,53.95
+ ,54.54
+ ,53.23
+ ,52.00
+ ,51.13
+ ,50.16
+ ,50.19
+ ,51.28
+ ,53.34
+ ,53.08
+ ,51.75
+ ,52.20
+ ,51.11
+ ,51.30
+ ,52.61
+ ,51.41
+ ,51.61
+ ,52.35
+ ,52.21
+ ,52.64
+ ,52.50
+ ,53.55
+ ,52.41
+ ,53.80
+ ,54.22
+ ,53.21
+ ,53.50
+ ,52.30
+ ,52.47
+ ,51.32
+ ,51.37
+ ,52.32
+ ,52.24
+ ,53.40
+ ,52.83
+ ,53.08
+ ,53.62
+ ,54.14
+ ,54.62
+ ,54.13
+ ,52.94
+ ,54.43
+ ,53.98
+ ,54.05
+ ,53.45
+ ,53.39
+ ,53.43
+ ,53.03
+ ,52.16
+ ,51.72
+ ,51.42
+ ,50.54
+ ,50.13
+ ,50.63
+ ,52.10
+ ,52.64
+ ,53.47
+ ,53.39
+ ,54.11
+ ,54.70
+ ,54.34
+ ,53.39
+ ,52.47
+ ,52.74
+ ,53.44
+ ,51.61
+ ,53.30
+ ,53.10
+ ,53.26
+ ,53.09
+ ,53.39
+ ,52.56
+ ,52.57
+ ,53.18
+ ,52.00
+ ,51.75
+ ,50.34
+ ,51.96
+ ,53.07
+ ,53.14
+ ,52.60
+ ,52.70
+ ,53.20
+ ,52.23
+ ,53.21
+ ,54.17
+ ,54.33
+ ,55.12
+ ,56.58
+ ,55.91
+ ,54.65
+ ,55.03
+ ,55.56
+ ,57.15
+ ,58.24
+ ,58.77
+ ,58.97
+ ,57.29
+ ,56.04
+ ,56.72
+ ,56.46
+ ,57.10
+ ,56.81
+ ,56.32
+ ,55.81
+ ,56.03
+ ,56.90
+ ,56.59
+ ,54.60
+ ,55.08
+ ,55.48
+ ,54.42
+ ,54.55
+ ,55.17
+ ,55.10
+ ,53.89
+ ,53.02
+ ,53.24
+ ,52.32
+ ,52.57
+ ,51.69
+ ,50.26
+ ,47.94
+ ,47.32
+ ,48.72
+ ,49.84
+ ,50.15
+ ,50.00
+ ,49.34)
> par1 = '50'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 224
> (np <- floor(n / par1))
[1] 4
> 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] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 4 4 4 4 4 4 4 4 4 4 4 4 4 4
[39] 4 4 4 4 4 4 4 4 4 4 4 4
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] 31.82 31.30 50.16 52.64 56.03
 [2,] 32.23 33.42 50.19 53.47 56.90
 [3,] 30.51 35.65 51.28 53.39 56.59
 [4,] 30.68 35.87 53.34 54.11 54.60
 [5,] 31.32 34.73 53.08 54.70 55.08
 [6,] 30.77 34.40 51.75 54.34 55.48
 [7,] 30.88 33.56 52.20 53.39 54.42
 [8,] 31.92 34.36 51.11 52.47 54.55
 [9,] 29.68 33.45 51.30 52.74 55.17
[10,] 28.56 32.98 52.61 53.44 55.10
[11,] 32.85 32.72 51.41 51.61 53.89
[12,] 36.35 32.75 51.61 53.30 53.02
[13,] 35.13 34.21 52.35 53.10 53.24
[14,] 30.05 36.69 52.21 53.26 52.32
[15,] 30.80 37.50 52.64 53.09 52.57
[16,] 32.00 37.52 52.50 53.39 51.69
[17,] 31.78 37.23 53.55 52.56 50.26
[18,] 29.92 35.72 52.41 52.57 47.94
[19,] 31.77 35.61 53.80 53.18 47.32
[20,] 29.93 39.34 54.22 52.00 48.72
[21,] 31.05 39.31 53.21 51.75 49.84
[22,] 32.25 40.62 53.50 50.34 50.15
[23,] 33.34 40.62 52.30 51.96 50.00
[24,] 35.35 41.80 52.47 53.07 49.34
[25,] 33.45 42.60 51.32 53.14    NA
[26,] 32.48 43.67 51.37 52.60    NA
[27,] 31.55 45.47 52.32 52.70    NA
[28,] 32.02 46.40 52.24 53.20    NA
[29,] 29.49 45.84 53.40 52.23    NA
[30,] 27.18 46.59 52.83 53.21    NA
[31,] 28.66 46.70 53.08 54.17    NA
[32,] 30.05 48.81 53.62 54.33    NA
[33,] 31.14 48.92 54.14 55.12    NA
[34,] 29.66 46.95 54.62 56.58    NA
[35,] 30.05 44.90 54.13 55.91    NA
[36,] 26.32 43.63 52.94 54.65    NA
[37,] 25.32 45.28 54.43 55.03    NA
[38,] 23.06 46.05 53.98 55.56    NA
[39,] 24.45 46.51 54.05 57.15    NA
[40,] 24.91 46.95 53.45 58.24    NA
[41,] 26.65 46.45 53.39 58.77    NA
[42,] 28.20 49.82 53.43 58.97    NA
[43,] 30.50 51.82 53.03 57.29    NA
[44,] 26.90 52.50 52.16 56.04    NA
[45,] 28.00 53.25 51.72 56.72    NA
[46,] 26.12 53.95 51.42 56.46    NA
[47,] 29.80 54.54 50.54 57.10    NA
[48,] 32.23 53.23 50.13 56.81    NA
[49,] 31.64 52.00 50.63 56.32    NA
[50,] 29.66 51.13 52.10 55.81    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,]  0.41  2.12  0.03  0.83  0.87
 [2,] -1.72  2.23  1.09 -0.08 -0.31
 [3,]  0.17  0.22  2.06  0.72 -1.99
 [4,]  0.64 -1.14 -0.26  0.59  0.48
 [5,] -0.55 -0.33 -1.33 -0.36  0.40
 [6,]  0.11 -0.84  0.45 -0.95 -1.06
 [7,]  1.04  0.80 -1.09 -0.92  0.13
 [8,] -2.24 -0.91  0.19  0.27  0.62
 [9,] -1.12 -0.47  1.31  0.70 -0.07
[10,]  4.29 -0.26 -1.20 -1.83 -1.21
[11,]  3.50  0.03  0.20  1.69 -0.87
[12,] -1.22  1.46  0.74 -0.20  0.22
[13,] -5.08  2.48 -0.14  0.16 -0.92
[14,]  0.75  0.81  0.43 -0.17  0.25
[15,]  1.20  0.02 -0.14  0.30 -0.88
[16,] -0.22 -0.29  1.05 -0.83 -1.43
[17,] -1.86 -1.51 -1.14  0.01 -2.32
[18,]  1.85 -0.11  1.39  0.61 -0.62
[19,] -1.84  3.73  0.42 -1.18  1.40
[20,]  1.12 -0.03 -1.01 -0.25  1.12
[21,]  1.20  1.31  0.29 -1.41  0.31
[22,]  1.09  0.00 -1.20  1.62 -0.15
[23,]  2.01  1.18  0.17  1.11 -0.66
[24,] -1.90  0.80 -1.15  0.07    NA
[25,] -0.97  1.07  0.05 -0.54    NA
[26,] -0.93  1.80  0.95  0.10    NA
[27,]  0.47  0.93 -0.08  0.50    NA
[28,] -2.53 -0.56  1.16 -0.97    NA
[29,] -2.31  0.75 -0.57  0.98    NA
[30,]  1.48  0.11  0.25  0.96    NA
[31,]  1.39  2.11  0.54  0.16    NA
[32,]  1.09  0.11  0.52  0.79    NA
[33,] -1.48 -1.97  0.48  1.46    NA
[34,]  0.39 -2.05 -0.49 -0.67    NA
[35,] -3.73 -1.27 -1.19 -1.26    NA
[36,] -1.00  1.65  1.49  0.38    NA
[37,] -2.26  0.77 -0.45  0.53    NA
[38,]  1.39  0.46  0.07  1.59    NA
[39,]  0.46  0.44 -0.60  1.09    NA
[40,]  1.74 -0.50 -0.06  0.53    NA
[41,]  1.55  3.37  0.04  0.20    NA
[42,]  2.30  2.00 -0.40 -1.68    NA
[43,] -3.60  0.68 -0.87 -1.25    NA
[44,]  1.10  0.75 -0.44  0.68    NA
[45,] -1.88  0.70 -0.30 -0.26    NA
[46,]  3.68  0.59 -0.88  0.64    NA
[47,]  2.43 -1.31 -0.41 -0.29    NA
[48,] -0.59 -1.23  0.50 -0.49    NA
[49,] -1.98 -0.87  1.47 -0.51    NA
[50,]  1.64 -0.97  0.54  0.22    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/16qst1321368623.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/2pm2l1321368623.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/3wkyh1321368623.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/4swm01321368623.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] [,10] [,11] [,12]
[1,] 31.30 32.23 30.51 30.68 31.32 30.77 30.88 31.92 29.68 28.56 32.72 32.75
[2,] 31.82 33.42 35.65 35.87 34.73 34.40 33.56 34.36 33.45 32.98 32.85 36.35
[3,] 50.16 50.19 51.28 53.34 53.08 51.75 52.20 51.11 51.30 52.61 51.41 51.61
[4,] 52.64 53.47 53.39 54.11 54.70 54.34 53.39 52.47 52.74 53.44 51.61 53.02
[5,] 56.03 56.90 56.59 54.60 55.08 55.48 54.42 54.55 55.17 55.10 53.89 53.30
     [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
[1,] 34.21 30.05 30.80 32.00 31.78 29.92 31.77 29.93 31.05 32.25 33.34 35.35
[2,] 35.13 36.69 37.50 37.52 37.23 35.72 35.61 39.34 39.31 40.62 40.62 41.80
[3,] 52.35 52.21 52.57 51.69 50.26 47.94 47.32 48.72 49.84 50.15 50.00 49.34
[4,] 53.10 52.32 52.64 52.50 52.56 52.41 53.18 52.00 51.75 50.34 51.96 52.47
[5,] 53.24 53.26 53.09 53.39 53.55 52.57 53.80 54.22 53.21 53.50 52.30 53.07
      [,25]  [,26]  [,27] [,28]  [,29]  [,30]  [,31]  [,32] [,33]  [,34]  [,35]
[1,] 33.450 32.480 31.550 32.02 29.490 27.180 28.660 30.050 31.14 29.660 30.050
[2,] 38.025 38.075 38.510 39.21 37.665 36.885 37.680 39.430 40.03 38.305 37.475
[3,] 46.960 47.520 48.895 49.32 49.035 49.710 49.890 51.215 51.53 50.785 49.515
[4,] 52.230 51.985 52.510 52.72 52.815 53.020 53.625 53.975 54.63 55.600 55.020
[5,] 53.140 52.600 52.700 53.20 53.400 53.210 54.170 54.330 55.12 56.580 55.910
      [,36]  [,37]  [,38] [,39]  [,40] [,41]  [,42]  [,43] [,44]  [,45]  [,46]
[1,] 26.320 25.320 23.060 24.45 24.910 26.65 28.200 30.500 26.90 28.000 26.120
[2,] 34.975 35.300 34.555 35.48 35.930 36.55 39.010 41.160 39.53 39.860 38.770
[3,] 48.285 49.855 50.015 50.28 50.200 49.92 51.625 52.425 52.33 52.485 52.685
[4,] 53.795 54.730 54.770 55.60 55.845 56.08 56.200 55.160 54.27 54.985 55.205
[5,] 54.650 55.030 55.560 57.15 58.240 58.77 58.970 57.290 56.04 56.720 56.460
     [,47] [,48]  [,49]  [,50]
[1,] 29.80 32.23 31.640 29.660
[2,] 40.17 41.18 41.135 40.395
[3,] 52.54 51.68 51.315 51.615
[4,] 55.82 55.02 54.160 53.955
[5,] 57.10 56.81 56.320 55.810

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 4 4 4 4 4 4 4 4 4 4 4 4 4 4
[39] 4 4 4 4 4 4 4 4 4 4 4 4

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 35.44864 36.02272 38.74496 40.45166 38.96925 37.66045 38.18817 38.31352
[2,] 64.87136 64.35728 63.81504 66.22834 67.19075 65.83955 66.21183 63.90648
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 37.66973 38.15302 38.15423 39.83102 39.65244 41.16588 41.87211 41.10517
[2,] 64.93027 67.06698 64.66577 63.38898 65.04756 63.25412 63.26789 62.27483
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 39.42786 36.14689 34.90508 39.77448 41.04993 43.28187 41.98718 41.80060
[2,] 61.09214 59.73311 59.73492 57.66552 58.63007 57.01813 58.01282 56.87940
        [,25]   [,26]  [,27]   [,28]   [,29]    [,30]    [,31]    [,32]  [,33]
[1,] 35.73805 36.5311 37.835 38.6471 37.0665 36.96335 37.29345 39.72445 39.996
[2,] 58.18195 58.5089 59.955 59.9929 61.0035 62.45665 62.48655 62.70555 63.064
        [,34]    [,35]   [,36]   [,37]    [,38]   [,39]    [,40]   [,41]
[1,] 37.12195 35.65445 33.4172 34.5053 34.04515 34.3852 34.46715 34.4913
[2,] 64.44805 63.37555 63.1528 65.2047 65.98485 66.1748 65.93285 65.3487
       [,42]  [,43]   [,44]    [,45]    [,46]   [,47]   [,48]    [,49]   [,50]
[1,] 38.0449 41.365 40.6854 40.53625 39.70135 40.1765 40.7464 41.02525 40.9026
[2,] 65.2051 63.485 63.9746 64.43375 65.66865 64.9035 62.6136 61.60475 62.3274

$out
numeric(0)

$group
numeric(0)

$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" "40" "41" "42" "43" "44" "45"
[46] "46" "47" "48" "49" "50"

Warning message:
In bxp(list(stats = c(31.3, 31.82, 50.16, 52.64, 56.03, 32.23, 33.42,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5pzh51321368623.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] [,11] [,12]
[1,] 0.03 -1.72 0.17 -1.14 -0.55 -1.06 -1.09 -2.24 -1.12 -1.83 -0.87 -1.22
[2,] 0.41 -0.31 0.17 -0.26 -0.55 -0.95 -0.92 -0.91 -0.47 -1.21  0.03 -0.20
[3,] 0.83 -0.08 0.22  0.48 -0.36 -0.84  0.13  0.19 -0.07 -1.20  0.20  0.22
[4,] 0.87  1.09 0.72  0.59 -0.33  0.11  0.80  0.27  0.70 -0.26  1.69  0.74
[5,] 0.87  2.23 0.72  0.64 -0.33  0.45  1.04  0.62  1.31 -0.26  3.50  1.46
     [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23]  [,24]
[1,] -0.92 -0.17 -0.14 -1.43 -2.32 -0.62 -1.84 -1.01  0.29 -1.20 -0.66 -1.900
[2,] -0.92  0.25 -0.14 -0.83 -1.86 -0.11 -1.18 -0.25  0.29 -0.15  0.17 -1.525
[3,] -0.14  0.43  0.02 -0.29 -1.51  0.61  0.42 -0.03  0.31  0.00  1.11 -0.540
[4,]  0.16  0.75  0.30 -0.22 -1.14  1.39  1.40  1.12  1.20  1.09  1.18  0.435
[5,]  0.16  0.81  0.30 -0.22 -1.14  1.85  3.73  1.12  1.31  1.62  2.01  0.800
      [,25]  [,26]  [,27]  [,28]  [,29] [,30] [,31] [,32]  [,33] [,34]  [,35]
[1,] -0.970 -0.930 -0.080 -2.530 -2.310 0.110 0.160 0.110 -1.970 -2.05 -3.730
[2,] -0.755 -0.415  0.195 -1.750 -1.440 0.180 0.350 0.315 -1.725 -1.36 -2.500
[3,] -0.245  0.525  0.485 -0.765  0.090 0.605 0.965 0.655 -0.500 -0.58 -1.265
[4,]  0.560  1.375  0.715  0.300  0.865 1.220 1.750 0.940  0.970 -0.05 -1.225
[5,]  1.070  1.800  0.930  1.160  0.980 1.480 2.110 1.090  1.460  0.39 -1.190
      [,36]  [,37] [,38]  [,39]  [,40] [,41] [,42]  [,43]  [,44] [,45]  [,46]
[1,] -1.000 -2.260 0.070 -0.600 -0.500 0.040 -1.68 -3.600 -0.440 -1.88 -0.880
[2,] -0.310 -1.355 0.265 -0.080 -0.280 0.120 -1.04 -2.425  0.120 -1.09 -0.145
[3,]  0.935  0.040 0.925  0.450  0.235 0.875  0.80 -1.060  0.715 -0.28  0.615
[4,]  1.570  0.650 1.490  0.775  1.135 2.460  2.15 -0.095  0.925  0.22  2.160
[5,]  1.650  0.770 1.590  1.090  1.740 3.370  2.30  0.680  1.100  0.70  3.680
     [,47]  [,48]  [,49]  [,50]
[1,] -1.31 -1.230 -1.980 -0.970
[2,] -0.86 -0.910 -1.425 -0.375
[3,] -0.35 -0.540 -0.690  0.380
[4,]  1.07  0.005  0.480  1.090
[5,]  2.43  0.500  1.470  1.640

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4
[39] 4 4 4 4 4 4 4 4 4 4 4 4

$conf
          [,1]       [,2]       [,3]       [,4]       [,5]        [,6]
[1,] 0.5049652 -1.0692365 -0.1686286 -0.1206079 -0.5154514 -1.58899333
[2,] 1.1550348  0.9092365  0.6086286  1.0806079 -0.2045486 -0.09100667
          [,7]      [,8]      [,9]      [,10]      [,11]      [,12]      [,13]
[1,] -1.085348 -0.643785 -0.896719 -1.8712676 -0.9729518 -0.4442016 -0.9031253
[2,]  1.345348  1.023785  0.756719 -0.5287324  1.3729518  0.8842016  0.6231253
          [,14]      [,15]      [,16]     [,17]      [,18]     [,19]      [,20]
[1,] 0.07670126 -0.2909029 -0.7210245 -2.018750 -0.4498962 -1.403022 -0.9980385
[2,] 0.78329874  0.3309029  0.1410245 -1.001250  1.6698962  2.243022  0.9380385
          [,21]      [,22]     [,23]   [,24]    [,25]   [,26]  [,27]   [,28]
[1,] -0.3330037 -0.8761809 0.3963365 -2.0884 -1.28385 -0.8891 0.0742 -2.3845
[2,]  0.9530037  0.8761809 1.8236635  1.0084  0.79385  1.9391 0.8958  0.8545
        [,29]   [,30]  [,31]   [,32]    [,33]   [,34]    [,35]   [,36]    [,37]
[1,] -1.73095 -0.2166 -0.141 0.16125 -2.62905 -1.6149 -2.27225 -0.5502 -1.54395
[2,]  1.91095  1.4266  2.071 1.14875  1.62905  0.4549 -0.25775  2.4202  1.62395
        [,38]    [,39]    [,40]   [,41]   [,42]   [,43]   [,44]   [,45]
[1,] -0.04275 -0.22545 -0.88285 -0.9736 -1.7201 -2.9007 0.07905 -1.3149
[2,]  1.89275  1.12545  1.35285  2.7236  3.3201  0.7807 1.35095  0.7549
        [,46]   [,47]    [,48]    [,49]    [,50]
[1,] -1.20595 -1.8747 -1.26285 -2.19495 -0.77735
[2,]  2.43595  1.1747  0.18285  0.81495  1.53735

$out
 [1]  2.12  2.06 -1.99 -1.33  0.40  4.29 -5.08  2.48  1.20 -0.88  1.05  0.01
[13] -1.41

$group
 [1]  1  3  3  5  5 10 13 13 15 15 16 17 21

$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" "40" "41" "42" "43" "44" "45"
[46] "46" "47" "48" "49" "50"

Warning message:
In bxp(list(stats = c(0.0300000000000011, 0.409999999999997, 0.829999999999998,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6x15v1321368623.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]
[1,] 24.450 31.300 50.130 50.340 47.320
[2,] 28.660 35.650 51.610 52.740 50.075
[3,] 30.595 43.115 52.485 53.415 53.130
[4,] 31.920 46.950 53.430 55.910 55.090
[5,] 36.350 54.540 54.620 58.970 56.900

$n
[1] 50 50 50 50 24

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]
[1,] 29.86657 40.59006 52.07833 52.70668 51.51258
[2,] 31.32343 45.63994 52.89167 54.12332 54.74742

$out
[1] 23.06

$group
[1] 1

$names
[1] "1" "2" "3" "4" NA 

Warning message:
In bxp(list(stats = c(24.45, 28.66, 30.595, 31.92, 36.35, 31.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7id171321368623.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,] 43.7120 47.3200 42.2300
[2,] 44.9525 49.8400 44.5200
[3,] 45.4520 50.9475 46.3125
[4,] 46.4060 51.7500 48.1175
[5,] 48.1600 53.3400 50.2925

$n
[1] 50 50 50

$conf
         [,1]     [,2]     [,3]
[1,] 45.12722 50.52072 45.50865
[2,] 45.77678 51.37428 47.11635

$out
[1] 46.96

$group
[1] 2

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

> dev.off()
null device 
          1 
> 
> try(system("convert tmp/16qst1321368623.ps tmp/16qst1321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/2pm2l1321368623.ps tmp/2pm2l1321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/3wkyh1321368623.ps tmp/3wkyh1321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/4swm01321368623.ps tmp/4swm01321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pzh51321368623.ps tmp/5pzh51321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/6x15v1321368623.ps tmp/6x15v1321368623.png",intern=TRUE))
character(0)
> try(system("convert tmp/7id171321368623.ps tmp/7id171321368623.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.680   0.220   2.908