R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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(10407.00,10463.00,10556.00,10646.00,10702.00,11353.00,11346.00,11451.00,11964.00,12574.00,13031.00,13812.00,14544.00,14931.00,14886.00,16005.00,17064.00,15168.00,16050.00,15839.00,15137.00,14954.00,15648.00,15305.00,15579.00,16348.00,15928.00,16171.00,15937.00,15713.00,15594.00,15683.00,16438.00,17032.00,17696.00,17745.00,19394.00,20148.00,20108.00,18584.00,18441.00,18391.00,19178.00,18079.00,18483.00,19644.00,19195.00,19650.00,20830.00,23595.00,22937.00,21814.00,21928.00,21777.00,21383.00,21467.00,22052.00,22680.00,24320.00,24977.00,25204.00)
> par1 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 61
> (np <- floor(n / par1))
[1] 5
> 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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]
 [1,] 10407 14544 15579 19394 20830 25204
 [2,] 10463 14931 16348 20148 23595    NA
 [3,] 10556 14886 15928 20108 22937    NA
 [4,] 10646 16005 16171 18584 21814    NA
 [5,] 10702 17064 15937 18441 21928    NA
 [6,] 11353 15168 15713 18391 21777    NA
 [7,] 11346 16050 15594 19178 21383    NA
 [8,] 11451 15839 15683 18079 21467    NA
 [9,] 11964 15137 16438 18483 22052    NA
[10,] 12574 14954 17032 19644 22680    NA
[11,] 13031 15648 17696 19195 24320    NA
[12,] 13812 15305 17745 19650 24977    NA
> darr
      [,1]  [,2] [,3]  [,4]  [,5] [,6]
 [1,]   56   387  769   754  2765   NA
 [2,]   93   -45 -420   -40  -658   NA
 [3,]   90  1119  243 -1524 -1123   NA
 [4,]   56  1059 -234  -143   114   NA
 [5,]  651 -1896 -224   -50  -151   NA
 [6,]   -7   882 -119   787  -394   NA
 [7,]  105  -211   89 -1099    84   NA
 [8,]  513  -702  755   404   585   NA
 [9,]  610  -183  594  1161   628   NA
[10,]  457   694  664  -449  1640   NA
[11,]  781  -343   49   455   657   NA
[12,]  732   274 1649  1180   227   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/www/html/rcomp/tmp/16pk31289902713.ps",horizontal=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/www/html/rcomp/tmp/26pk31289902713.ps",horizontal=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/www/html/rcomp/tmp/36pk31289902713.ps",horizontal=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/www/html/rcomp/tmp/4yh1o1289902713.ps",horizontal=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] [,11] [,12]
[1,] 10407.0 10463 10556 16005 15937 11353 11346 15683 11964 12574 13031 13812
[2,] 14544.0 14931 14886 16005 15937 15168 15594 15683 15137 14954 15648 15305
[3,] 17486.5 16348 15928 16171 17064 15713 16050 15839 16438 17032 17696 17745
[4,] 20830.0 20148 20108 18584 18441 18391 19178 18079 18483 19644 19195 19650
[5,] 25204.0 23595 22937 21814 21928 21777 21383 21467 22052 22680 24320 24977

$n
 [1] 6 5 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 13431.83 12661.68 12238.15 14348.69 15294.68 13435.64 13517.55 14145.99
[2,] 21541.17 20034.32 19617.85 17993.31 18833.32 17990.36 18582.45 17532.01
         [,9]    [,10]   [,11]    [,12]
[1,] 14073.72 13718.06 15189.7 14674.83
[2,] 18802.28 20345.94 20202.3 20815.17

$out
[1] 10646 10702 11451

$group
[1] 4 5 8

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

Warning message:
In bxp(list(stats = c(10407, 14544, 17486.5, 20830, 25204, 10463,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5yh1o1289902713.ps",horizontal=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] [,12]
[1,]   56 -658 -1524 -234 -224 -394 -211  404  594   457  -343   227
[2,]  387 -420 -1123 -143 -224 -119 -211  404  594   457    49   274
[3,]  754  -45    90   56 -151   -7   84  513  610   664   455   732
[4,]  769  -40   243  114  -50  787   89  585  628   694   657  1180
[5,]  769   93  1119  114  -50  882  105  755  628   694   781  1649

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
          [,1]      [,2]      [,3]      [,4]       [,5]      [,6]      [,7]
[1,]  484.0798 -313.5070 -875.2122 -125.5956 -273.94796 -647.1773 -127.9792
[2,] 1023.9202  223.5070 1055.2122  237.5956  -28.05204  633.1773  295.9792
         [,8]     [,9]    [,10]     [,11]      [,12]
[1,] 385.1059 585.9757 496.5364  25.38873   91.82268
[2,] 640.8941 634.0243 831.4636 884.61127 1372.17732

$out
 [1]  2765  1059   651 -1896 -1099  -702  -183  1161  -449  1640

$group
 [1]  1  4  5  5  7  8  9  9 10 10

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

Warning message:
In bxp(list(stats = c(56, 387, 754, 769, 769, -658, -420, -45, -40,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/698191289902713.ps",horizontal=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]
[1,] 10407.0 14544.0 15579 18079.0 20830 25204
[2,] 10601.0 14942.5 15698 18462.0 21622 25204
[3,] 11349.5 15236.5 16054 19186.5 21990 25204
[4,] 12269.0 15922.0 16735 19647.0 23266 25204
[5,] 13812.0 17064.0 17745 20148.0 24977 25204

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]  [,6]
[1,] 10588.71 14789.74 15581.02 18646.01 21240.16 25204
[2,] 12110.29 15683.26 16526.98 19726.99 22739.84 25204

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(10407, 10601, 11349.5, 12269, 13812, 14544,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/798191289902713.ps",horizontal=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,] 16480.40 15713.00 16779.50
[2,] 16677.10 15989.00 17035.00
[3,] 16848.90 16393.00 17342.50
[4,] 17518.23 17275.25 17487.25
[5,] 18297.80 17745.00 17636.88

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 16465.25 15806.33 17136.23
[2,] 17232.55 16979.67 17548.77

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(16480.4, 16677.1, 16848.9, 17518.2333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/16pk31289902713.ps tmp/16pk31289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/26pk31289902713.ps tmp/26pk31289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/36pk31289902713.ps tmp/36pk31289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/4yh1o1289902713.ps tmp/4yh1o1289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/5yh1o1289902713.ps tmp/5yh1o1289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/698191289902713.ps tmp/698191289902713.png",intern=TRUE))
character(0)
> try(system("convert tmp/798191289902713.ps tmp/798191289902713.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.206   0.938   1.562