R version 3.1.0 (2014-04-10) -- "Spring Dance"
Copyright (C) 2014 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 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 216
> (np <- floor(n / par1))
[1] 18
> 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] 18 18 18 18 18 18 18 18 18 18 18 18
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]   [,6]   [,7]   [,8]   [,9]  [,10]  [,11]
 [1,] 71.97 76.23 76.32 79.37 79.18  81.36  96.90  95.17  94.12  88.02  93.03
 [2,] 72.32 76.97 78.92 82.89 80.91  84.81 101.58  95.19  97.86  92.95  95.81
 [3,] 74.07 78.40 80.47 84.43 83.91  90.30 105.46  99.44  98.83  97.05  99.10
 [4,] 77.95 78.60 83.14 85.32 86.30  95.43 110.06 109.18 102.29 101.44 100.12
 [5,] 81.75 80.08 85.42 87.71 89.76  97.59 107.90 109.15 104.49 100.34 100.67
 [6,] 80.81 81.12 81.53 84.68 85.11  97.80 102.87 109.72 102.11  99.98 103.87
 [7,] 74.10 80.31 87.31 80.62 83.81  99.48  96.28 108.41 102.14  94.17 102.39
 [8,] 71.37 84.59 86.01 84.79 85.36  97.52  98.59 102.96 101.28  94.54 107.21
 [9,] 75.21 81.34 85.10 85.49 85.89 104.39 103.22 107.64 101.21  95.12 105.71
[10,] 76.90 80.95 79.91 81.68 82.59  97.74  98.60  97.28  94.20  98.04  99.79
[11,] 74.44 80.48 78.60 77.69 80.87  91.37  91.79  97.25  88.47  93.72  96.12
[12,] 74.76 75.26 78.60 78.31 80.27  92.42  93.83  91.84  88.08  93.83  96.17
       [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18] [,19]
 [1,]  97.23 103.03 104.74 101.70 108.50 106.69 118.40    NA
 [2,]  98.08 104.70 107.50 102.94 113.84 109.87 123.32    NA
 [3,]  99.84 107.45 109.53 104.37 114.59 110.82 123.39    NA
 [4,]  99.72 109.67 109.42 106.93 116.27 114.14 127.04    NA
 [5,]  99.92 110.54 108.60 107.82 113.63 113.31 129.35    NA
 [6,] 102.70 112.05 110.72 110.83 112.29 115.16 127.12    NA
 [7,] 102.06 113.19 105.10 106.86 110.31 111.06 122.10    NA
 [8,] 102.36 114.20 105.19 109.46 108.47 111.13 120.22    NA
 [9,] 102.43 112.56 102.55 108.80 110.67 115.96 121.53    NA
[10,] 100.60 107.36 101.25 108.69 109.10 117.57 119.01    NA
[11,]  98.40 103.93 101.56 107.77 107.02 114.69 114.27    NA
[12,]  98.61 103.83 101.62 108.64 108.12 119.42 114.46    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]   [,8]  [,9] [,10] [,11] [,12]
 [1,]  0.35  0.74  2.60  3.52  1.73  3.45  4.68   0.02  3.74  4.93  2.78  0.85
 [2,]  1.75  1.43  1.55  1.54  3.00  5.49  3.88   4.25  0.97  4.10  3.29  1.76
 [3,]  3.88  0.20  2.67  0.89  2.39  5.13  4.60   9.74  3.46  4.39  1.02 -0.12
 [4,]  3.80  1.48  2.28  2.39  3.46  2.16 -2.16  -0.03  2.20 -1.10  0.55  0.20
 [5,] -0.94  1.04 -3.89 -3.03 -4.65  0.21 -5.03   0.57 -2.38 -0.36  3.20  2.78
 [6,] -6.71 -0.81  5.78 -4.06 -1.30  1.68 -6.59  -1.31  0.03 -5.81 -1.48 -0.64
 [7,] -2.73  4.28 -1.30  4.17  1.55 -1.96  2.31  -5.45 -0.86  0.37  4.82  0.30
 [8,]  3.84 -3.25 -0.91  0.70  0.53  6.87  4.63   4.68 -0.07  0.58 -1.50  0.07
 [9,]  1.69 -0.39 -5.19 -3.81 -3.30 -6.65 -4.62 -10.36 -7.01  2.92 -5.92 -1.83
[10,] -2.46 -0.47 -1.31 -3.99 -1.72 -6.37 -6.81  -0.03 -5.73 -4.32 -3.67 -2.20
[11,]  0.32 -5.22  0.00  0.62 -0.60  1.05  2.04  -5.41 -0.39  0.11  0.05  0.21
[12,]  1.47  1.06  0.77  0.87  1.09  4.48  1.34   2.28 -0.06 -0.80  1.06  4.42
      [,13] [,14] [,15] [,16] [,17] [,18] [,19]
 [1,]  1.67  2.76  1.24  5.34  3.18  4.92    NA
 [2,]  2.75  2.03  1.43  0.75  0.95  0.07    NA
 [3,]  2.22 -0.11  2.56  1.68  3.32  3.65    NA
 [4,]  0.87 -0.82  0.89 -2.64 -0.83  2.31    NA
 [5,]  1.51  2.12  3.01 -1.34  1.85 -2.23    NA
 [6,]  1.14 -5.62 -3.97 -1.98 -4.10 -5.02    NA
 [7,]  1.01  0.09  2.60 -1.84  0.07 -1.88    NA
 [8,] -1.64 -2.64 -0.66  2.20  4.83  1.31    NA
 [9,] -5.20 -1.30 -0.11 -1.57  1.61 -2.52    NA
[10,] -3.43  0.31 -0.92 -2.08 -2.88 -4.74    NA
[11,] -0.10  0.06  0.87  1.10  4.73  0.19    NA
[12,]  0.91  0.08 -0.14 -1.43 -1.02    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/wessaorg/rcomp/tmp/14fvx1400524845.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/2nn0j1400524845.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/3m7841400524845.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/4xybc1400524845.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,]  71.970  72.320  74.07  77.950  80.08  80.810  74.10  71.37  75.210  76.90
[2,]  79.370  82.890  84.43  86.300  89.76  85.110  87.31  86.01  85.890  82.59
[3,]  94.645  96.835  99.27 101.865 102.58 102.785 102.10 101.82 102.885  98.32
[4,] 103.030 104.700 107.45 109.670 109.15 110.830 108.41 108.47 108.800 107.36
[5,] 118.400 123.320 123.39 127.040 129.35 127.120 122.10 120.22 121.530 119.01
      [,11]  [,12]
[1,]  74.44  74.76
[2,]  80.87  80.27
[3,]  94.92  93.83
[4,] 103.93 103.83
[5,] 114.69 119.42

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  85.83379  88.71275  90.69713  93.16179  95.35898  93.20663  94.24216
[2,] 103.45621 104.95725 107.84287 110.56821 109.80102 112.36337 109.95784
          [,8]     [,9]     [,10]     [,11]     [,12]
[1,]  93.45568  94.3531  89.09541  86.33224  85.05603
[2,] 110.18432 111.4169 107.54459 103.50776 102.60397

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(71.97, 79.37, 94.645, 103.03, 118.4, 72.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5pvvc1400524845.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,] 0.02 0.070 -0.120 -2.64 -5.030 -6.71 -5.450 -3.250 -10.36 -6.81 -0.60
[2,] 1.24 1.430  1.020 -0.82 -2.380 -5.02 -1.840 -0.910  -5.20 -4.32 -0.10
[3,] 2.77 1.755  2.615  0.88 -0.075 -1.73  0.195  0.555  -2.91 -2.67  0.15
[4,] 3.74 3.290  3.880  2.28  1.850 -0.64  2.310  3.840  -0.39 -1.31  0.87
[5,] 5.34 5.490  5.130  3.80  3.200  5.78  4.820  6.870   2.92  0.31  2.04
     [,12]
[1,] -1.43
[2,] -0.06
[3,]  0.91
[4,]  1.34
[5,]  2.28

$n
 [1] 18 18 18 18 18 18 18 18 18 18 18 17

$conf
         [,1]     [,2]     [,3]       [,4]      [,5]        [,6]    [,7]
[1,] 1.838976 1.062318 1.549909 -0.2744697 -1.650292 -3.36115392 -1.3505
[2,] 3.701024 2.447682 3.680091  2.0344697  1.500292 -0.09884608  1.7405
          [,8]     [,9]     [,10]      [,11]     [,12]
[1,] -1.213945 -4.70129 -3.790953 -0.2112373 0.3735112
[2,]  2.323945 -1.11871 -1.549047  0.5112373 1.4464888

$out
[1]  9.74 -5.22 -5.41  4.73  4.48  4.42

$group
[1]  3 11 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(0.019999999999996, 1.23999999999999, 2.77,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/631pm1400524845.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 76.32 77.690 79.180  81.360  91.79  91.84  88.080  88.020
[2,] 73.195 77.685 78.76 79.995 80.890  90.835  96.59  96.22  94.160  93.775
[3,] 74.600 80.195 81.00 83.660 83.860  96.475 100.09 101.20 100.020  94.830
[4,] 77.425 81.035 85.26 85.055 85.625  97.770 104.34 108.78 102.125  99.010
[5,] 81.750 84.590 87.31 87.710 89.760 104.390 110.06 109.72 104.490 101.440
       [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18] [,19]
[1,]  93.030  97.230 103.030 101.250 101.700 107.020 106.690 114.270    NA
[2,]  96.145  98.505 104.315 102.085 105.615 108.485 110.940 118.705    NA
[3,]  99.955  99.880 108.560 105.145 107.795 110.490 113.725 121.815    NA
[4,] 103.130 102.210 112.305 109.010 108.745 113.735 115.560 125.215    NA
[5,] 107.210 102.700 114.200 110.720 110.830 116.270 119.420 129.350    NA

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]    [,6]      [,7]     [,8]
[1,] 72.67067 78.66704 78.03531 81.3521 81.70033 93.3119  96.55517  95.4713
[2,] 76.52933 81.72296 83.96469 85.9679 86.01967 99.6381 103.62483 106.9287
          [,9]    [,10]     [,11]     [,12]    [,13]    [,14]    [,15]    [,16]
[1,]  96.38711 92.44228  96.76909  98.19012 104.9157 101.9865 106.3674 108.0954
[2,] 103.65289 97.21772 103.14091 101.56988 112.2043 108.3035 109.2226 112.8846
        [,17]    [,18] [,19]
[1,] 111.6178 118.8457    NA
[2,] 115.8322 124.7843    NA

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

Warning message:
In bxp(list(stats = c(71.37, 73.195, 74.6, 77.425, 81.75, 75.26,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7dfjx1400524845.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,]  92.88667  93.8300  91.28250
[2,]  94.96944  95.8775  93.43438
[3,]  98.36528 100.5450  97.28688
[4,] 100.49556 102.3400  98.80938
[5,] 101.55722 102.8850 100.36500

$n
[1] 12 12 12

$conf
          [,1]      [,2]     [,3]
[1,]  95.84478  97.59741 94.83530
[2,] 100.88577 103.49259 99.73845

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(92.8866666666667, 94.9694444444444, 98.3652777777778,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/14fvx1400524845.ps tmp/14fvx1400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nn0j1400524845.ps tmp/2nn0j1400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/3m7841400524845.ps tmp/3m7841400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/4xybc1400524845.ps tmp/4xybc1400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/5pvvc1400524845.ps tmp/5pvvc1400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/631pm1400524845.ps tmp/631pm1400524845.png",intern=TRUE))
character(0)
> try(system("convert tmp/7dfjx1400524845.ps tmp/7dfjx1400524845.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.282   0.605   3.944