R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(26862,28261,31761,32624,32950,34711,34124,31860,32404,33027,28944,29702,29144,29715,35710,37324,37687,39071,38986,37702,36444,36423,35266,35461,33073,35721,44670,43461,42998,45039,41555,43850,42253,42036,39552,38572,36862,37883,47264,43698,48335,50368,45633,48708,44588,44021,41601,39761,40848,45488,51253,49495,53348,50286,50191,50596,46689,49249,43597,41862)
> par1 = '12'
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 60
> (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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,] 26862 29144 33073 36862 40848   NA
 [2,] 28261 29715 35721 37883 45488   NA
 [3,] 31761 35710 44670 47264 51253   NA
 [4,] 32624 37324 43461 43698 49495   NA
 [5,] 32950 37687 42998 48335 53348   NA
 [6,] 34711 39071 45039 50368 50286   NA
 [7,] 34124 38986 41555 45633 50191   NA
 [8,] 31860 37702 43850 48708 50596   NA
 [9,] 32404 36444 42253 44588 46689   NA
[10,] 33027 36423 42036 44021 49249   NA
[11,] 28944 35266 39552 41601 43597   NA
[12,] 29702 35461 38572 39761 41862   NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,]  1399   571  2648  1021  4640   NA
 [2,]  3500  5995  8949  9381  5765   NA
 [3,]   863  1614 -1209 -3566 -1758   NA
 [4,]   326   363  -463  4637  3853   NA
 [5,]  1761  1384  2041  2033 -3062   NA
 [6,]  -587   -85 -3484 -4735   -95   NA
 [7,] -2264 -1284  2295  3075   405   NA
 [8,]   544 -1258 -1597 -4120 -3907   NA
 [9,]   623   -21  -217  -567  2560   NA
[10,] -4083 -1157 -2484 -2420 -5652   NA
[11,]   758   195  -980 -1840 -1735   NA
[12,]  -558 -2388 -1710  1087    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/www/rcomp/tmp/1xsfl1321391106.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/www/rcomp/tmp/29nw71321391106.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/www/rcomp/tmp/3qoln1321391106.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/www/rcomp/tmp/41pn51321391106.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] [,11] [,12]
[1,] 26862 28261 31761 32624 32950 34711 34124 31860 32404 33027 28944 29702
[2,] 29144 29715 35710 37324 37687 39071 38986 37702 36444 36423 35266 35461
[3,] 33073 35721 44670 43461 42998 45039 41555 43850 42253 42036 39552 38572
[4,] 36862 37883 47264 43698 48335 50286 45633 48708 44588 44021 41601 39761
[5,] 40848 45488 51253 49495 53348 50368 50191 50596 46689 49249 43597 41862

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 27619.48 29949.51 36505.97 38957.15 35474.15 37114.51 36858.25 36073.19
[2,] 38526.52 41492.49 52834.03 47964.85 50521.85 52963.49 46251.75 51626.81
         [,9]    [,10]    [,11]    [,12]
[1,] 36498.47 36667.27 35075.70 35533.63
[2,] 48007.53 47404.73 44028.30 41610.37

$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(26862, 29144, 33073, 36862, 40848, 28261,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5e96h1321391106.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]   [,12]
[1,]  571 3500 -3566 -463 1384 -4735 -2264 -4120 -567 -5652 -1840 -2388.0
[2,] 1021 5765 -1758  326 1384 -3484 -1284 -3907 -217 -4083 -1735 -2049.0
[3,] 1399 5995 -1209  363 1761  -587   405 -1597  -21 -2484  -980 -1134.0
[4,] 2648 8949   863 3853 2033   -95  2295 -1258  623 -2420   195   264.5
[5,] 4640 9381  1614 4637 2041   -85  3075   544  623 -1157   758  1087.0

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

$conf
          [,1]     [,2]      [,3]      [,4]     [,5]      [,6]      [,7]
[1,]  249.3659 3745.194 -3060.992 -2129.169 1302.418 -2981.659 -2123.912
[2,] 2548.6341 8244.806   642.992  2855.169 2219.582  1807.659  2933.912
           [,8]      [,9]     [,10]      [,11]     [,12]
[1,] -3468.7767 -614.5419 -3659.072 -2343.7331 -2961.665
[2,]   274.7767  572.5419 -1308.928   383.7331   693.665

$out
[1] -3062  2560

$group
[1] 5 9

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

Warning message:
In bxp(list(stats = c(571, 1021, 1399, 2648, 4640, 3500, 5765, 5995,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6ht7d1321391106.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]
[1,] 26862.0 35266.0 33073.0 36862.0 40848.0   NA
[2,] 29323.0 35363.5 39062.0 40681.0 44542.5   NA
[3,] 32132.0 36433.5 42144.5 44304.5 49372.0   NA
[4,] 32988.5 37694.5 43655.5 47799.5 50441.0   NA
[5,] 34711.0 39071.0 45039.0 50368.0 53348.0   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 30460.14 35370.32 40049.37 41057.70 46681.65   NA
[2,] 33803.86 37496.68 44239.63 47551.30 52062.35   NA

$out
[1] 29144 29715

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(26862, 29323, 32132, 32988.5, 34711, 35266,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7sgkt1321391106.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,] 33357.8 33073.0 33003.00
[2,] 37431.8 39062.0 38022.25
[3,] 41135.8 42144.5 40513.50
[4,] 42337.4 43655.5 42660.25
[5,] 43895.0 45039.0 44678.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 38898.32 40049.37 38398.08
[2,] 43373.28 44239.63 42628.92

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(33357.8, 37431.8, 41135.8, 42337.4, 43895,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1xsfl1321391106.ps tmp/1xsfl1321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/29nw71321391106.ps tmp/29nw71321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/3qoln1321391106.ps tmp/3qoln1321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/41pn51321391106.ps tmp/41pn51321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/5e96h1321391106.ps tmp/5e96h1321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ht7d1321391106.ps tmp/6ht7d1321391106.png",intern=TRUE))
character(0)
> try(system("convert tmp/7sgkt1321391106.ps tmp/7sgkt1321391106.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.836   0.412   2.276