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(103.48,103.93,103.89,104.40,104.79,104.77,105.13,105.26,104.96,104.75,105.01,105.15,105.20,105.77,105.78,106.26,106.13,106.12,106.57,106.44,106.54,107.10,108.10,108.40,108.84,109.62,110.42,110.67,111.66,112.28,112.87,112.18,112.36,112.16,111.49,111.25,111.36,111.74,111.10,111.33,111.25,111.04,110.97,111.31,111.02,111.07,111.36,111.54,112.05,112.52,112.94,113.33,113.78,113.77,113.82,113.89,114.25,114.41,114.55,115.00,115.66,116.33,116.91,117.20,117.59,117.95,118.09,117.99,118.31,118.49,118.96,119.01,119.88,120.59,120.85,120.93,120.89,120.61,120.83,121.36,121.57,121.79)
> par1 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 82
> (np <- floor(n / par1))
[1] 6
> 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] 7 7 7 7 7 7 7 7 7 7 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]
 [1,] 103.48 105.20 108.84 111.36 112.05 115.66 119.88
 [2,] 103.93 105.77 109.62 111.74 112.52 116.33 120.59
 [3,] 103.89 105.78 110.42 111.10 112.94 116.91 120.85
 [4,] 104.40 106.26 110.67 111.33 113.33 117.20 120.93
 [5,] 104.79 106.13 111.66 111.25 113.78 117.59 120.89
 [6,] 104.77 106.12 112.28 111.04 113.77 117.95 120.61
 [7,] 105.13 106.57 112.87 110.97 113.82 118.09 120.83
 [8,] 105.26 106.44 112.18 111.31 113.89 117.99 121.36
 [9,] 104.96 106.54 112.36 111.02 114.25 118.31 121.57
[10,] 104.75 107.10 112.16 111.07 114.41 118.49 121.79
[11,] 105.01 108.10 111.49 111.36 114.55 118.96     NA
[12,] 105.15 108.40 111.25 111.54 115.00 119.01     NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,]  0.45  0.57  0.78  0.38  0.47  0.67  0.71
 [2,] -0.04  0.01  0.80 -0.64  0.42  0.58  0.26
 [3,]  0.51  0.48  0.25  0.23  0.39  0.29  0.08
 [4,]  0.39 -0.13  0.99 -0.08  0.45  0.39 -0.04
 [5,] -0.02 -0.01  0.62 -0.21 -0.01  0.36 -0.28
 [6,]  0.36  0.45  0.59 -0.07  0.05  0.14  0.22
 [7,]  0.13 -0.13 -0.69  0.34  0.07 -0.10  0.53
 [8,] -0.30  0.10  0.18 -0.29  0.36  0.32  0.21
 [9,] -0.21  0.56 -0.20  0.05  0.16  0.18  0.22
[10,]  0.26  1.00 -0.67  0.29  0.14  0.47    NA
[11,]  0.14  0.30 -0.24  0.18  0.45  0.05    NA
[12,]  0.05  0.44  0.11  0.51  0.66  0.87    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/1t1v21353454059.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/29dx51353454059.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/3sxd71353454059.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/4mpnb1353454059.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,] 103.480 103.930 103.890 104.400 104.790 104.77 105.130 105.260 104.96
[2,] 107.020 107.695 108.100 108.465 108.690 108.58 108.770 108.875 108.78
[3,] 111.360 111.740 111.100 111.330 111.660 112.28 112.870 112.180 112.36
[4,] 113.855 114.425 114.925 115.265 115.685 115.86 115.955 115.940 116.28
[5,] 119.880 120.590 120.850 120.930 120.890 120.61 120.830 121.360 121.57
       [,10]   [,11]   [,12]
[1,] 104.750 105.010 105.150
[2,] 109.085 108.100 108.400
[3,] 112.160 111.425 111.395
[4,] 116.450 114.550 115.000
[5,] 121.790 118.960 119.010

$n
 [1] 7 7 7 7 7 7 7 7 7 7 6 6

$conf
         [,1]    [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 107.2782 107.721 107.0242 107.2691 107.4827 107.9325 108.5792 107.9609
[2,] 115.4418 115.759 115.1758 115.3909 115.8373 116.6275 117.1608 116.3991
         [,9]    [,10]    [,11]    [,12]
[1,] 107.8811 107.7617 107.2645 107.1378
[2,] 116.8389 116.5583 115.5855 115.6522

$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(103.48, 107.02, 111.36, 113.855, 119.88, 103.93,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/5gh9z1353454059.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.38 -0.640 0.080 -0.13 -0.280 -0.070 -0.130 -0.300 -0.210 0.140 -0.24
[2,] 0.46 -0.015 0.240 -0.06 -0.115  0.095 -0.115 -0.095 -0.075 0.140  0.05
[3,] 0.57  0.260 0.290  0.39 -0.010  0.220  0.070  0.180  0.160 0.275  0.16
[4,] 0.69  0.500 0.435  0.42  0.175  0.405  0.235  0.265  0.200 0.470  0.30
[5,] 0.78  0.800 0.510  0.99  0.360  0.590  0.530  0.360  0.560 0.470  0.45
     [,12]
[1,] 0.050
[2,] 0.110
[3,] 0.475
[4,] 0.660
[5,] 0.870

$n
 [1] 7 7 7 7 7 7 7 7 7 6 6 6

$conf
          [,1]        [,2]      [,3]      [,4]       [,5]     [,6]       [,7]
[1,] 0.4326477 -0.04754969 0.1735491 0.1033517 -0.1831833 0.034873 -0.1390144
[2,] 0.7073523  0.56754969 0.4064509 0.6766483  0.1631833 0.405127  0.2790144
            [,8]         [,9]      [,10]        [,11]     [,12]
[1,] -0.03498619 -0.004225564 0.06213934 -0.001258075 0.1202322
[2,]  0.39498619  0.324225564 0.48786066  0.321258075 0.8297678

$out
[1]  0.62 -0.69  1.00 -0.67

$group
[1]  5  7 10 10

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

Warning message:
In bxp(list(stats = c(0.379999999999995, 0.460000000000001, 0.569999999999993,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/6lrtp1353454059.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]
[1,] 103.480 105.200 108.840 110.970 112.050 115.660 119.88
[2,] 104.165 105.950 110.545 111.055 113.135 117.055 120.61
[3,] 104.780 106.350 111.575 111.280 113.800 117.970 120.87
[4,] 105.070 106.835 112.230 111.360 114.330 118.400 121.36
[5,] 105.260 108.100 112.870 111.740 115.000 119.010 121.79

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]
[1,] 104.3672 105.9463 110.8065 111.1409 113.255 117.3565 120.4953
[2,] 105.1928 106.7537 112.3435 111.4191 114.345 118.5835 121.2447

$out
[1] 108.4

$group
[1] 2

$names
[1] "1" "2" "3" "4" "5" "6" NA 

Warning message:
In bxp(list(stats = c(103.48, 104.165, 104.78, 105.07, 105.26, 105.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/fisher/rcomp/tmp/7odqz1353454059.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,] 110.9243 111.1000 110.4375
[2,] 111.6385 111.3775 111.4312
[3,] 112.1579 111.7000 112.0263
[4,] 112.6221 112.2300 112.3850
[5,] 112.8243 112.8700 112.7675

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 111.7092 111.3112 111.5912
[2,] 112.6065 112.0888 112.4613

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(110.924285714286, 111.638452380952, 112.157857142857,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1t1v21353454059.ps tmp/1t1v21353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/29dx51353454059.ps tmp/29dx51353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/3sxd71353454059.ps tmp/3sxd71353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/4mpnb1353454059.ps tmp/4mpnb1353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/5gh9z1353454059.ps tmp/5gh9z1353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/6lrtp1353454059.ps tmp/6lrtp1353454059.png",intern=TRUE))
character(0)
> try(system("convert tmp/7odqz1353454059.ps tmp/7odqz1353454059.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.451   0.811   4.257