R version 2.13.0 (2011-04-13)
Copyright (C) 2011 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(61.04
+ ,64.21
+ ,63.92
+ ,68.48
+ ,67.76
+ ,65.37
+ ,61.84
+ ,59.88
+ ,62.91
+ ,63.67
+ ,59.96
+ ,62.80
+ ,64.26
+ ,66.35
+ ,67.43
+ ,69.45
+ ,72.98
+ ,71.58
+ ,69.97
+ ,71.93
+ ,68.53
+ ,71.29
+ ,72.75
+ ,71.64
+ ,69.97
+ ,70.78
+ ,73.52
+ ,73.86
+ ,76.80
+ ,74.57
+ ,73.36
+ ,74.77
+ ,74.35
+ ,71.14
+ ,72.51
+ ,69.19
+ ,66.25
+ ,63.33
+ ,60.55
+ ,62.91
+ ,59.76
+ ,58.57
+ ,59.33
+ ,60.75
+ ,59.05
+ ,59.59
+ ,58.97
+ ,59.24
+ ,63.28
+ ,62.03
+ ,63.43
+ ,62.41
+ ,61.05
+ ,56.31
+ ,53.30
+ ,51.99
+ ,55.42
+ ,59.02
+ ,59.89
+ ,59.39
+ ,61.44
+ ,61.64
+ ,59.69
+ ,57.11
+ ,62.28
+ ,65.87
+ ,64.28
+ ,63.63
+ ,64.11
+ ,66.46
+ ,61.93
+ ,62.37
+ ,64.94
+ ,65.20
+ ,65.08
+ ,64.76
+ ,68.00
+ ,69.14
+ ,70.68
+ ,72.81
+ ,73.93
+ ,75.79
+ ,77.02
+ ,75.48
+ ,71.47
+ ,71.98
+ ,71.09
+ ,74.04
+ ,76.70
+ ,79.10
+ ,81.62
+ ,81.66
+ ,81.22
+ ,83.69
+ ,88.60
+ ,91.86
+ ,95.93
+ ,96.32
+ ,93.84
+ ,98.18
+ ,88.71
+ ,88.28
+ ,91.27
+ ,93.31
+ ,96.00
+ ,95.98
+ ,97.91
+ ,92.69
+ ,90.57
+ ,90.71
+ ,88.96
+ ,91.77
+ ,95.50
+ ,98.81
+ ,01.84
+ ,05.15
+ ,10.21
+ ,01.84
+ ,05.62
+ ,06.23
+ ,10.14
+ ,16.69
+ ,18.52
+ ,16.32
+ ,25.96
+ ,26.29
+ ,32.19
+ ,27.35
+ ,38.54
+ ,34.86
+ ,35.36
+ ,40.21
+ ,45.29
+ ,45.08
+ ,28.88
+ ,23.26
+ ,25.10
+ ,15.20
+ ,13.77
+ ,14.59
+ ,15.46
+ ,06.23
+ ,01.18
+ ,04.55
+ ,06.89
+ ,93.88
+ ,77.70
+ ,71.85
+ ,64.15
+ ,67.81
+ ,61.04
+ ,57.04
+ ,49.93
+ ,54.43
+ ,40.81
+ ,46.28
+ ,42.36
+ ,37.71
+ ,44.60
+ ,46.34
+ ,40.83
+ ,36.51
+ ,46.47
+ ,41.68
+ ,40.17
+ ,37.51
+ ,40.03
+ ,44.76
+ ,45.52
+ ,46.25
+ ,52.07
+ ,52.38
+ ,52.51
+ ,52.24
+ ,50.33
+ ,51.55
+ ,53.20
+ ,58.63
+ ,56.34
+ ,61.67
+ ,66.31
+ ,68.44
+ ,72.04
+ ,69.55
+ ,69.16
+ ,65.63
+ ,59.89
+ ,63.56
+ ,68.05
+ ,69.45
+ ,70.93
+ ,67.51
+ ,73.89
+ ,72.74
+ ,68.02
+ ,69.29
+ ,72.04
+ ,66.02
+ ,69.95
+ ,71.77
+ ,78.53
+ ,80.50
+ ,77.00
+ ,77.43
+ ,76.35
+ ,77.47
+ ,76.05
+ ,75.47
+ ,69.87
+ ,73.36
+ ,78.05
+ ,79.36
+ ,82.75
+ ,78.00
+ ,74.54
+ ,72.89
+ ,71.19
+ ,74.13
+ ,79.81
+ ,79.66
+ ,81.50
+ ,81.24
+ ,80.68
+ ,80.00
+ ,84.87
+ ,84.92
+ ,83.24
+ ,85.12
+ ,86.15
+ ,75.11
+ ,71.61
+ ,70.04
+ ,73.97
+ ,71.51
+ ,73.78
+ ,77.18
+ ,78.86
+ ,72.14
+ ,76.09
+ ,76.01
+ ,78.98
+ ,78.95
+ ,80.70
+ ,75.39
+ ,73.82
+ ,75.17
+ ,74.60
+ ,76.45
+ ,73.66
+ ,76.49
+ ,81.58
+ ,82.66
+ ,81.25
+ ,81.69
+ ,81.43
+ ,86.85
+ ,84.88
+ ,81.98
+ ,83.76
+ ,89.19
+ ,87.79
+ ,88.02
+ ,91.51
+ ,91.38
+ ,88.03
+ ,91.54
+ ,89.11
+ ,89.34
+ ,89.03
+ ,85.58
+ ,86.20
+ ,97.88
+ ,04.42
+ ,01.16
+ ,01.07
+ ,05.40
+ ,07.94
+ ,12.79
+ ,09.66
+ ,12.29
+ ,13.93
+ ,97.18
+ ,99.65
+ ,00.10
+ ,00.59
+ ,00.22
+ ,99.29
+ ,93.01
+ ,91.16
+ ,94.94
+ ,96.20
+ ,97.24
+ ,99.87
+ ,95.70
+ ,86.88
+ ,85.38
+ ,82.26
+ ,85.37
+ ,86.45
+ ,87.24
+ ,87.96
+ ,79.85
+ ,79.20
+ ,82.98
+ ,86.80
+ ,87.40
+ ,93.32
+ ,94.26
+ ,98.99)
> 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] 309
> (np <- floor(n / par1))
[1] 25
> 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] 26 26 26 26 26 26 26 26 26 25 25 25
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 61.04 64.26 69.97 66.25 63.28 61.44 64.94 71.47 95.93 90.57 10.14 45.29
 [2,] 64.21 66.35 70.78 63.33 62.03 61.64 65.20 71.98 96.32 90.71 16.69 45.08
 [3,] 63.92 67.43 73.52 60.55 63.43 59.69 65.08 71.09 93.84 88.96 18.52 28.88
 [4,] 68.48 69.45 73.86 62.91 62.41 57.11 64.76 74.04 98.18 91.77 16.32 23.26
 [5,] 67.76 72.98 76.80 59.76 61.05 62.28 68.00 76.70 88.71 95.50 25.96 25.10
 [6,] 65.37 71.58 74.57 58.57 56.31 65.87 69.14 79.10 88.28 98.81 26.29 15.20
 [7,] 61.84 69.97 73.36 59.33 53.30 64.28 70.68 81.62 91.27  1.84 32.19 13.77
 [8,] 59.88 71.93 74.77 60.75 51.99 63.63 72.81 81.66 93.31  5.15 27.35 14.59
 [9,] 62.91 68.53 74.35 59.05 55.42 64.11 73.93 81.22 96.00 10.21 38.54 15.46
[10,] 63.67 71.29 71.14 59.59 59.02 66.46 75.79 83.69 95.98  1.84 34.86  6.23
[11,] 59.96 72.75 72.51 58.97 59.89 61.93 77.02 88.60 97.91  5.62 35.36  1.18
[12,] 62.80 71.64 69.19 59.24 59.39 62.37 75.48 91.86 92.69  6.23 40.21  4.55
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,]  6.89 42.36 45.52 66.31 73.89 76.35 71.19 86.15 78.98 81.25 88.03  7.94
 [2,] 93.88 37.71 46.25 68.44 72.74 77.47 74.13 75.11 78.95 81.69 91.54 12.79
 [3,] 77.70 44.60 52.07 72.04 68.02 76.05 79.81 71.61 80.70 81.43 89.11  9.66
 [4,] 71.85 46.34 52.38 69.55 69.29 75.47 79.66 70.04 75.39 86.85 89.34 12.29
 [5,] 64.15 40.83 52.51 69.16 72.04 69.87 81.50 73.97 73.82 84.88 89.03 13.93
 [6,] 67.81 36.51 52.24 65.63 66.02 73.36 81.24 71.51 75.17 81.98 85.58 97.18
 [7,] 61.04 46.47 50.33 59.89 69.95 78.05 80.68 73.78 74.60 83.76 86.20 99.65
 [8,] 57.04 41.68 51.55 63.56 71.77 79.36 80.00 77.18 76.45 89.19 97.88  0.10
 [9,] 49.93 40.17 53.20 68.05 78.53 82.75 84.87 78.86 73.66 87.79  4.42  0.59
[10,] 54.43 37.51 58.63 69.45 80.50 78.00 84.92 72.14 76.49 88.02  1.16  0.22
[11,] 40.81 40.03 56.34 70.93 77.00 74.54 83.24 76.09 81.58 91.51  1.07 99.29
[12,] 46.28 44.76 61.67 67.51 77.43 72.89 85.12 76.01 82.66 91.38  5.40 93.01
      [,25] [,26]
 [1,] 91.16 87.96
 [2,] 94.94 79.85
 [3,] 96.20 79.20
 [4,] 97.24 82.98
 [5,] 99.87 86.80
 [6,] 95.70 87.40
 [7,] 86.88 93.32
 [8,] 85.38 94.26
 [9,] 82.26 98.99
[10,] 85.37    NA
[11,] 86.45    NA
[12,] 87.24    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]  [,10] [,11]  [,12]
 [1,]  3.17  2.09  0.81 -2.92 -1.25  0.20  0.26  0.51  0.39   0.14  6.55  -0.21
 [2,] -0.29  1.08  2.74 -2.78  1.40 -1.95 -0.12 -0.89 -2.48  -1.75  1.83 -16.20
 [3,]  4.56  2.02  0.34  2.36 -1.02 -2.58 -0.32  2.95  4.34   2.81 -2.20  -5.62
 [4,] -0.72  3.53  2.94 -3.15 -1.36  5.17  3.24  2.66 -9.47   3.73  9.64   1.84
 [5,] -2.39 -1.40 -2.23 -1.19 -4.74  3.59  1.14  2.40 -0.43   3.31  0.33  -9.90
 [6,] -3.53 -1.61 -1.21  0.76 -3.01 -1.59  1.54  2.52  2.99 -96.97  5.90  -1.43
 [7,] -1.96  1.96  1.41  1.42 -1.31 -0.65  2.13  0.04  2.04   3.31 -4.84   0.82
 [8,]  3.03 -3.40 -0.42 -1.70  3.43  0.48  1.12 -0.44  2.69   5.06 11.19   0.87
 [9,]  0.76  2.76 -3.21  0.54  3.60  2.35  1.86  2.47 -0.02  -8.37 -3.68  -9.23
[10,] -3.71  1.46  1.37 -0.62  0.87 -4.53  1.23  4.91  1.93   3.78  0.50  -5.05
[11,]  2.84 -1.11 -3.32  0.27 -0.50  0.44 -1.54  3.26 -5.22   0.61  4.85   3.37
[12,]  1.46 -1.67 -2.94  4.04  2.05  2.57 -4.01  4.07 -2.12   3.91  5.08   2.34
       [,13] [,14] [,15] [,16] [,17] [,18] [,19]  [,20] [,21] [,22]  [,23]
 [1,]  86.99 -4.65  0.73  2.13 -1.15  1.12  2.94 -11.04 -0.03  0.44   3.51
 [2,] -16.18  6.89  5.82  3.60 -4.72 -1.42  5.68  -3.50  1.75 -0.26  -2.43
 [3,]  -5.85  1.74  0.31 -2.49  1.27 -0.58 -0.15  -1.57 -5.31  5.42   0.23
 [4,]  -7.70 -5.51  0.13 -0.39  2.75 -5.60  1.84   3.93 -1.57 -1.97  -0.31
 [5,]   3.66 -4.32 -0.27 -3.53 -6.02  3.49 -0.26  -2.46  1.35 -2.90  -3.45
 [6,]  -6.77  9.96 -1.91 -5.74  3.93  4.69 -0.56   2.27 -0.57  1.78   0.62
 [7,]  -4.00 -4.79  1.22  3.67  1.82  1.31 -0.68   3.40  1.85  5.43  11.68
 [8,]  -7.11 -1.51  1.65  4.49  6.76  3.39  4.87   1.68 -2.79 -1.40 -93.46
 [9,]   4.50 -2.66  5.43  1.40  1.97 -4.75  0.05  -6.72  2.83  0.23  -3.26
[10,] -13.62  2.52 -2.29  1.48 -3.50 -3.46 -1.68   3.95  5.09  3.49  -0.09
[11,]   5.47  4.73  5.33 -3.42  0.43 -1.65  1.88  -0.08  1.08 -0.13   4.33
[12,]  -3.92  0.76  4.64  6.38 -1.08 -1.70  1.03   2.97 -1.41 -3.35   2.54
       [,24] [,25] [,26]
 [1,]   4.85  3.78 -8.11
 [2,]  -3.13  1.26 -0.65
 [3,]   2.63  1.04  3.78
 [4,]   1.64  2.63  3.82
 [5,]  83.25 -4.17  0.60
 [6,]   2.47 -8.82  5.92
 [7,] -99.55 -1.50  0.94
 [8,]   0.49 -3.12  4.73
 [9,]  -0.37  3.11    NA
[10,]  99.07  1.08    NA
[11,]  -6.28  0.79    NA
[12,]  -1.85  0.72    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/wessaorg/rcomp/tmp/1boi71321378043.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/wessaorg/rcomp/tmp/2l0zn1321378043.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/wessaorg/rcomp/tmp/31ute1321378043.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/wessaorg/rcomp/tmp/4akl81321378043.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]
[1,] 42.36 37.71 44.60 46.340 40.830 52.240 32.190 14.59  4.42 34.86 35.36
[2,] 61.04 62.03 60.55 62.410 61.050 65.370 59.330 51.99 49.93 54.43 56.34
[3,] 68.14 71.38 71.35 69.795 70.955 71.545 70.325 71.85 68.29 69.45 72.51
[4,] 81.25 79.85 79.81 79.660 81.500 81.980 81.620 80.00 81.22 78.00 81.58
[5,] 95.93 96.32 96.20 98.180 99.870 98.810 99.650 97.88 98.99 95.98 99.29
     [,12]
[1,] 40.21
[2,] 59.24
[3,] 69.19
[4,] 82.66
[5,] 93.01

$n
 [1] 26 26 26 26 26 26 26 26 26 25 25 25

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 61.87766 65.85823 65.38203 64.44985 64.61829 66.39817 63.41814 63.17072
[2,] 74.40234 76.90177 77.31797 75.14015 77.29171 76.69183 77.23186 80.52928
         [,9]    [,10]    [,11]    [,12]
[1,] 58.59437 62.00188 64.53416 61.78928
[2,] 77.98563 76.89812 80.48584 76.59072

$out
 [1] 10.14  6.89  7.94 16.69 12.79 18.52 28.88  9.66 16.32 23.26 12.29 25.96
[13] 25.10 13.93 26.29 15.20 36.51  1.84 13.77  5.15  0.10  0.59  1.84  6.23
[25]  1.16  0.22  5.62  1.18  1.07  6.23  4.55  5.40

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

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

Warning message:
In bxp(list(stats = c(42.36, 61.04, 68.14, 81.25, 95.93, 37.71,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5fg7i1321378043.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]
[1,] -4.650 -4.72 -5.850 -7.70 -9.90 -6.77 -4.840 -7.110 -9.23 -5.05 -6.28
[2,] -0.210 -2.48 -1.570 -1.57 -3.45 -1.91 -1.310 -1.510 -3.21 -2.29 -1.11
[3,]  0.475 -0.47  0.325  1.74 -0.81  0.03  1.265  0.995  0.54  1.08  0.44
[4,]  2.940  1.75  2.630  3.24  1.35  2.52  2.040  3.430  2.47  2.52  3.26
[5,]  6.550  6.89  5.420  9.64  3.66  5.92  5.430  6.760  5.43  5.09  5.47
     [,12]
[1,] -4.01
[2,] -1.70
[3,]  1.03
[4,]  2.97
[5,]  6.38

$n
 [1] 26 26 26 26 26 26 26 26 25 25 25 25

$conf
         [,1]       [,2]       [,3]      [,4]       [,5]      [,6]      [,7]
[1,] -0.50107 -1.7807226 -0.9764267 0.2495566 -2.2973448 -1.342695 0.2269573
[2,]  1.45107  0.8407226  1.6264267 3.2304434  0.6773448  1.402695 2.3030427
           [,8]     [,9]    [,10]    [,11]    [,12]
[1,] -0.5357257 -1.25488 -0.43996 -0.94092 -0.44572
[2,]  2.5257257  2.33488  2.59996  1.82092  2.50572

$out
 [1]  86.99 -11.04  -8.11 -16.20 -16.18  -9.47  83.25 -96.97   9.96  -8.82
[11]  11.68 -99.55  11.19 -93.46 -13.62  99.07

$group
 [1]  1  1  1  2  2  4  5  6  6  6  7  7  8  8 10 10

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

Warning message:
In bxp(list(stats = c(-4.65, -0.210000000000001, 0.475000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6sgh11321378043.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]  [,10]  [,11]
[1,] 59.88 64.260 69.19 58.570 51.990 59.690 64.760 71.090 88.280  1.840 10.140
[2,] 61.44 67.980 70.96 59.145 55.865 61.540 65.140 73.010 91.980  5.385 17.605
[3,] 63.29 70.630 73.44 59.675 59.640 62.325 69.910 80.160 94.885 49.585 26.820
[4,] 64.79 71.785 74.46 61.830 62.220 64.195 74.705 82.675 96.160 91.240 35.110
[5,] 68.48 72.980 76.80 63.330 63.430 66.460 77.020 91.860 98.180 98.810 40.210
     [,12]  [,13]  [,14] [,15]  [,16]  [,17]  [,18]  [,19]  [,20]  [,21]  [,22]
[1,]  1.18 40.810 36.510 45.52 63.560 66.020 69.870 74.130 70.040 73.660 81.250
[2,] 10.00 48.105 38.870 50.94 65.970 69.620 73.950 79.735 71.875 74.885 81.835
[3,] 15.33 59.040 41.255 52.31 68.245 72.390 76.200 80.960 74.540 76.470 85.865
[4,] 26.99 69.830 44.680 54.77 69.500 77.215 78.025 84.055 76.635 79.840 88.605
[5,] 45.29 93.880 46.470 58.63 72.040 80.500 82.750 85.120 78.860 82.660 91.510
      [,23]  [,24]  [,25] [,26]
[1,]  1.070  0.100 82.260 79.20
[2,]  4.910  4.265 85.915 82.98
[3,] 87.115 12.540 89.200 87.40
[4,] 89.225 95.095 95.950 93.32
[5,] 97.880 99.650 99.870 98.99

$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]  9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 61.76204 68.89451 71.84363 58.45035 56.74144 61.11404 65.54734 75.75173
[2,] 64.81796 72.36549 75.03637 60.89965 62.53856 63.53596 74.27266 84.56827
         [,9]    [,10]    [,11]     [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 92.97847 10.42596 18.83585  7.580747 49.13108 38.60502 50.56311 66.63494
[2,] 96.79153 88.74404 34.80415 23.079253 68.94892 43.90498 54.05689 69.85506
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]     [,23]     [,24]
[1,] 68.92587 74.34137 78.98962 72.36893 74.20999 82.77716  48.65836 -28.88817
[2,] 75.85413 78.05863 82.93038 76.71107 78.73001 88.95284 125.57164  53.96817
        [,25]    [,26]
[1,] 84.62297 81.95427
[2,] 93.77703 92.84573

$out
[1] 66.25 57.11  6.89 61.67 59.89 71.19 86.15

$group
[1]  4  6 13 15 16 19 20

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

Warning message:
In bxp(list(stats = c(59.88, 61.44, 63.29, 64.79, 68.48, 64.26,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7c3um1321378043.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,] 59.05600 68.1400 65.68875
[2,] 63.01198 69.3200 67.75313
[3,] 65.12712 70.6400 70.51875
[4,] 67.19577 71.4625 70.93062
[5,] 69.47769 72.5100 73.61500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 63.21886 69.66279 69.06947
[2,] 67.03537 71.61721 71.96803

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(59.056, 63.0119846153846, 65.1271153846154,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1boi71321378043.ps tmp/1boi71321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/2l0zn1321378043.ps tmp/2l0zn1321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/31ute1321378043.ps tmp/31ute1321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/4akl81321378043.ps tmp/4akl81321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/5fg7i1321378043.ps tmp/5fg7i1321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sgh11321378043.ps tmp/6sgh11321378043.png",intern=TRUE))
character(0)
> try(system("convert tmp/7c3um1321378043.ps tmp/7c3um1321378043.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.578   0.352   2.296