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(-9,-13,-18,-11,-9,-10,-13,-11,-5,-15,-6,-6,-3,-1,-3,-4,-6,0,-4,-2,-2,-6,-7,-6,-6,-3,-2,-5,-11,-11,-11,-10,-14,-8,-9,-5,-1,-2,-5,-4,-6,-2,-2,-2,-2,2,1,-8,-1,1,-1,2,2,1,-1,-2,-2,-1,-8,-4,-6,-3,-3,-7,-9,-11,-13,-11,-9,-17,-22,-25,-20,-24,-24,-22,-19,-18,-17,-11,-11)
> 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] 81
> (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 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   -9   -3   -6   -1   -1   -6  -20
 [2,]  -13   -1   -3   -2    1   -3  -24
 [3,]  -18   -3   -2   -5   -1   -3  -24
 [4,]  -11   -4   -5   -4    2   -7  -22
 [5,]   -9   -6  -11   -6    2   -9  -19
 [6,]  -10    0  -11   -2    1  -11  -18
 [7,]  -13   -4  -11   -2   -1  -13  -17
 [8,]  -11   -2  -10   -2   -2  -11  -11
 [9,]   -5   -2  -14   -2   -2   -9  -11
[10,]  -15   -6   -8    2   -1  -17   NA
[11,]   -6   -7   -9    1   -8  -22   NA
[12,]   -6   -6   -5   -8   -4  -25   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   -4    2    3   -1    2    3   -4
 [2,]   -5   -2    1   -3   -2    0    0
 [3,]    7   -1   -3    1    3   -4    2
 [4,]    2   -2   -6   -2    0   -2    3
 [5,]   -1    6    0    4   -1   -2    1
 [6,]   -3   -4    0    0   -2   -2    1
 [7,]    2    2    1    0   -1    2    6
 [8,]    6    0   -4    0    0    2    0
 [9,]  -10   -4    6    4    1   -8   NA
[10,]    9   -1   -1   -1   -7   -5   NA
[11,]    0    1    4   -9    4   -3   NA
[12,]    3    0    4    7   -2    5   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/18cz71259748988.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/2duom1259748988.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/3cr6v1259748988.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/4p6e41259748988.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,] -9.0 -13.0 -24.0  -11  -11  -18  -17  -11  -14   -17  -9.0    -8
[2,] -7.5  -8.0 -11.5   -9  -10  -11  -13  -11  -10   -15  -9.0    -8
[3,] -6.0  -3.0  -3.0   -5   -9  -10  -11  -10   -5    -7  -7.5    -6
[4,] -2.0  -1.5  -2.5   -4   -6   -1   -3   -2   -2    -1  -6.0    -5
[5,] -1.0   1.0  -1.0    2   -6    1   -1   -2   -2     2  -6.0    -4

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

$conf
          [,1]       [,2]      [,3]      [,4]       [,5]       [,6]       [,7]
[1,] -9.284511 -6.8816951 -8.374655 -7.985919 -11.388735 -15.971839 -16.971839
[2,] -2.715489  0.8816951  2.374655 -2.014081  -6.611265  -4.028161  -5.028161
           [,8]       [,9]      [,10]     [,11]     [,12]
[1,] -15.374655 -9.7774709 -16.030452 -9.435097 -7.935097
[2,]  -4.625345 -0.2225291   2.030452 -5.564903 -4.064903

$out
[1] -20 -24 -22   2 -19   1 -22 -25

$group
[1]  1  2  4  5  5 11 11 12

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

Warning message:
In bxp(list(stats = c(-9, -7.5, -6, -2, -1, -13, -8, -3, -1.5, 1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5v82z1259748988.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,] -4.0 -5.0 -4.0   -6 -2.0 -4.0 -1.0    0 -10.0    -7  -9.0  -2.0
[2,] -2.5 -2.5 -2.0   -2 -1.0 -2.5  0.5    0  -8.0    -5  -3.0   0.0
[3,]  2.0 -2.0  1.0   -2  0.0 -2.0  2.0    0  -1.5    -1   0.5   3.5
[4,]  2.5  0.0  2.5    1  2.5  0.0  2.0    1   4.0    -1   4.0   5.0
[5,]  3.0  1.0  7.0    3  6.0  1.0  2.0    2   6.0    -1   4.0   7.0

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

$conf
           [,1]       [,2]      [,3]       [,4]      [,5]       [,6]     [,7]
[1,] -0.9859193 -3.4929597 -1.687327 -3.7915516 -2.090144 -3.4929597 1.104224
[2,]  4.9859193 -0.5070403  3.687327 -0.2084484  2.090144 -0.5070403 2.895776
           [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -0.5971839 -9.240388 -3.580129 -4.015226 0.2748385
[2,]  0.5971839  6.240388  1.580129  5.015226 6.7251615

$out
[1]  6  6 -4  9

$group
[1]  7  8  8 10

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

Warning message:
In bxp(list(stats = c(-4, -2.5, 2, 2.5, 3, -5, -2.5, -2, 0, 1, -4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/694xs1259748988.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]
[1,] -18.0 -7.0 -14.0 -8.0   -4 -25.0  -24
[2,] -13.0 -6.0 -11.0 -4.5   -2 -15.0  -22
[3,] -10.5 -3.5  -8.5 -2.0   -1 -10.0  -19
[4,]  -7.5 -2.0  -5.0 -1.5    1  -6.5  -17
[5,]  -5.0  0.0  -2.0  2.0    2  -3.0  -11

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

$conf
           [,1]      [,2]      [,3]       [,4]       [,5]       [,6]      [,7]
[1,] -13.008587 -5.324427 -11.23664 -3.3683201 -2.3683201 -13.876907 -21.63333
[2,]  -7.991413 -1.675573  -5.76336 -0.6316799  0.3683201  -6.123093 -16.36667

$out
[1] -8

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(-18, -13, -10.5, -7.5, -5, -7, -6, -3.5, -2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7o7wv1259748988.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,] -9.000000 -11.0 -8.000
[2,] -8.392857  -9.5 -7.625
[3,] -7.392857  -6.5 -6.500
[4,] -6.785714  -5.0 -6.000
[5,] -6.428571  -3.0 -4.750

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] -8.125886 -8.55248 -7.241173
[2,] -6.659828 -4.44752 -5.758827

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(-9, -8.39285714285714, -7.39285714285714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/18cz71259748988.ps tmp/18cz71259748988.png")
> system("convert tmp/2duom1259748988.ps tmp/2duom1259748988.png")
> system("convert tmp/3cr6v1259748988.ps tmp/3cr6v1259748988.png")
> system("convert tmp/4p6e41259748988.ps tmp/4p6e41259748988.png")
> system("convert tmp/5v82z1259748988.ps tmp/5v82z1259748988.png")
> system("convert tmp/694xs1259748988.ps tmp/694xs1259748988.png")
> system("convert tmp/7o7wv1259748988.ps tmp/7o7wv1259748988.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.194   0.940   1.553