R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
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(71.97
+ ,72.32
+ ,74.07
+ ,77.95
+ ,81.75
+ ,80.81
+ ,74.1
+ ,71.37
+ ,75.21
+ ,76.9
+ ,74.44
+ ,74.76
+ ,76.23
+ ,76.97
+ ,78.4
+ ,78.6
+ ,80.08
+ ,81.12
+ ,80.31
+ ,84.59
+ ,81.34
+ ,80.95
+ ,80.48
+ ,75.26
+ ,76.32
+ ,78.92
+ ,80.47
+ ,83.14
+ ,85.42
+ ,81.53
+ ,87.31
+ ,86.01
+ ,85.1
+ ,79.91
+ ,78.6
+ ,78.6
+ ,79.37
+ ,82.89
+ ,84.43
+ ,85.32
+ ,87.71
+ ,84.68
+ ,80.62
+ ,84.79
+ ,85.49
+ ,81.68
+ ,77.69
+ ,78.31
+ ,79.18
+ ,80.91
+ ,83.91
+ ,86.3
+ ,89.76
+ ,85.11
+ ,83.81
+ ,85.36
+ ,85.89
+ ,82.59
+ ,80.87
+ ,80.27
+ ,81.36
+ ,84.81
+ ,90.3
+ ,95.43
+ ,97.59
+ ,97.8
+ ,99.48
+ ,97.52
+ ,104.39
+ ,97.74
+ ,91.37
+ ,92.42
+ ,96.9
+ ,101.58
+ ,105.46
+ ,110.06
+ ,107.9
+ ,102.87
+ ,96.28
+ ,98.59
+ ,103.22
+ ,98.6
+ ,91.79
+ ,93.83
+ ,95.17
+ ,95.19
+ ,99.44
+ ,109.18
+ ,109.15
+ ,109.72
+ ,108.41
+ ,102.96
+ ,107.64
+ ,97.28
+ ,97.25
+ ,91.84
+ ,94.12
+ ,97.86
+ ,98.83
+ ,102.29
+ ,104.49
+ ,102.11
+ ,102.14
+ ,101.28
+ ,101.21
+ ,94.2
+ ,88.47
+ ,88.08
+ ,88.02
+ ,92.95
+ ,97.05
+ ,101.44
+ ,100.34
+ ,99.98
+ ,94.17
+ ,94.54
+ ,95.12
+ ,98.04
+ ,93.72
+ ,93.83
+ ,93.03
+ ,95.81
+ ,99.1
+ ,100.12
+ ,100.67
+ ,103.87
+ ,102.39
+ ,107.21
+ ,105.71
+ ,99.79
+ ,96.12
+ ,96.17
+ ,97.23
+ ,98.08
+ ,99.84
+ ,99.72
+ ,99.92
+ ,102.7
+ ,102.06
+ ,102.36
+ ,102.43
+ ,100.6
+ ,98.4
+ ,98.61
+ ,103.03
+ ,104.7
+ ,107.45
+ ,109.67
+ ,110.54
+ ,112.05
+ ,113.19
+ ,114.2
+ ,112.56
+ ,107.36
+ ,103.93
+ ,103.83
+ ,104.74
+ ,107.5
+ ,109.53
+ ,109.42
+ ,108.6
+ ,110.72
+ ,105.1
+ ,105.19
+ ,102.55
+ ,101.25
+ ,101.56
+ ,101.62
+ ,101.7
+ ,102.94
+ ,104.37
+ ,106.93
+ ,107.82
+ ,110.83
+ ,106.86
+ ,109.46
+ ,108.8
+ ,108.69
+ ,107.77
+ ,108.64
+ ,108.5
+ ,113.84
+ ,114.59
+ ,116.27
+ ,113.63
+ ,112.29
+ ,110.31
+ ,108.47
+ ,110.67
+ ,109.1
+ ,107.02
+ ,108.12
+ ,106.69
+ ,109.87
+ ,110.82
+ ,114.14
+ ,113.31
+ ,115.16
+ ,111.06
+ ,111.13
+ ,115.96
+ ,117.57
+ ,114.69
+ ,119.42
+ ,118.4
+ ,123.32
+ ,123.39
+ ,127.04
+ ,129.35
+ ,127.12
+ ,122.1
+ ,120.22
+ ,121.53
+ ,119.01
+ ,114.27
+ ,114.46)
> par1 = '18'
> par1 <- '18'
> #'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] 216
> (np <- floor(n / par1))
[1] 12
> 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] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
> arr
       [,1]  [,2]  [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]  [,11]
 [1,] 71.97 80.31 79.37  83.81  96.90 108.41  88.02 102.39 103.03 105.10 108.50
 [2,] 72.32 84.59 82.89  85.36 101.58 102.96  92.95 107.21 104.70 105.19 113.84
 [3,] 74.07 81.34 84.43  85.89 105.46 107.64  97.05 105.71 107.45 102.55 114.59
 [4,] 77.95 80.95 85.32  82.59 110.06  97.28 101.44  99.79 109.67 101.25 116.27
 [5,] 81.75 80.48 87.71  80.87 107.90  97.25 100.34  96.12 110.54 101.56 113.63
 [6,] 80.81 75.26 84.68  80.27 102.87  91.84  99.98  96.17 112.05 101.62 112.29
 [7,] 74.10 76.32 80.62  81.36  96.28  94.12  94.17  97.23 113.19 101.70 110.31
 [8,] 71.37 78.92 84.79  84.81  98.59  97.86  94.54  98.08 114.20 102.94 108.47
 [9,] 75.21 80.47 85.49  90.30 103.22  98.83  95.12  99.84 112.56 104.37 110.67
[10,] 76.90 83.14 81.68  95.43  98.60 102.29  98.04  99.72 107.36 106.93 109.10
[11,] 74.44 85.42 77.69  97.59  91.79 104.49  93.72  99.92 103.93 107.82 107.02
[12,] 74.76 81.53 78.31  97.80  93.83 102.11  93.83 102.70 103.83 110.83 108.12
[13,] 76.23 87.31 79.18  99.48  95.17 102.14  93.03 102.06 104.74 106.86 106.69
[14,] 76.97 86.01 80.91  97.52  95.19 101.28  95.81 102.36 107.50 109.46 109.87
[15,] 78.40 85.10 83.91 104.39  99.44 101.21  99.10 102.43 109.53 108.80 110.82
[16,] 78.60 79.91 86.30  97.74 109.18  94.20 100.12 100.60 109.42 108.69 114.14
[17,] 80.08 78.60 89.76  91.37 109.15  88.47 100.67  98.40 108.60 107.77 113.31
[18,] 81.12 78.60 85.11  92.42 109.72  88.08 103.87  98.61 110.72 108.64 115.16
       [,12] [,13]
 [1,] 111.06    NA
 [2,] 111.13    NA
 [3,] 115.96    NA
 [4,] 117.57    NA
 [5,] 114.69    NA
 [6,] 119.42    NA
 [7,] 118.40    NA
 [8,] 123.32    NA
 [9,] 123.39    NA
[10,] 127.04    NA
[11,] 129.35    NA
[12,] 127.12    NA
[13,] 122.10    NA
[14,] 120.22    NA
[15,] 121.53    NA
[16,] 119.01    NA
[17,] 114.27    NA
[18,] 114.46    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]   [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]  0.35  4.28  3.52  1.55  4.68  -5.45  4.93  4.82  1.67  0.09  5.34  0.07
 [2,]  1.75 -3.25  1.54  0.53  3.88   4.68  4.10 -1.50  2.75 -2.64  0.75  4.83
 [3,]  3.88 -0.39  0.89 -3.30  4.60 -10.36  4.39 -5.92  2.22 -1.30  1.68  1.61
 [4,]  3.80 -0.47  2.39 -1.72 -2.16  -0.03 -1.10 -3.67  0.87  0.31 -2.64 -2.88
 [5,] -0.94 -5.22 -3.03 -0.60 -5.03  -5.41 -0.36  0.05  1.51  0.06 -1.34  4.73
 [6,] -6.71  1.06 -4.06  1.09 -6.59   2.28 -5.81  1.06  1.14  0.08 -1.98 -1.02
 [7,] -2.73  2.60  4.17  3.45  2.31   3.74  0.37  0.85  1.01  1.24 -1.84  4.92
 [8,]  3.84  1.55  0.70  5.49  4.63   0.97  0.58  1.76 -1.64  1.43  2.20  0.07
 [9,]  1.69  2.67 -3.81  5.13 -4.62   3.46  2.92 -0.12 -5.20  2.56 -1.57  3.65
[10,] -2.46  2.28 -3.99  2.16 -6.81   2.20 -4.32  0.20 -3.43  0.89 -2.08  2.31
[11,]  0.32 -3.89  0.62  0.21  2.04  -2.38  0.11  2.78 -0.10  3.01  1.10 -2.23
[12,]  1.47  5.78  0.87  1.68  1.34   0.03 -0.80 -0.64  0.91 -3.97 -1.43 -5.02
[13,]  0.74 -1.30  1.73 -1.96  0.02  -0.86  2.78  0.30  2.76  2.60  3.18 -1.88
[14,]  1.43 -0.91  3.00  6.87  4.25  -0.07  3.29  0.07  2.03 -0.66  0.95  1.31
[15,]  0.20 -5.19  2.39 -6.65  9.74  -7.01  1.02 -1.83 -0.11 -0.11  3.32 -2.52
[16,]  1.48 -1.31  3.46 -6.37 -0.03  -5.73  0.55 -2.20 -0.82 -0.92 -0.83 -4.74
[17,]  1.04  0.00 -4.65  1.05  0.57  -0.39  3.20  0.21  2.12  0.87  1.85  0.19
[18,] -0.81  0.77 -1.30  4.48 -1.31  -0.06 -1.48  4.42 -5.62 -0.14 -4.10    NA
      [,13]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    NA
[13,]    NA
[14,]    NA
[15,]    NA
[16,]    NA
[17,]    NA
[18,]    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/1w3e61393876425.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/283h11393876425.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/3j95f1393876425.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/4rkz51393876425.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]
[1,]  71.970  72.320  74.070  77.950  80.480  75.260  74.100  71.370  75.210
[2,]  82.060  84.975  85.160  83.955  84.730  82.745  80.990  84.800  87.895
[3,]  99.645 102.270 104.005 100.520  98.795  98.075  95.225  97.970  99.335
[4,] 106.755 106.200 107.545 109.865 109.220 107.460 106.005 105.705 107.520
[5,] 111.060 113.840 115.960 117.570 114.690 119.420 118.400 123.320 123.390
       [,10]   [,11]   [,12]   [,13]  [,14]   [,15]  [,16]   [,17]   [,18]
[1,]  76.900  74.440  74.760  76.230  76.97  78.400  78.60  78.600  78.600
[2,]  89.285  88.605  87.680  90.170  90.60  92.100  90.25  89.115  86.595
[3,]  99.160  98.755  99.955 100.770  99.40 101.820 100.36  99.535 101.240
[4,] 107.145 105.755 105.975 105.715 108.48 109.165 109.30 108.875 110.220
[5,] 127.040 129.350 127.120 122.100 120.22 121.530 119.01 114.270 115.160

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  88.38144  92.58914  93.79505  88.70228  87.62495  86.80232  83.81549
[2,] 110.90856 111.95086 114.21495 112.33772 109.96505 109.34768 106.63451
          [,8]      [,9]     [,10]     [,11]     [,12]     [,13]     [,14]
[1,]  88.43509  90.38391  91.01393  90.93277  91.61053  93.67982  91.24481
[2,] 107.50491 108.28609 107.30607 106.57723 108.29947 107.86018 107.55519
         [,15]     [,16]     [,17]     [,18]
[1,]  94.03654  91.67117  90.52233  90.46448
[2,] 109.60346 109.04883 108.54767 112.01552

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

Warning message:
In bxp(list(stats = c(71.97, 82.06, 99.645, 106.755, 111.06, 72.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5qco81393876425.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,] -5.450 -3.250 -5.92 -3.670 -5.410 -6.710 -2.730 -1.64 -5.200 -6.81 -3.890
[2,]  0.220 -0.485 -2.30 -2.400 -4.030 -4.935  0.610  0.64 -2.690 -3.71 -1.165
[3,]  2.595  1.645  1.25 -0.785 -0.770 -0.470  1.775  1.49  2.125 -0.94  0.265
[4,]  4.750  3.990  3.05  0.590  0.055  1.075  3.595  3.02  3.190  2.18  1.570
[5,]  5.340  4.830  4.60  3.800  4.730  2.280  4.920  5.49  5.130  2.31  3.010
      [,12] [,13]  [,14]  [,15]  [,16]  [,17]  [,18]
[1,] -3.970 -1.96 -0.910 -7.010 -6.370 -0.390 -1.480
[2,] -1.115 -1.08  0.000 -3.855 -3.470  0.095 -1.395
[3,]  0.450  0.52  1.370 -0.110 -0.875  0.720 -0.810
[4,]  1.405  2.68  3.145  1.705  0.260  1.450  0.355
[5,]  1.680  3.18  6.870  9.740  3.460  3.200  0.770

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 11

$conf
          [,1]       [,2]      [,3]       [,4]      [,5]      [,6]      [,7]
[1,] 0.5288366 -0.3960775 -1.190171 -2.1487591 -2.633196 -3.211201 0.4135215
[2,] 4.6611634  3.6860775  3.690171  0.5787591  1.093196  2.271201 3.1364785
         [,8]       [,9]     [,10]      [,11]      [,12]     [,13]       [,14]
[1,] 0.404466 -0.5569075 -3.626469 -0.9824519 -0.6993889 -1.194961 -0.06445561
[2,] 2.575534  4.8069075  1.746469  1.5124519  1.5993889  2.234961  2.80445561
         [,15]     [,16]     [,17]       [,18]
[1,] -2.645953 -2.576278 0.1019754 -1.64367887
[2,]  2.425953  0.826278 1.3380246  0.02367887

$out
[1] -10.36   5.78  -5.02  -4.65   4.48   4.42  -5.62  -4.10

$group
[1]  3 12 12 17 18 18 18 18

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18"

Warning message:
In bxp(list(stats = c(-5.45, 0.219999999999999, 2.595, 4.74999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6bl891393876425.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]
[1,] 71.370 75.260 77.69  80.270  91.79  88.08  88.02  96.12 103.03 101.250
[2,] 74.100 78.920 80.62  83.810  96.28  94.20  93.83  98.40 104.74 102.550
[3,] 76.565 80.715 84.17  90.835 100.51 100.02  96.43  99.88 109.01 106.025
[4,] 78.600 84.590 85.32  97.590 107.90 102.29 100.12 102.39 110.72 108.640
[5,] 81.750 87.310 89.76 104.390 110.06 108.41 103.87 107.21 114.20 110.830
       [,11]   [,12] [,13]
[1,] 106.690 111.060    NA
[2,] 108.500 114.690    NA
[3,] 110.745 119.215    NA
[4,] 113.840 123.320    NA
[5,] 116.270 129.350    NA

$n
 [1] 18 18 18 18 18 18 18 18 18 18 18 18  0

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]      [,6]     [,7]      [,8]
[1,] 74.88916 78.60344 82.41968 85.7032  96.1826  97.00721 94.08754  98.39409
[2,] 78.24084 82.82656 85.92032 95.9668 104.8374 103.03279 98.77246 101.36591
        [,9]   [,10]    [,11]    [,12] [,13]
[1,] 106.783 103.757 108.7563 116.0011    NA
[2,] 111.237 108.293 112.7337 122.4289    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(71.37, 74.1, 76.565, 78.6, 81.75, 75.26, 78.92,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7l3pp1393876425.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,]  94.81667  97.970  92.51375
[2,]  97.06000  98.795  95.43125
[3,]  98.10250  99.590  97.59063
[4,]  98.59167 100.770  99.08750
[5,] 100.38833 102.270 102.29125

$n
[1] 18 18 18

$conf
         [,1]      [,2]     [,3]
[1,] 97.53209  98.85449 96.22900
[2,] 98.67291 100.32551 98.95225

$out
[1] 104.005  95.225

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(94.8166666666667, 97.06, 98.1025, 98.5916666666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1w3e61393876425.ps tmp/1w3e61393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/283h11393876425.ps tmp/283h11393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/3j95f1393876425.ps tmp/3j95f1393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/4rkz51393876425.ps tmp/4rkz51393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/5qco81393876425.ps tmp/5qco81393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bl891393876425.ps tmp/6bl891393876425.png",intern=TRUE))
character(0)
> try(system("convert tmp/7l3pp1393876425.ps tmp/7l3pp1393876425.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.122   1.109   6.263