R version 2.10.1 (2009-12-14)
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(132.8,132.5,131.4,131.4,130.7,131.5,131.2,130.1,130.5,129,128.2,128.4,127.3,127.7,127,123.9,125.4,124.6,124.5,124.8,124.1,124.2,122.8,122.3,121.1,121.7,122.2,122.2,122.7,121.7,121,119.8,120.2,116.6,116,118,117.1,116.2,113.3,114.3,113.6,113,112.9,112.7,112.5,113,111.9,110.9,109.8,108.3,109.2,109.2,108.7,109.8,110.8,110,109.6,109.5,110.8,111.6,113.1,114.3,114.1,113.8,112.6,112.7,111.5,110.7,110.4,109.7,110,111.3,109,108.2,107.2,108.7,110.3,110.3,109.5,109.5,109.4,109.6,111.3,110,109.5,110.693,109.195,108.095,108.199,106.87,105.278,108.711,111.192,109.641,109.42,109.935,111.126,110.733,110.34,111.766,111.294,111.54,112.008,111.007,114.963,112.045,110.703,108.894,107.51,111.35,112.964,115.203,115.182,115.191,112.346,110.774,113.07,111.138,109.092,107.971,107.051)
> par1 = '12'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 121
> (np <- floor(n / par1))
[1] 10
> 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] 11 10 10 10 10 10 10 10 10 10 10 10
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]    [,8]    [,9]   [,10]   [,11]
 [1,] 132.8 127.3 121.1 117.1 109.8 113.1 109.0 109.500 111.126 107.510 107.051
 [2,] 132.5 127.7 121.7 116.2 108.3 114.3 108.2 110.693 110.733 111.350      NA
 [3,] 131.4 127.0 122.2 113.3 109.2 114.1 107.2 109.195 110.340 112.964      NA
 [4,] 131.4 123.9 122.2 114.3 109.2 113.8 108.7 108.095 111.766 115.203      NA
 [5,] 130.7 125.4 122.7 113.6 108.7 112.6 110.3 108.199 111.294 115.182      NA
 [6,] 131.5 124.6 121.7 113.0 109.8 112.7 110.3 106.870 111.540 115.191      NA
 [7,] 131.2 124.5 121.0 112.9 110.8 111.5 109.5 105.278 112.008 112.346      NA
 [8,] 130.1 124.8 119.8 112.7 110.0 110.7 109.5 108.711 111.007 110.774      NA
 [9,] 130.5 124.1 120.2 112.5 109.6 110.4 109.4 111.192 114.963 113.070      NA
[10,] 129.0 124.2 116.6 113.0 109.5 109.7 109.6 109.641 112.045 111.138      NA
[11,] 128.2 122.8 116.0 111.9 110.8 110.0 111.3 109.420 110.703 109.092      NA
[12,] 128.4 122.3 118.0 110.9 111.6 111.3 110.0 109.935 108.894 107.971      NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]   [,8]   [,9]  [,10] [,11]
 [1,] -0.3  0.4  0.6 -0.9 -1.5  1.2 -0.8  1.193 -0.393  3.840    NA
 [2,] -1.1 -0.7  0.5 -2.9  0.9 -0.2 -1.0 -1.498 -0.393  1.614    NA
 [3,]  0.0 -3.1  0.0  1.0  0.0 -0.3  1.5 -1.100  1.426  2.239    NA
 [4,] -0.7  1.5  0.5 -0.7 -0.5 -1.2  1.6  0.104 -0.472 -0.021    NA
 [5,]  0.8 -0.8 -1.0 -0.6  1.1  0.1  0.0 -1.329  0.246  0.009    NA
 [6,] -0.3 -0.1 -0.7 -0.1  1.0 -1.2 -0.8 -1.592  0.468 -2.845    NA
 [7,] -1.1  0.3 -1.2 -0.2 -0.8 -0.8  0.0  3.433 -1.001 -1.572    NA
 [8,]  0.4 -0.7  0.4 -0.2 -0.4 -0.3 -0.1  2.481  3.956  2.296    NA
 [9,] -1.5  0.1 -3.6  0.5 -0.1 -0.7  0.2 -1.551 -2.918 -1.932    NA
[10,] -0.8 -1.4 -0.6 -1.1  1.3  0.3  1.7 -0.221 -1.342 -2.046    NA
[11,]  0.2 -0.5  2.0 -1.0  0.8  1.3 -1.3  0.515 -1.809 -1.121    NA
[12,] -1.1 -1.2 -0.9 -1.1  1.5 -2.3 -0.5  1.191 -1.384 -0.920    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/1yhfs1272317599.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/rcomp/tmp/29qwv1272317599.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/rcomp/tmp/39qwv1272317599.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/rcomp/tmp/49qwv1272317599.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]
[1,] 107.051 108.200 107.200 108.095 108.199 106.87 105.278 108.7110 109.400
[2,] 109.250 110.693 109.200 109.200 110.300 110.30 110.800 110.0000 110.400
[3,] 111.126 112.825 113.132 114.050 113.100 112.85 112.177 110.8905 112.785
[4,] 119.100 121.700 122.200 122.200 122.700 121.70 121.000 119.8000 120.200
[5,] 132.800 132.500 131.400 131.400 130.700 131.50 131.200 130.1000 130.500
        [,10]   [,11]   [,12]
[1,] 109.5000 109.092 107.971
[2,] 109.6410 110.000 109.935
[3,] 111.5915 111.050 111.100
[4,] 116.6000 116.000 118.000
[5,] 124.2000 122.800 128.400

$n
 [1] 11 10 10 10 10 10 10 10 10 10 10 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 106.4336 107.3255 106.6367 107.5547 106.9045 107.1541 107.0807 105.9940
[2,] 115.8184 118.3245 119.6273 120.5453 119.2955 118.5459 117.2733 115.7870
         [,9]    [,10]    [,11]    [,12]
[1,] 107.8885 108.1145 108.0522 107.0704
[2,] 117.6815 115.0685 114.0478 115.1296

$out
[1] 129.0 128.2

$group
[1] 10 11

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

Warning message:
In bxp(list(stats = c(107.051, 109.25, 111.126, 119.1, 132.8, 108.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5cagt1272317600.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]
[1,] -1.500 -2.9000 -1.100 -1.2000 -1.3290 -2.845 -1.572 -0.700 -3.600 -2.046
[2,] -0.800 -1.1000 -0.300 -0.7000 -0.8000 -1.200 -1.100 -0.300 -1.932 -1.342
[3,]  0.050 -0.5465  0.000 -0.2465  0.0045 -0.500 -0.800  0.150 -1.100 -0.700
[4,]  1.193  0.5000  1.426  0.5000  0.2460 -0.100  0.000  2.296  0.100  0.300
[5,]  3.840  1.6140  2.239  1.6000  1.1000  1.000  0.300  3.956  0.500  1.700
      [,11]  [,12]
[1,] -1.809 -1.384
[2,] -1.121 -1.200
[3,] -0.150 -1.010
[4,]  0.800 -0.500
[5,]  2.000 -0.500

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

$conf
           [,1]       [,2]       [,3]       [,4]       [,5]        [,6]
[1,] -0.9457823 -1.3459238 -0.8623784 -0.8460678 -0.5181233 -1.04960386
[2,]  1.0457823  0.2529238  0.8623784  0.3530678  0.5271233  0.04960386
           [,7]      [,8]        [,9]      [,10]      [,11]      [,12]
[1,] -1.3496039 -1.147065 -2.11526822 -1.5204087 -1.1098082 -1.3597479
[2,] -0.2503961  1.447065 -0.08473178  0.1204087  0.8098082 -0.6602521

$out
[1] -3.100  3.433  1.500 -2.300  1.191

$group
[1]  3  7 12 12 12

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

Warning message:
In bxp(list(stats = c(-1.5, -0.799999999999997, 0.0499999999999972,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6cagt1272317600.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]   [,7]     [,8]    [,9]   [,10]
[1,] 128.20 122.30 116.00 110.90 108.3 109.70 107.20 106.8700 110.340 107.510
[2,] 129.55 124.00 118.90 112.60 109.2 110.55 108.85 108.1470 110.718 109.933
[3,] 130.95 124.55 121.05 113.00 109.7 112.05 109.50 109.3075 111.210 111.848
[4,] 131.45 126.20 121.95 113.95 110.4 113.45 110.15 109.7880 111.887 114.126
[5,] 132.80 127.70 122.70 114.30 111.6 114.30 111.30 111.1920 112.045 115.203
       [,11]
[1,] 107.051
[2,] 107.051
[3,] 107.051
[4,] 107.051
[5,] 107.051

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 130.0834 123.5466 119.6589 112.3843 109.1527 110.7273 108.9071 108.5590
[2,] 131.8166 125.5534 122.4411 113.6157 110.2473 113.3727 110.0929 110.0560
         [,9]    [,10]   [,11]
[1,] 110.6768 109.9355 107.051
[2,] 111.7432 113.7605 107.051

$out
[1] 117.100 116.200 105.278 114.963 108.894

$group
[1] 4 4 8 9 9

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

Warning message:
In bxp(list(stats = c(128.2, 129.55, 130.95, 131.45, 132.8, 122.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7cagt1272317600.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,] 113.9300 110.8905 112.5754
[2,] 114.6258 111.1130 113.6378
[3,] 115.3478 112.4810 114.7872
[4,] 115.7883 112.9750 115.2438
[5,] 116.1676 114.0500 115.6845

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 114.8176 111.6317 114.0547
[2,] 115.8781 113.3303 115.5197

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(113.93, 114.6258, 115.34785, 115.78825, 116.1676,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1yhfs1272317599.ps tmp/1yhfs1272317599.png",intern=TRUE))
character(0)
> try(system("convert tmp/29qwv1272317599.ps tmp/29qwv1272317599.png",intern=TRUE))
character(0)
> try(system("convert tmp/39qwv1272317599.ps tmp/39qwv1272317599.png",intern=TRUE))
character(0)
> try(system("convert tmp/49qwv1272317599.ps tmp/49qwv1272317599.png",intern=TRUE))
character(0)
> try(system("convert tmp/5cagt1272317600.ps tmp/5cagt1272317600.png",intern=TRUE))
character(0)
> try(system("convert tmp/6cagt1272317600.ps tmp/6cagt1272317600.png",intern=TRUE))
character(0)
> try(system("convert tmp/7cagt1272317600.ps tmp/7cagt1272317600.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.620   1.400   2.745