R version 2.7.2 (2008-08-25)
Copyright (C) 2008 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(110.40,96.40,101.90,106.20,81.00,94.70,101.00,109.40,102.30,90.70,96.20,96.10,106.00,103.10,102.00,104.70,86.00,92.10,106.90,112.60,101.70,92.00,97.40,97.00,105.40,102.70,98.10,104.50,87.40,89.90,109.80,111.70,98.60,96.90,95.10,97.00,112.70,102.90,97.40,111.40,87.40,96.80,114.10,110.30,103.90,101.60,94.60,95.90,104.70,102.80,98.10,113.90,80.90,95.70,113.20,105.90,108.80,102.30,99.00,100.70,115.50)
> par1 = '12'
> x <- x[x>quantile(x,0.05) & x<quantile(x,0.95)] 
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 52
> (np <- floor(n / par1))
[1] 4
> 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 4 4 4 4 4 4 4 4
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] 110.4 103.1  98.1 111.4 108.8
 [2,]  96.4 102.0 104.5  96.8 102.3
 [3,] 101.9 104.7  89.9 110.3  99.0
 [4,] 106.2  92.1 109.8 103.9 100.7
 [5,]  94.7 106.9 111.7 101.6    NA
 [6,] 101.0 112.6  98.6  94.6    NA
 [7,] 109.4 101.7  96.9  95.9    NA
 [8,] 102.3  92.0  95.1 104.7    NA
 [9,]  90.7  97.4  97.0 102.8    NA
[10,]  96.2  97.0 112.7  98.1    NA
[11,]  96.1 105.4 102.9  95.7    NA
[12,] 106.0 102.7  97.4 105.9    NA
> darr
       [,1]  [,2]  [,3]  [,4] [,5]
 [1,] -14.0  -1.1   6.4 -14.6 -6.5
 [2,]   5.5   2.7 -14.6  13.5 -3.3
 [3,]   4.3 -12.6  19.9  -6.4  1.7
 [4,] -11.5  14.8   1.9  -2.3   NA
 [5,]   6.3   5.7 -13.1  -7.0   NA
 [6,]   8.4 -10.9  -1.7   1.3   NA
 [7,]  -7.1  -9.7  -1.8   8.8   NA
 [8,] -11.6   5.4   1.9  -1.9   NA
 [9,]   5.5  -0.4  15.7  -4.7   NA
[10,]  -0.1   8.4  -9.8  -2.4   NA
[11,]   9.9  -2.7  -5.5  10.2   NA
[12,]  -2.9  -4.6  14.0   2.9   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/1p60l1226047046.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/2tgnk1226047046.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/30wok1226047046.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/4yrqi1226047046.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]
[1,]  98.1  96.4  99.0 100.7  94.70  94.6  95.90  92.00  90.70  96.20  95.70
[2,] 103.1  96.8  99.0 100.7  98.15  96.6  96.40  93.55  93.85  96.60  95.90
[3,] 108.8 102.0 101.9 103.9 104.25  99.8  99.30  98.70  97.20  97.55  99.50
[4,] 110.4 102.3 104.7 106.2 109.30 106.8 105.55 103.50 100.10 105.40 104.15
[5,] 111.4 104.5 110.3 109.8 111.70 112.6 109.40 104.70 102.80 112.70 105.40
      [,12]
[1,]  97.40
[2,] 100.05
[3,] 104.30
[4,] 105.95
[5,] 106.00

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

$conf
         [,1]      [,2]     [,3]     [,4]     [,5]    [,6]     [,7]     [,8]
[1,] 103.6418  98.11371  97.8724 100.0137  95.4415  91.742  92.0715  90.8395
[2,] 113.9582 105.88629 105.9276 107.7863 113.0585 107.858 106.5285 106.5605
         [,9]   [,10]    [,11]   [,12]
[1,]  92.2625  90.598  92.9825  99.639
[2,] 102.1375 104.502 106.0175 108.961

$out
[1] 89.9 92.1

$group
[1] 3 4

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

Warning message:
In bxp(list(stats = c(98.1, 103.1, 108.8, 110.4, 111.4, 96.4, 96.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/55dur1226047046.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]
[1,] -14.6 -14.6 -12.6 -11.50 -13.10 -10.90 -9.70 -11.60 -4.70 -9.80 -5.50
[2,] -14.0  -3.3  -6.4  -6.90 -10.05  -6.30 -8.40  -6.75 -2.55 -6.10 -4.10
[3,]  -6.5   2.7   1.7  -0.20  -0.65  -0.20 -4.45   0.00  2.55 -1.25  3.60
[4,]  -1.1   5.5   4.3   8.35   6.00   4.85  3.50   3.65 10.60  4.15 10.05
[5,]   6.4  13.5  19.9  14.80   6.30   8.40  8.80   5.40 15.70  8.40 10.20
             [,12]
[1,] -4.600000e+00
[2,] -3.750000e+00
[3,] -7.105427e-15
[4,]  8.450000e+00
[5,]  1.400000e+01

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

$conf
           [,1]      [,2]      [,3]     [,4]     [,5]    [,6]    [,7]   [,8]
[1,] -15.615108 -3.518058 -5.860593 -12.2475 -13.3295 -9.0085 -13.851 -8.216
[2,]   2.615108  8.918058  9.260593  11.8475  12.0295  8.6085   4.951  8.216
        [,9]   [,10]   [,11]  [,12]
[1,] -7.8385 -9.3475 -7.5785 -9.638
[2,] 12.9385  6.8475 14.7785  9.638

$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(-14.6, -14, -6.5, -1.09999999999999, 6.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6vdkh1226047046.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]
[1,]  90.70  92.00  89.90  94.60  99.00
[2,]  96.15  97.20  96.95  96.35  99.85
[3,] 101.45 102.35  98.35 102.20 101.50
[4,] 106.10 105.05 107.15 105.30 105.55
[5,] 110.40 112.60 112.70 111.40 108.80

$n
[1] 12 12 12 12  4

$conf
          [,1]      [,2]      [,3]      [,4]    [,5]
[1,]  96.91174  98.76956  93.69771  98.11784  96.997
[2,] 105.98826 105.93044 103.00229 106.28216 106.003

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(90.7, 96.15, 101.45, 106.1, 110.4, 92, 97.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7pfdm1226047046.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,]  96.9750  97.200  97.0875
[2,] 100.2125  99.000  99.4125
[3,] 101.0800 100.850 100.4437
[4,] 102.7700 104.075 103.5500
[5,] 106.3600 108.800 106.7500

$n
[1] 12 12 12

$conf
         [,1]      [,2]     [,3]
[1,]  99.9135  98.53526  98.5566
[2,] 102.2465 103.16474 102.3309

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(96.975, 100.2125, 101.08, 102.77, 106.36,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1p60l1226047046.ps tmp/1p60l1226047046.png")
> system("convert tmp/2tgnk1226047046.ps tmp/2tgnk1226047046.png")
> system("convert tmp/30wok1226047046.ps tmp/30wok1226047046.png")
> system("convert tmp/4yrqi1226047046.ps tmp/4yrqi1226047046.png")
> system("convert tmp/55dur1226047046.ps tmp/55dur1226047046.png")
> system("convert tmp/6vdkh1226047046.ps tmp/6vdkh1226047046.png")
> system("convert tmp/7pfdm1226047046.ps tmp/7pfdm1226047046.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.390   0.927   1.718