R version 2.8.0 (2008-10-20)
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(129.5,103.7,114.2,123.5,111.3,126.5,119.6,116.5,116.7,119.4,124,130.6,120.1,113.2,111.1,126,115.8,111,128.7,112.6,114.7,118.5,124.8,128.6,127,111.8,100.6,122.9,117.8,108.1,129.6,111.4,110,115.2,118.8,116.2,126.3,106.7,96.5,119.1,109.6,110.3,118.8,104.5,107.7,127.7,118.5,120.1,127.4,107.8,106.5,124.6,101.9,106.5,119.4,103.3,99.6,120.9,111.7,123.9)
> 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] 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,] 129.5 120.1 127.0 126.3 127.4   NA
 [2,] 103.7 113.2 111.8 106.7 107.8   NA
 [3,] 114.2 111.1 100.6  96.5 106.5   NA
 [4,] 123.5 126.0 122.9 119.1 124.6   NA
 [5,] 111.3 115.8 117.8 109.6 101.9   NA
 [6,] 126.5 111.0 108.1 110.3 106.5   NA
 [7,] 119.6 128.7 129.6 118.8 119.4   NA
 [8,] 116.5 112.6 111.4 104.5 103.3   NA
 [9,] 116.7 114.7 110.0 107.7  99.6   NA
[10,] 119.4 118.5 115.2 127.7 120.9   NA
[11,] 124.0 124.8 118.8 118.5 111.7   NA
[12,] 130.6 128.6 116.2 120.1 123.9   NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,] -25.8  -6.9 -15.2 -19.6 -19.6   NA
 [2,]  10.5  -2.1 -11.2 -10.2  -1.3   NA
 [3,]   9.3  14.9  22.3  22.6  18.1   NA
 [4,] -12.2 -10.2  -5.1  -9.5 -22.7   NA
 [5,]  15.2  -4.8  -9.7   0.7   4.6   NA
 [6,]  -6.9  17.7  21.5   8.5  12.9   NA
 [7,]  -3.1 -16.1 -18.2 -14.3 -16.1   NA
 [8,]   0.2   2.1  -1.4   3.2  -3.7   NA
 [9,]   2.7   3.8   5.2  20.0  21.3   NA
[10,]   4.6   6.3   3.6  -9.2  -9.2   NA
[11,]   6.6   3.8  -2.6   1.6  12.2   NA
[12,] -10.5  -1.6  10.1   7.3    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/html/rcomp/tmp/1nwzc1229416536.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/2wd9h1229416536.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/3r0be1229416536.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/40l0w1229416536.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,] 126.3 103.7  96.5 122.9 101.9 106.5 118.8 103.3  99.6 115.2 111.7 116.2
[2,] 126.3 106.7 100.6 122.9 109.6 108.1 119.4 104.5 107.7 118.5 118.5 120.1
[3,] 127.0 107.8 106.5 123.5 111.3 110.3 119.6 111.4 110.0 119.4 118.8 123.9
[4,] 127.4 111.8 111.1 124.6 115.8 111.0 128.7 112.6 114.7 120.9 124.0 128.6
[5,] 127.4 113.2 114.2 126.0 117.8 111.0 129.6 116.5 116.7 120.9 124.8 130.6

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

$conf
         [,1]     [,2]      [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 126.2227 104.1964  99.08073 122.2988 106.9191 108.2509 113.0286 105.6766
[2,] 127.7773 111.4036 113.91927 124.7012 115.6809 112.3491 126.1714 117.1234
         [,9]    [,10]    [,11]    [,12]
[1,] 105.0538 117.7042 114.9137 117.8939
[2,] 114.9462 121.0958 122.6863 129.9061

$out
[1] 129.5 120.1 119.1 126.5 127.7

$group
[1]  1  1  4  6 10

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

Warning message:
In bxp(list(stats = c(126.3, 126.3, 127, 127.4, 127.4, 103.7, 106.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5vw5a1229416536.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,] -25.8 -11.2  9.3 -12.2 -9.7  8.5 -18.2 -3.7  2.7  -9.2  -2.6 -10.50
[2,] -19.6 -10.2 14.9 -12.2 -4.8  8.5 -16.1 -1.4  3.8  -9.2   1.6  -6.05
[3,] -19.6  -2.1 18.1 -10.2  0.7 12.9 -16.1  0.2  5.2   3.6   3.8   2.85
[4,] -15.2  -1.3 22.3  -9.5  4.6 17.7 -14.3  2.1 20.0   4.6   6.6   8.70
[5,] -15.2  10.5 22.6  -9.5 15.2 21.5 -14.3  3.2 21.3   6.3  12.2  10.10

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

$conf
          [,1]      [,2]     [,3]       [,4]      [,5]      [,6]      [,7]
[1,] -22.70903 -8.388718 12.87118 -12.107813 -5.942016  6.399303 -17.37188
[2,] -16.49097  4.188718 23.32882  -8.292187  7.342016 19.400697 -14.82812
          [,8]      [,9]     [,10]     [,11]   [,12]
[1,] -2.273091 -6.246879 -6.151045 0.2670126 -8.8025
[2,]  2.673091 16.646879 13.351045 7.3329874 14.5025

$out
[1]  -6.9  -5.1 -22.7  -6.9  -3.1

$group
[1] 1 4 4 6 7

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

Warning message:
In bxp(list(stats = c(-25.8, -19.6, -19.6, -15.2, -15.2, -11.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6i8871229416536.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,] 103.70 111.00 100.60  96.5  99.60   NA
[2,] 115.35 112.90 110.70 107.2 104.90   NA
[3,] 119.50 117.15 115.70 114.4 109.75   NA
[4,] 125.25 125.40 120.85 119.6 122.40   NA
[5,] 130.60 128.70 129.60 127.7 127.40   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 114.9845 111.4487 111.0705 108.7443 101.7681   NA
[2,] 124.0155 122.8513 120.3295 120.0557 117.7319   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(103.7, 115.35, 119.5, 125.25, 130.6, 111,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7u5gt1229416536.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,] 105.78 106.50 105.85
[2,] 109.70 110.15 109.40
[3,] 116.02 115.10 116.20
[4,] 123.22 121.55 123.90
[5,] 126.06 127.00 126.85

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 109.8534 109.9004 109.5865
[2,] 122.1866 120.2996 122.8135

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(105.78, 109.7, 116.02, 123.22, 126.06, 106.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1nwzc1229416536.ps tmp/1nwzc1229416536.png")
> system("convert tmp/2wd9h1229416536.ps tmp/2wd9h1229416536.png")
> system("convert tmp/3r0be1229416536.ps tmp/3r0be1229416536.png")
> system("convert tmp/40l0w1229416536.ps tmp/40l0w1229416536.png")
> system("convert tmp/5vw5a1229416536.ps tmp/5vw5a1229416536.png")
> system("convert tmp/6i8871229416536.ps tmp/6i8871229416536.png")
> system("convert tmp/7u5gt1229416536.ps tmp/7u5gt1229416536.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.186   0.934   1.649