R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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(115.01
+ ,124.56
+ ,128.60
+ ,131.54
+ ,124.93
+ ,126.77
+ ,128.97
+ ,149.91
+ ,149.55
+ ,149.91
+ ,159.47
+ ,162.04
+ ,163.88
+ ,166.82
+ ,173.80
+ ,172.06
+ ,178.21
+ ,174.53
+ ,176.37
+ ,168.65
+ ,166.82
+ ,164.24
+ ,163.51
+ ,164.61
+ ,163.51
+ ,164.61
+ ,167.92
+ ,160.94
+ ,156.53
+ ,159.47
+ ,149.18
+ ,135.58
+ ,124.93
+ ,116.84
+ ,114.27
+ ,114.64
+ ,112.80
+ ,113.17
+ ,111.33
+ ,113.91
+ ,119.42
+ ,120.89
+ ,122.04
+ ,120.15
+ ,118.21
+ ,125.33
+ ,135.37
+ ,146.89
+ ,149.54
+ ,160.84
+ ,169.95
+ ,177.13
+ ,169.35
+ ,159.93
+ ,149.69
+ ,148.67
+ ,136.32
+ ,128.17
+ ,138.74
+ ,140.58
+ ,147.57
+ ,147.83
+ ,155.65
+ ,148.88
+ ,147.90
+ ,141.99
+ ,136.58
+ ,121.82
+ ,127.52
+ ,129.80
+ ,131.29
+ ,135.96
+ ,146.50
+ ,158.65
+ ,153.21
+ ,147.04
+ ,141.04
+ ,140.45
+ ,140.15
+ ,139.30
+ ,137.60
+ ,146.02
+ ,158.79
+ ,167.19
+ ,161.99
+ ,164.62
+ ,156.21
+ ,154.42
+ ,150.39
+ ,148.98
+ ,158.61
+ ,169.98
+ ,190.09
+ ,184.39
+ ,193.67
+ ,203.79
+ ,204.07
+ ,208.92
+ ,206.88
+ ,218.89
+ ,215.52
+ ,251.66
+ ,262.11
+ ,227.27
+ ,202.60
+ ,191.63
+ ,178.71
+ ,178.32
+ ,176.41
+ ,175.70
+ ,175.73
+ ,172.35
+ ,176.61
+ ,183.49
+ ,172.59
+ ,148.39
+ ,138.31
+ ,150.61
+ ,151.74
+ ,151.66
+ ,149.88
+ ,144.62
+ ,137.10
+ ,140.05
+ ,138.92
+ ,130.15
+ ,128.92
+ ,120.64
+ ,118.54
+ ,107.95
+ ,107.93
+ ,126.54
+ ,130.21
+ ,126.21
+ ,125.29
+ ,117.03
+ ,117.34
+ ,113.87
+ ,113.00
+ ,111.41
+ ,103.02
+ ,111.41
+ ,113.19
+ ,108.10
+ ,108.80
+ ,102.16
+ ,105.83
+ ,108.41
+ ,105.70
+ ,105.11
+ ,110.78
+ ,113.51
+ ,108.98
+ ,108.28
+ ,117.49
+ ,128.22
+ ,127.73
+ ,128.01
+ ,132.84
+ ,128.12
+ ,130.28
+ ,129.30
+ ,135.00
+ ,127.23
+ ,123.79
+ ,121.92
+ ,122.03
+ ,123.34
+ ,125.27
+ ,122.53
+ ,125.31
+ ,123.28
+ ,122.56
+ ,123.72
+ ,121.46
+ ,132.03
+ ,149.30
+ ,161.26
+ ,187.84
+ ,190.32
+ ,176.26
+ ,168.98
+ ,149.60
+ ,150.84
+ ,141.81
+ ,138.62
+ ,141.96
+ ,131.35
+ ,131.62
+ ,148.72
+ ,145.62
+ ,147.46
+ ,160.55
+ ,165.57
+ ,166.33
+ ,161.39
+ ,166.28
+ ,166.58
+ ,163.73
+ ,154.74
+ ,150.60
+ ,141.29
+ ,151.03
+ ,150.15
+ ,156.57
+ ,153.87
+ ,153.59
+ ,151.30
+ ,150.99
+ ,140.88
+ ,144.25
+ ,141.93
+ ,143.87
+ ,149.36
+ ,159.71
+ ,167.83
+ ,161.12
+ ,164.44
+ ,167.16
+ ,179.84
+ ,174.44
+ ,180.35
+ ,193.17
+ ,195.16
+ ,202.43
+ ,189.91
+ ,195.98
+ ,212.09
+ ,205.81
+ ,204.31
+ ,196.07
+ ,199.98
+ ,199.10
+ ,198.31
+ ,195.72
+ ,223.04
+ ,238.41
+ ,259.73
+ ,326.54
+ ,335.15
+ ,321.81
+ ,368.62
+ ,369.59
+ ,425.00
+ ,439.72
+ ,362.23
+ ,328.76
+ ,348.55
+ ,328.18
+ ,329.34
+ ,295.55
+ ,237.38
+ ,226.85
+ ,220.14
+ ,239.36
+ ,224.69
+ ,230.98
+ ,233.47
+ ,256.70
+ ,253.41
+ ,224.95
+ ,210.37
+ ,191.09
+ ,198.85
+ ,211.04
+ ,206.25
+ ,201.51
+ ,194.54
+ ,191.07
+ ,192.82
+ ,181.88
+ ,157.67
+ ,195.82
+ ,246.25
+ ,271.69
+ ,270.23
+ ,274.08
+ ,306.53
+ ,326.55
+ ,348.15
+ ,316.75
+ ,336.12
+ ,354.47
+ ,326.43
+ ,303.88
+ ,327.07
+ ,315.92
+ ,289.01
+ ,281.01
+ ,269.03
+ ,274.89
+ ,277.77
+ ,283.88
+ ,266.32
+ ,264.36
+ ,276.19
+ ,345.69
+ ,349.40
+ ,353.42
+ ,358.20
+ ,361.00)
> par1 = '12'
> par1 <- '12'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Mean Plot (v1.0.4) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanplot.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 301
> (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 25 25 25 25 25 25 25 25 25 25 25
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 115.01 163.88 163.51 112.80 149.54 147.57 146.50 161.99 204.07 176.41
 [2,] 124.56 166.82 164.61 113.17 160.84 147.83 158.65 164.62 208.92 175.70
 [3,] 128.60 173.80 167.92 111.33 169.95 155.65 153.21 156.21 206.88 175.73
 [4,] 131.54 172.06 160.94 113.91 177.13 148.88 147.04 154.42 218.89 172.35
 [5,] 124.93 178.21 156.53 119.42 169.35 147.90 141.04 150.39 215.52 176.61
 [6,] 126.77 174.53 159.47 120.89 159.93 141.99 140.45 148.98 251.66 183.49
 [7,] 128.97 176.37 149.18 122.04 149.69 136.58 140.15 158.61 262.11 172.59
 [8,] 149.91 168.65 135.58 120.15 148.67 121.82 139.30 169.98 227.27 148.39
 [9,] 149.55 166.82 124.93 118.21 136.32 127.52 137.60 190.09 202.60 138.31
[10,] 149.91 164.24 116.84 125.33 128.17 129.80 146.02 184.39 191.63 150.61
[11,] 159.47 163.51 114.27 135.37 138.74 131.29 158.79 193.67 178.71 151.74
[12,] 162.04 164.61 114.64 146.89 140.58 135.96 167.19 203.79 178.32 151.66
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 149.88 130.21 108.80 127.73 125.27 176.26 160.55 156.57 161.12 205.81
 [2,] 144.62 126.21 102.16 128.01 122.53 168.98 165.57 153.87 164.44 204.31
 [3,] 137.10 125.29 105.83 132.84 125.31 149.60 166.33 153.59 167.16 196.07
 [4,] 140.05 117.03 108.41 128.12 123.28 150.84 161.39 151.30 179.84 199.98
 [5,] 138.92 117.34 105.70 130.28 122.56 141.81 166.28 150.99 174.44 199.10
 [6,] 130.15 113.87 105.11 129.30 123.72 138.62 166.58 140.88 180.35 198.31
 [7,] 128.92 113.00 110.78 135.00 121.46 141.96 163.73 144.25 193.17 195.72
 [8,] 120.64 111.41 113.51 127.23 132.03 131.35 154.74 141.93 195.16 223.04
 [9,] 118.54 103.02 108.98 123.79 149.30 131.62 150.60 143.87 202.43 238.41
[10,] 107.95 111.41 108.28 121.92 161.26 148.72 141.29 149.36 189.91 259.73
[11,] 107.93 113.19 117.49 122.03 187.84 145.62 151.03 159.71 195.98 326.54
[12,] 126.54 108.10 128.22 123.34 190.32 147.46 150.15 167.83 212.09 335.15
       [,21]  [,22]  [,23]  [,24]  [,25] [,26]
 [1,] 321.81 226.85 211.04 274.08 281.01   361
 [2,] 368.62 220.14 206.25 306.53 269.03    NA
 [3,] 369.59 239.36 201.51 326.55 274.89    NA
 [4,] 425.00 224.69 194.54 348.15 277.77    NA
 [5,] 439.72 230.98 191.07 316.75 283.88    NA
 [6,] 362.23 233.47 192.82 336.12 266.32    NA
 [7,] 328.76 256.70 181.88 354.47 264.36    NA
 [8,] 348.55 253.41 157.67 326.43 276.19    NA
 [9,] 328.18 224.95 195.82 303.88 345.69    NA
[10,] 329.34 210.37 246.25 327.07 349.40    NA
[11,] 295.55 191.09 271.69 315.92 353.42    NA
[12,] 237.38 198.85 270.23 289.01 358.20    NA
> darr
       [,1]  [,2]   [,3]  [,4]   [,5]   [,6]  [,7]  [,8]   [,9]  [,10]  [,11]
 [1,]  9.55  2.94   1.10  0.37  11.30   0.26 12.15  2.63   4.85  -0.71  -5.26
 [2,]  4.04  6.98   3.31 -1.84   9.11   7.82 -5.44 -8.41  -2.04   0.03  -7.52
 [3,]  2.94 -1.74  -6.98  2.58   7.18  -6.77 -6.17 -1.79  12.01  -3.38   2.95
 [4,] -6.61  6.15  -4.41  5.51  -7.78  -0.98 -6.00 -4.03  -3.37   4.26  -1.13
 [5,]  1.84 -3.68   2.94  1.47  -9.42  -5.91 -0.59 -1.41  36.14   6.88  -8.77
 [6,]  2.20  1.84 -10.29  1.15 -10.24  -5.41 -0.30  9.63  10.45 -10.90  -1.23
 [7,] 20.94 -7.72 -13.60 -1.89  -1.02 -14.76 -0.85 11.37 -34.84 -24.20  -8.28
 [8,] -0.36 -1.83 -10.65 -1.94 -12.35   5.70 -1.70 20.11 -24.67 -10.08  -2.10
 [9,]  0.36 -2.58  -8.09  7.12  -8.15   2.28  8.42 -5.70 -10.97  12.30 -10.59
[10,]  9.56 -0.73  -2.57 10.04  10.57   1.49 12.77  9.28 -12.92   1.13  -0.02
[11,]  2.57  1.10   0.37 11.52   1.84   4.67  8.40 10.12  -0.39  -0.08  18.61
[12,]  1.84 -1.10  -1.84  2.65   6.99  10.54 -5.20  0.28  -1.91  -1.78   3.67
      [,12] [,13] [,14]  [,15]  [,16] [,17]  [,18]  [,19]  [,20]  [,21]  [,22]
 [1,] -4.00 -6.64  0.28  -2.74  -7.28  5.02  -2.70   3.32  -1.50  46.81  -6.71
 [2,] -0.92  3.67  4.83   2.78 -19.38  0.76  -0.28   2.72  -8.24   0.97  19.22
 [3,] -8.26  2.58 -4.72  -2.03   1.24 -4.94  -2.29  12.68   3.91  55.41 -14.67
 [4,]  0.31 -2.71  2.16  -0.72  -9.03  4.89  -0.31  -5.40  -0.88  14.72   6.29
 [5,] -3.47 -0.59 -0.98   1.16  -3.19  0.30 -10.11   5.91  -0.79 -77.49   2.49
 [6,] -0.87  5.67  5.70  -2.26   3.34 -2.85   3.37  12.82  -2.59 -33.47  23.23
 [7,] -1.59  2.73 -7.77  10.57 -10.61 -8.99  -2.32   1.99  27.32  19.79  -3.29
 [8,] -8.39 -4.53 -3.44  17.27   0.27 -4.14   1.94   7.27  15.37 -20.37 -28.46
 [9,]  8.39 -0.70 -1.87  11.96  17.10 -9.31   5.49 -12.52  21.32   1.16 -14.58
[10,]  1.78  9.21  0.11  26.58  -3.10  9.74  10.35   6.07  66.81 -33.79 -19.28
[11,] -5.09 10.73  1.31   2.48   1.84 -0.88   8.12  16.11   8.61 -58.17   7.76
[12,]  0.70 -0.49  1.93 -14.06  13.09  6.42  -6.71  -6.28 -13.34 -10.53  12.19
       [,23]  [,24]  [,25] [,26]
 [1,]  -4.79  32.45 -11.98    NA
 [2,]  -4.74  20.02   5.86    NA
 [3,]  -6.97  21.60   2.88    NA
 [4,]  -3.47 -31.40   6.11    NA
 [5,]   1.75  19.37 -17.56    NA
 [6,] -10.94  18.35  -1.96    NA
 [7,] -24.21 -28.04  11.83    NA
 [8,]  38.15 -22.55  69.50    NA
 [9,]  50.43  23.19   3.71    NA
[10,]  25.44 -11.15   4.02    NA
[11,]  -1.46 -26.91   4.78    NA
[12,]   3.85  -8.00   2.80    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/1t21o1355070836.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/2y8pk1355070836.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/3ucxe1355070836.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/4aww01355070836.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]
[1,] 108.800 102.16 105.83 108.41 105.70 105.11 110.78 111.41 103.02 107.95
[2,] 146.500 144.62 137.10 140.05 138.92 130.15 135.00 131.35 127.52 128.17
[3,] 161.555 164.61 166.33 160.94 156.53 159.47 149.69 148.67 149.30 149.91
[4,] 205.810 204.31 196.07 194.54 191.07 192.82 193.17 195.16 202.43 191.63
[5,] 281.010 269.03 274.89 224.69 230.98 266.32 264.36 276.19 303.88 259.73
      [,11]  [,12]
[1,] 107.93 108.10
[2,] 135.37 140.58
[3,] 159.47 164.61
[4,] 193.67 203.79
[5,] 271.69 289.01

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

$conf
        [,1]    [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 143.177 145.748 147.6955 143.7212 140.0506 139.6663 131.3083 128.506
[2,] 179.933 183.472 184.9645 178.1588 173.0094 179.2737 168.0717 168.834
         [,9]    [,10]    [,11]    [,12]
[1,] 125.6284 129.8566 141.0472 144.6356
[2,] 172.9716 169.9634 177.8928 184.5844

$out
 [1] 321.81 361.00 368.62 306.53 369.59 326.55 425.00 348.15 277.77 439.72
[11] 316.75 283.88 362.23 336.12 328.76 354.47 348.55 326.43 328.18 345.69
[21] 329.34 327.07 349.40 326.54 295.55 315.92 353.42 335.15 358.20

$group
 [1]  1  1  2  2  3  3  4  4  4  5  5  5  6  6  7  7  8  8  9  9 10 10 10 11 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(108.8, 146.5, 161.555, 205.81, 281.01, 102.16,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5s4xn1355070836.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,] -11.98 -8.41 -14.67 -9.03 -10.11 -10.94 -28.04 -28.46 -14.58 -12.92 -5.09
[2,]  -4.00 -2.04  -4.94 -4.41  -3.68  -2.85 -10.61 -10.08  -8.09  -0.73 -0.08
[3,]   0.28  0.97  -1.74 -0.98  -0.59  -0.30  -2.32  -1.94   1.16   4.02  2.48
[4,]   4.85  4.83   2.95  4.26   1.84   5.67   2.73   5.70   8.42  10.04  8.40
[5,]  12.15  9.11  12.68 14.72   6.88  18.35  20.94  20.11  23.19  25.44 18.61
      [,12]
[1,] -14.06
[2,]  -5.20
[3,]   0.28
[4,]   3.67
[5,]  13.09

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] -2.5166 -1.20092 -4.23324 -3.71972 -2.33432 -2.99232 -6.53544 -6.92648
[2,]  3.0766  3.14092  0.75324  1.75972  1.15432  2.39232  1.89544  3.04648
         [,9]   [,10]    [,11]    [,12]
[1,] -4.05716 0.61668 -0.19968 -2.52292
[2,]  6.37716 7.42332  5.15968  3.08292

$out
 [1]  46.81  32.45 -19.38  19.22  20.02  55.41  21.60 -31.40  36.14 -77.49
[11]  19.37 -17.56 -33.47  23.23 -34.84  27.32  38.15  69.50  50.43  26.58
[21]  66.81 -33.79 -19.28 -58.17 -26.91

$group
 [1]  1  1  2  2  2  3  3  4  5  5  5  5  6  6  7  7  8  8  9 10 10 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(-11.98, -4.00000000000001, 0.279999999999987,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/656rk1355070836.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,] 115.010 163.510 114.270 111.330 128.170 121.820 137.60 148.980 178.320
[2,] 125.850 164.425 120.885 113.540 139.660 130.545 140.30 155.315 197.115
[3,] 130.255 167.735 152.855 119.785 149.615 139.285 146.26 163.305 207.900
[4,] 149.910 174.165 162.225 123.685 165.095 147.865 155.93 187.240 223.080
[5,] 162.040 178.210 167.920 135.370 177.130 155.650 167.19 203.790 251.660
       [,10]   [,11]   [,12]   [,13]   [,14]  [,15]   [,16]   [,17]   [,18]
[1,] 138.310 107.930 103.020 102.160 121.920 121.46 131.350 141.290 140.880
[2,] 151.135 119.590 111.410 105.765 123.565 122.92 140.215 150.815 144.060
[3,] 172.470 129.535 113.530 108.605 127.870 125.29 146.540 160.970 151.145
[4,] 176.070 139.485 121.315 112.145 129.790 155.28 150.220 165.925 155.220
[5,] 183.490 149.880 130.210 117.490 135.000 190.32 150.840 166.580 167.830
      [,19]   [,20]   [,21]   [,22]   [,23]   [,24]   [,25] [,26]
[1,] 161.12 195.720 295.550 191.090 157.670 274.080 264.360   361
[2,] 170.80 198.705 324.995 215.255 191.945 305.205 271.960   361
[3,] 185.13 205.060 338.945 225.900 198.665 321.590 279.390   361
[4,] 195.57 249.070 369.105 236.415 228.645 331.595 347.545   361
[5,] 212.09 259.730 425.000 256.700 271.690 354.470 358.200   361

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 119.2811 163.2925 133.9995 115.1578 138.0139 131.3852 139.1311 148.7438
[2,] 141.2289 172.1775 171.7105 124.4122 161.2161 147.1848 153.3889 177.8662
         [,9]   [,10]    [,11]    [,12]   [,13]    [,14]    [,15]    [,16]
[1,] 196.0572 161.097 120.4608 109.0123 105.695 125.0307 110.5304 141.9767
[2,] 219.7428 183.843 138.6092 118.0477 111.515 130.7093 140.0496 151.1033
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 154.0782 146.0548 173.8322 182.0882 318.8261 216.2488 181.9259 309.5533
[2,] 167.8618 156.2352 196.4278 228.0318 359.0639 235.5512 215.4041 333.6267
        [,25] [,26]
[1,] 244.9152   361
[2,] 313.8648   361

$out
[1] 146.89 262.11 128.22 176.26 168.98 326.54 335.15 439.72 237.38

$group
[1]  4  9 13 16 16 20 20 21 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" NA  

Warning message:
In bxp(list(stats = c(115.01, 125.85, 130.255, 149.91, 162.04, 163.51,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7o2kl1355070836.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,] 173.7204 148.6700 159.9000
[2,] 177.1292 149.8000 163.6700
[3,] 178.3900 159.4700 164.9850
[4,] 181.3062 163.0825 169.7400
[5,] 184.3420 166.3300 176.0712

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 176.4848 153.4118 162.2164
[2,] 180.2952 165.5282 167.7536

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(173.7204, 177.1292, 178.39, 181.306192307692,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1t21o1355070836.ps tmp/1t21o1355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/2y8pk1355070836.ps tmp/2y8pk1355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ucxe1355070836.ps tmp/3ucxe1355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aww01355070836.ps tmp/4aww01355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/5s4xn1355070836.ps tmp/5s4xn1355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/656rk1355070836.ps tmp/656rk1355070836.png",intern=TRUE))
character(0)
> try(system("convert tmp/7o2kl1355070836.ps tmp/7o2kl1355070836.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.739   0.666   4.412