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(167.72
+ ,192.92
+ ,212.81
+ ,252.28
+ ,286.61
+ ,302.78
+ ,298.51
+ ,291.40
+ ,285.65
+ ,280.06
+ ,268.99
+ ,264.28
+ ,259.43
+ ,256.08
+ ,253.90
+ ,242.28
+ ,237.24
+ ,228.98
+ ,223.90
+ ,209.33
+ ,201.66
+ ,198.37
+ ,198.35
+ ,197.61
+ ,194.02
+ ,191.54
+ ,193.93
+ ,191.13
+ ,187.94
+ ,188.75
+ ,208.90
+ ,233.51
+ ,240.53
+ ,243.93
+ ,246.96
+ ,248.20
+ ,236.15
+ ,242.74
+ ,236.37
+ ,228.26
+ ,231.40
+ ,230.30
+ ,237.94
+ ,235.72
+ ,225.62
+ ,220.83
+ ,208.16
+ ,195.57
+ ,192.71
+ ,178.10
+ ,173.42
+ ,171.84
+ ,170.02
+ ,172.08
+ ,169.63
+ ,169.04
+ ,169.12
+ ,165.73
+ ,160.10
+ ,169.89
+ ,169.25
+ ,181.62
+ ,179.36
+ ,179.72
+ ,181.27
+ ,180.59
+ ,178.59
+ ,181.00
+ ,180.53
+ ,172.27
+ ,173.03
+ ,175.13
+ ,178.08
+ ,180.94
+ ,182.60
+ ,176.81
+ ,176.76
+ ,177.35
+ ,177.41
+ ,180.41
+ ,175.61
+ ,176.09
+ ,178.86
+ ,178.12
+ ,182.85
+ ,189.78
+ ,188.52
+ ,189.55
+ ,191.19
+ ,195.34
+ ,199.33
+ ,197.14
+ ,197.90
+ ,200.22
+ ,199.88
+ ,200.11
+ ,201.55
+ ,205.29
+ ,204.26
+ ,206.65
+ ,208.64
+ ,207.46
+ ,209.78
+ ,206.26
+ ,203.24
+ ,206.21
+ ,201.83
+ ,194.05
+ ,208.76
+ ,283.04
+ ,288.76
+ ,280.55
+ ,275.65
+ ,273.49
+ ,267.79
+ ,255.84
+ ,238.04
+ ,237.52
+ ,240.07
+ ,241.26
+ ,239.75
+ ,233.27
+ ,229.20
+ ,221.25
+ ,212.90
+ ,212.89
+ ,211.20
+ ,202.66
+ ,207.60
+ ,212.51
+ ,198.86
+ ,196.75
+ ,206.73
+ ,205.11
+ ,200.31
+ ,207.19
+ ,206.05
+ ,206.21
+ ,213.89
+ ,211.03
+ ,213.41
+ ,213.54
+ ,212.91
+ ,217.64
+ ,221.47
+ ,233.30
+ ,232.89
+ ,234.03
+ ,240.52
+ ,240.10
+ ,246.43
+ ,242.53
+ ,236.78
+ ,231.15
+ ,217.14
+ ,222.17
+ ,221.92
+ ,228.52
+ ,219.14
+ ,216.34
+ ,223.33
+ ,228.84
+ ,238.03
+ ,236.83
+ ,230.34
+ ,231.66
+ ,224.57
+ ,225.66
+ ,227.34
+ ,231.09
+ ,222.15
+ ,221.49
+ ,218.15
+ ,218.01
+ ,213.02
+ ,206.09
+ ,204.39
+ ,204.02
+ ,201.33
+ ,197.80
+ ,202.10
+ ,207.86
+ ,204.02
+ ,202.78
+ ,208.23
+ ,205.39
+ ,211.89
+ ,211.34
+ ,204.30
+ ,200.36
+ ,190.80
+ ,188.66
+ ,187.96
+ ,198.44
+ ,199.91
+ ,198.80
+ ,202.17
+ ,204.77
+ ,206.77
+ ,195.66
+ ,195.51
+ ,188.26
+ ,190.34
+ ,195.83
+ ,201.30
+ ,194.35
+ ,173.30
+ ,174.57
+ ,177.31
+ ,182.65
+ ,258.23
+ ,352.53
+ ,336.13
+ ,332.72
+ ,336.19
+ ,330.27
+ ,333.84
+ ,318.34
+ ,307.82
+ ,306.20
+ ,301.47
+ ,304.05
+ ,302.81
+ ,308.87
+ ,304.84
+ ,294.99
+ ,289.77
+ ,283.66
+ ,272.47
+ ,271.17
+ ,267.62
+ ,268.26
+ ,259.21
+ ,259.77
+ ,251.27
+ ,236.48
+ ,228.85
+ ,227.42
+ ,221.28
+ ,231.44
+ ,235.55
+ ,218.16
+ ,193.87
+ ,242.02
+ ,239.01
+ ,270.85
+ ,330.09
+ ,340.88
+ ,340.52
+ ,323.75
+ ,319.82
+ ,298.25
+ ,284.87
+ ,272.65
+ ,293.71
+ ,284.63
+ ,288.02
+ ,280.39
+ ,277.90
+ ,272.43
+ ,272.45
+ ,262.39
+ ,254.86
+ ,250.18
+ ,250.30
+ ,252.51
+ ,256.14
+ ,256.62
+ ,261.15
+ ,253.69
+ ,259.56
+ ,262.91
+ ,263.61
+ ,272.54
+ ,249.20
+ ,237.97
+ ,266.52
+ ,285.31
+ ,279.82
+ ,279.09
+ ,274.94
+ ,270.33
+ ,272.15
+ ,275.94
+ ,266.28
+ ,262.94
+ ,254.23
+ ,252.65
+ ,257.10
+ ,262.74
+ ,261.33
+ ,250.43
+ ,243.42
+ ,259.98
+ ,254.47
+ ,248.75
+ ,245.75
+ ,250.61
+ ,244.92
+ ,238.35
+ ,238.52
+ ,240.99
+ ,227.12
+ ,235.75
+ ,245.19
+ ,240.48
+ ,251.29
+ ,248.33
+ ,242.10
+ ,237.20
+ ,234.84
+ ,237.57
+ ,249.45
+ ,243.33
+ ,247.45
+ ,252.35
+ ,243.74
+ ,246.57
+ ,243.89
+ ,244.52
+ ,243.85
+ ,235.66
+ ,230.23
+ ,238.34
+ ,233.66
+ ,229.78
+ ,234.98
+ ,234.61
+ ,225.33
+ ,211.81
+ ,233.04
+ ,232.06
+ ,224.86
+ ,222.61
+ ,208.20
+ ,198.37
+ ,224.90
+ ,200.32
+ ,193.27
+ ,223.82
+ ,235.60
+ ,251.15
+ ,252.76
+ ,267.62
+ ,268.88
+ ,259.11
+ ,239.95
+ ,234.08
+ ,218.47
+ ,222.08
+ ,240.83
+ ,238.76
+ ,247.93
+ ,235.82
+ ,234.18
+ ,232.25
+ ,233.66
+ ,231.58
+ ,221.40
+ ,209.12)
> 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] 360
> (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] 30 30 30 30 30 30 30 30 30 30 30 30
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 167.72 259.43 194.02 236.15 192.71 169.25 178.08 182.85 201.55 208.76
 [2,] 192.92 256.08 191.54 242.74 178.10 181.62 180.94 189.78 205.29 283.04
 [3,] 212.81 253.90 193.93 236.37 173.42 179.36 182.60 188.52 204.26 288.76
 [4,] 252.28 242.28 191.13 228.26 171.84 179.72 176.81 189.55 206.65 280.55
 [5,] 286.61 237.24 187.94 231.40 170.02 181.27 176.76 191.19 208.64 275.65
 [6,] 302.78 228.98 188.75 230.30 172.08 180.59 177.35 195.34 207.46 273.49
 [7,] 298.51 223.90 208.90 237.94 169.63 178.59 177.41 199.33 209.78 267.79
 [8,] 291.40 209.33 233.51 235.72 169.04 181.00 180.41 197.14 206.26 255.84
 [9,] 285.65 201.66 240.53 225.62 169.12 180.53 175.61 197.90 203.24 238.04
[10,] 280.06 198.37 243.93 220.83 165.73 172.27 176.09 200.22 206.21 237.52
[11,] 268.99 198.35 246.96 208.16 160.10 173.03 178.86 199.88 201.83 240.07
[12,] 264.28 197.61 248.20 195.57 169.89 175.13 178.12 200.11 194.05 241.26
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 239.75 206.73 221.47 221.92 227.34 202.10 187.96 201.30 333.84 272.47
 [2,] 233.27 205.11 233.30 228.52 231.09 207.86 198.44 194.35 318.34 271.17
 [3,] 229.20 200.31 232.89 219.14 222.15 204.02 199.91 173.30 307.82 267.62
 [4,] 221.25 207.19 234.03 216.34 221.49 202.78 198.80 174.57 306.20 268.26
 [5,] 212.90 206.05 240.52 223.33 218.15 208.23 202.17 177.31 301.47 259.21
 [6,] 212.89 206.21 240.10 228.84 218.01 205.39 204.77 182.65 304.05 259.77
 [7,] 211.20 213.89 246.43 238.03 213.02 211.89 206.77 258.23 302.81 251.27
 [8,] 202.66 211.03 242.53 236.83 206.09 211.34 195.66 352.53 308.87 236.48
 [9,] 207.60 213.41 236.78 230.34 204.39 204.30 195.51 336.13 304.84 228.85
[10,] 212.51 213.54 231.15 231.66 204.02 200.36 188.26 332.72 294.99 227.42
[11,] 198.86 212.91 217.14 224.57 201.33 190.80 190.34 336.19 289.77 221.28
[12,] 196.75 217.64 222.17 225.66 197.80 188.66 195.83 330.27 283.66 231.44
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30]
 [1,] 235.55 284.87 250.30 266.52 257.10 238.52 249.45 233.66 224.90 218.47
 [2,] 218.16 272.65 252.51 285.31 262.74 240.99 243.33 229.78 200.32 222.08
 [3,] 193.87 293.71 256.14 279.82 261.33 227.12 247.45 234.98 193.27 240.83
 [4,] 242.02 284.63 256.62 279.09 250.43 235.75 252.35 234.61 223.82 238.76
 [5,] 239.01 288.02 261.15 274.94 243.42 245.19 243.74 225.33 235.60 247.93
 [6,] 270.85 280.39 253.69 270.33 259.98 240.48 246.57 211.81 251.15 235.82
 [7,] 330.09 277.90 259.56 272.15 254.47 251.29 243.89 233.04 252.76 234.18
 [8,] 340.88 272.43 262.91 275.94 248.75 248.33 244.52 232.06 267.62 232.25
 [9,] 340.52 272.45 263.61 266.28 245.75 242.10 243.85 224.86 268.88 233.66
[10,] 323.75 262.39 272.54 262.94 250.61 237.20 235.66 222.61 259.11 231.58
[11,] 319.82 254.86 249.20 254.23 244.92 234.84 230.23 208.20 239.95 221.40
[12,] 298.25 250.18 237.97 252.65 238.35 237.57 238.34 198.37 234.08 209.12
      [,31]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]  [,6]  [,7]  [,8]  [,9]  [,10]  [,11]
 [1,]  25.20  -3.35  -2.48   6.59 -14.61 12.37  2.86  6.93  3.74  74.28  -6.48
 [2,]  19.89  -2.18   2.39  -6.37  -4.68 -2.26  1.66 -1.26 -1.03   5.72  -4.07
 [3,]  39.47 -11.62  -2.80  -8.11  -1.58  0.36 -5.79  1.03  2.39  -8.21  -7.95
 [4,]  34.33  -5.04  -3.19   3.14  -1.82  1.55 -0.05  1.64  1.99  -4.90  -8.35
 [5,]  16.17  -8.26   0.81  -1.10   2.06 -0.68  0.59  4.15 -1.18  -2.16  -0.01
 [6,]  -4.27  -5.08  20.15   7.64  -2.45 -2.00  0.06  3.99  2.32  -5.70  -1.69
 [7,]  -7.11 -14.57  24.61  -2.22  -0.59  2.41  3.00 -2.19 -3.52 -11.95  -8.54
 [8,]  -5.75  -7.67   7.02 -10.10   0.08 -0.47 -4.80  0.76 -3.02 -17.80   4.94
 [9,]  -5.59  -3.29   3.40  -4.79  -3.39 -8.26  0.48  2.32  2.97  -0.52   4.91
[10,] -11.07  -0.02   3.03 -12.67  -5.63  0.76  2.77 -0.34 -4.38   2.55 -13.65
[11,]  -4.71  -0.74   1.24 -12.59   9.79  2.10 -0.74  0.23 -7.78   1.19  -2.11
[12,]  -4.85  -3.59 -12.05  -2.86  -0.64  2.95  4.73  1.44 14.71  -1.51   9.98
      [,12]  [,13] [,14] [,15] [,16]  [,17]  [,18]  [,19]  [,20]  [,21]  [,22]
 [1,] -1.62  11.83  6.60  3.75  5.76  10.48  -6.95 -15.50  -1.30 -17.39 -12.22
 [2,] -4.80  -0.41 -9.38 -8.94 -3.84   1.47 -21.05 -10.52  -3.55 -24.29  21.06
 [3,]  6.88   1.14 -2.80 -0.66 -1.24  -1.11   1.27  -1.62   0.64  48.15  -9.08
 [4,] -1.14   6.49  6.99 -3.34  5.45   3.37   2.74  -4.73  -9.05  -3.01   3.39
 [5,]  0.16  -0.42  5.51 -0.14 -2.84   2.60   5.34   2.58   0.56  31.84  -7.63
 [6,]  7.68   6.33  9.19 -4.99  6.50   2.00  75.58  -1.24  -8.50  59.24  -2.49
 [7,] -2.86  -3.90 -1.20 -6.93 -0.55 -11.11  94.30   6.06 -14.79  10.79  -5.47
 [8,]  2.38  -5.75 -6.49 -1.70 -7.04  -0.15 -16.40  -4.03  -7.63  -0.36   0.02
 [9,]  0.13  -5.63  1.32 -0.37 -3.94  -7.25  -3.41  -9.85  -1.43 -16.77 -10.06
[10,] -0.63 -14.01 -7.09 -2.69 -9.56   2.08   3.47  -5.22  -6.14  -3.93  -7.53
[11,]  4.73   5.03  1.09 -3.53 -2.14   5.49  -5.92  -6.11  10.16 -21.57  -4.68
[12,]  3.83  -0.25  1.68  4.30 -0.70   5.47   3.57 -11.19   4.11 -13.38   0.12
       [,23] [,24]  [,25]  [,26] [,27]  [,28]  [,29]  [,30] [,31]
 [1,]   2.21 18.79   5.64   2.47 -6.12  -3.88 -24.58   3.61    NA
 [2,]   3.63 -5.49  -1.41 -13.87  4.12   5.20  -7.05  18.75    NA
 [3,]   0.48 -0.73 -10.90   8.63  4.90  -0.37  30.55  -2.07    NA
 [4,]   4.53 -4.15  -7.01   9.44 -8.61  -9.28  11.78   9.17    NA
 [5,]  -7.46 -4.61  16.56  -4.71  2.83 -13.52  15.55 -12.11    NA
 [6,]   5.87  1.82  -5.51  10.81 -2.68  21.23   1.61  -1.64    NA
 [7,]   3.35  3.79  -5.72  -2.96  0.63  -0.98  14.86  -1.93    NA
 [8,]   0.70 -9.66  -3.00  -6.23 -0.67  -7.20   1.26   1.41    NA
 [9,]   8.93 -3.34   4.86  -4.90 -8.19  -2.25  -9.77  -2.08    NA
[10,] -23.34 -8.71  -5.69  -2.36 -5.43 -14.41 -19.16 -10.18    NA
[11,] -11.23 -1.58  -6.57   2.73  8.11  -9.83  -5.87 -12.28    NA
[12,]  28.55  4.45   0.17  11.88 -4.68  26.53 -15.61     NA    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/fisher/rcomp/tmp/1956c1352751729.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/fisher/rcomp/tmp/2ckfg1352751729.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/fisher/rcomp/tmp/3hgtj1352751729.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/fisher/rcomp/tmp/4fcje1352751729.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,] 167.72 178.10 173.300 171.840 170.02 172.08 169.630 169.04 169.12 165.730
[2,] 201.30 198.44 193.930 202.780 206.05 205.39 211.200 206.26 204.30 204.020
[3,] 223.41 229.15 224.635 231.145 233.50 229.64 237.985 236.10 232.00 231.365
[4,] 249.45 252.51 253.900 252.280 247.93 259.77 258.230 262.91 263.61 259.110
[5,] 284.87 318.34 307.820 306.200 301.47 304.05 302.810 340.88 340.52 332.720
      [,11]   [,12]
[1,] 160.10 169.890
[2,] 199.88 196.750
[3,] 221.34 223.915
[4,] 246.96 241.260
[5,] 289.77 298.250

$n
 [1] 30 30 30 30 30 30 30 30 30 30 30 30

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 209.5203 213.5526 207.3356 216.8659 221.419 213.9532 224.4184 219.7583
[2,] 237.2997 244.7474 241.9344 245.4241 245.581 245.3268 251.5516 252.4417
        [,9]    [,10]   [,11]    [,12]
[1,] 214.891 215.4733 207.759 211.0753
[2,] 249.109 247.2567 234.921 236.7547

$out
[1] 333.84 330.09 352.53 336.19 319.82 330.27

$group
[1]  1  7  8 11 11 12

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

> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/5w21c1352751729.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]
[1,] -24.580 -13.87 -8.210 -9.28 -8.260 -8.500 -14.79 -16.40 -10.060 -23.34
[2,]  -6.120  -6.37 -2.800 -4.73 -2.840 -2.490  -5.72  -7.04  -5.630 -10.18
[3,]   2.665  -2.22 -0.695  0.75  0.075  1.715  -2.06  -3.01  -3.315  -5.53
[4,]   6.600   2.39  1.270  4.53  2.830  7.640   3.00   0.08   0.480  -0.34
[5,]  25.200   5.72  6.880 11.78  5.510 21.230  14.86   7.02   8.930   3.47
       [,11]  [,12]
[1,] -12.590 -13.38
[2,]  -6.110  -2.86
[3,]  -1.845   1.44
[4,]   2.100   4.45
[5,]  10.160  14.71

$n
 [1] 30 30 30 30 30 30 30 30 30 30 30 29

$conf
          [,1]       [,2]       [,3]      [,4]      [,5]      [,6]       [,7]
[1,] -1.004303 -4.7469728 -1.8690616 -1.921206 -1.560609 -1.207173 -4.5754341
[2,]  6.334303  0.3069728  0.4790616  3.421206  1.710609  4.637173  0.4554341
           [,8]      [,9]     [,10]      [,11]     [,12]
[1,] -5.0638866 -5.077535 -8.368517 -4.2133158 -0.704744
[2,] -0.9561134 -1.552465 -2.691483  0.5233158  3.584744

$out
 [1]  74.28  19.89 -21.05 -24.29  21.06  18.75  39.47 -11.62  48.15  -9.08
[11] -10.90   8.63  30.55  34.33  16.17  31.84  16.56 -13.52  15.55 -12.11
[21]  75.58  59.24  24.61  94.30 -17.80 -16.77 -21.57  28.55  26.53 -15.61

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

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

> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/6tmac1352751729.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,] 167.720 197.610 187.940 208.160 165.730 169.250 175.610 182.850 201.550
[2,] 232.545 200.015 191.335 223.225 169.080 174.080 176.785 189.665 202.535
[3,] 274.525 226.440 201.460 230.850 169.955 179.540 177.745 196.240 205.750
[4,] 289.005 248.090 242.230 236.260 172.750 180.795 179.635 199.605 207.055
[5,] 302.780 259.430 248.200 242.740 178.100 181.620 182.600 200.220 209.780
       [,10]   [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18]
[1,] 208.760 196.750 200.310 217.140 216.340 197.800 200.360 187.960 173.300
[2,] 239.055 205.130 206.130 226.660 222.625 204.205 201.230 192.925 179.980
[3,] 261.815 212.700 209.110 233.665 227.090 215.515 204.160 197.135 229.765
[4,] 278.100 225.225 213.475 240.310 231.000 221.820 208.045 201.040 334.425
[5,] 288.760 239.750 217.640 246.430 238.030 231.090 211.890 206.770 352.530
       [,19]   [,20]  [,21]   [,22]   [,23]   [,24]   [,25]   [,26]   [,27]
[1,] 283.660 221.280 193.87 250.180 237.970 252.650 238.350 227.120 235.660
[2,] 298.230 230.145 237.28 267.410 251.405 264.610 245.335 236.475 240.835
[3,] 304.445 255.240 284.55 275.275 256.380 271.240 250.520 239.500 243.870
[4,] 308.345 267.940 326.92 284.750 262.030 277.515 258.540 243.645 247.010
[5,] 318.340 272.470 340.88 293.710 272.540 285.310 262.740 251.290 252.350
       [,28]   [,29]   [,30] [,31]
[1,] 198.370 193.270 209.120    NA
[2,] 217.210 224.360 221.740    NA
[3,] 227.555 237.775 232.955    NA
[4,] 233.350 255.935 237.290    NA
[5,] 234.980 268.880 247.930    NA

$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  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 248.7732 204.5127 178.2464 224.9046 168.2811 176.4772 176.4451 191.7063
[2,] 300.2768 248.3673 224.6736 236.7954 171.6289 182.6028 179.0449 200.7737
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 203.6884 244.0063 203.5345 205.7599 227.4391 223.2701 207.4807 201.0516
[2,] 207.8116 279.6237 221.8655 212.4601 239.8909 230.9099 223.5493 207.2684
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 193.4337 159.3216 299.8315 238.0014 243.6646 267.3661 251.5339 265.3539
[2,] 200.8363 300.2084 309.0585 272.4786 325.4354 283.1839 261.2261 277.1261
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30] [,31]
[1,] 244.4971 236.2297 241.0535 220.1934 223.3734 225.8625    NA
[2,] 256.5429 242.7703 246.6865 234.9166 252.1766 240.0475    NA

$out
[1] 195.57 192.71 160.10 194.05 190.80 188.66 333.84 230.23

$group
[1]  4  5  5  9 16 16 19 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(167.72, 232.545, 274.525, 289.005, 302.78,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/7lism1352751729.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,] 224.9660 221.3400 218.7488
[2,] 226.9313 224.2750 224.0250
[3,] 229.4742 230.3925 227.3512
[4,] 234.6377 232.7500 231.8281
[5,] 239.6453 237.9850 234.3313

$n
[1] 12 12 12

$conf
         [,1]    [,2]     [,3]
[1,] 225.9593 226.527 223.7922
[2,] 232.9891 234.258 230.9103

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(224.966, 226.931333333333, 229.474166666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1956c1352751729.ps tmp/1956c1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ckfg1352751729.ps tmp/2ckfg1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/3hgtj1352751729.ps tmp/3hgtj1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/4fcje1352751729.ps tmp/4fcje1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/5w21c1352751729.ps tmp/5w21c1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/6tmac1352751729.ps tmp/6tmac1352751729.png",intern=TRUE))
character(0)
> try(system("convert tmp/7lism1352751729.ps tmp/7lism1352751729.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.580   0.781   4.361