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(109.6,103,111.6,106.3,97.9,108.8,103.9,101.2,122.9,123.9,111.7,120.9,99.6,103.3,119.4,106.5,101.9,124.6,106.5,107.8,127.4,120.1,118.5,127.7,107.7,104.5,118.8,110.3,109.6,119.1,96.5,106.7,126.3,116.2,118.8,115.2,110,111.4,129.6,108.1,117.8,122.9,100.6,111.8,127,128.6,124.8,118.5,114.7,112.6,128.7,111,115.8,126,111.1,113.2,120.1,130.6,124,119.4,116.7)
> 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] 61
> (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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]
 [1,] 109.6  99.6 107.7 110.0 114.7 116.7
 [2,] 103.0 103.3 104.5 111.4 112.6    NA
 [3,] 111.6 119.4 118.8 129.6 128.7    NA
 [4,] 106.3 106.5 110.3 108.1 111.0    NA
 [5,]  97.9 101.9 109.6 117.8 115.8    NA
 [6,] 108.8 124.6 119.1 122.9 126.0    NA
 [7,] 103.9 106.5  96.5 100.6 111.1    NA
 [8,] 101.2 107.8 106.7 111.8 113.2    NA
 [9,] 122.9 127.4 126.3 127.0 120.1    NA
[10,] 123.9 120.1 116.2 128.6 130.6    NA
[11,] 111.7 118.5 118.8 124.8 124.0    NA
[12,] 120.9 127.7 115.2 118.5 119.4    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,]  -6.6   3.7  -3.2   1.4  -2.1   NA
 [2,]   8.6  16.1  14.3  18.2  16.1   NA
 [3,]  -5.3 -12.9  -8.5 -21.5 -17.7   NA
 [4,]  -8.4  -4.6  -0.7   9.7   4.8   NA
 [5,]  10.9  22.7   9.5   5.1  10.2   NA
 [6,]  -4.9 -18.1 -22.6 -22.3 -14.9   NA
 [7,]  -2.7   1.3  10.2  11.2   2.1   NA
 [8,]  21.7  19.6  19.6  15.2   6.9   NA
 [9,]   1.0  -7.3 -10.1   1.6  10.5   NA
[10,] -12.2  -1.6   2.6  -3.8  -6.6   NA
[11,]   9.2   9.2  -3.6  -6.3  -4.6   NA
[12,] -21.3 -20.0  -5.2  -3.8  -2.7   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/1vdgr1225742004.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/21x371225742004.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/31bmc1225742004.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/4wg9w1225742004.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,]  99.6 103.0 111.6 106.3  97.9 119.1  96.5 101.2 120.1 116.2 111.7 115.2
[2,] 107.7 103.3 118.8 106.5 101.9 119.1 100.6 106.7 122.9 120.1 118.5 118.5
[3,] 109.8 104.5 119.4 108.1 109.6 122.9 103.9 107.8 126.3 123.9 118.8 119.4
[4,] 114.7 111.4 128.7 110.3 115.8 124.6 106.5 111.8 127.0 128.6 124.0 120.9
[5,] 116.7 112.6 129.6 111.0 117.8 126.0 111.1 113.2 127.4 130.6 124.8 120.9

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

$conf
         [,1]      [,2]     [,3]     [,4]     [,5]     [,6]      [,7]     [,8]
[1,] 105.2848  98.77656 112.4047 105.4149  99.7783 119.0137  99.73107 104.1964
[2,] 114.3152 110.22344 126.3953 110.7851 119.4217 126.7863 108.06893 111.4036
         [,9]    [,10]    [,11]    [,12]
[1,] 123.4030 117.8939 114.9137 117.7042
[2,] 129.1970 129.9061 122.6863 121.0958

$out
[1] 108.8 127.7

$group
[1]  6 12

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

Warning message:
In bxp(list(stats = c(99.6, 107.7, 109.8, 114.7, 116.7, 103, 103.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/578ss1225742004.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,] -6.6 14.3 -21.5 -8.4  9.5 -22.6 -2.7 15.2 -10.1 -12.2  -6.3 -21.3
[2,] -3.2 14.3 -17.7 -4.6  9.5 -22.3  1.3 15.2  -7.3  -6.6  -4.6 -20.0
[3,] -2.1 16.1 -12.9 -0.7 10.2 -18.1  2.1 19.6   1.0  -3.8  -3.6  -5.2
[4,]  1.4 16.1  -8.5  4.8 10.9 -14.9 10.2 19.6   1.6  -1.6   9.2  -3.8
[5,]  3.7 18.2  -5.3  9.7 10.9  -4.9 11.2 21.7  10.5   2.6   9.2  -2.7

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

$conf
          [,1]     [,2]       [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -5.350348 14.82812 -19.400697 -7.342016  9.210764 -23.32882 -4.188718
[2,]  1.150348 17.37188  -6.399303  5.942016 11.189236 -12.87118  8.388718
         [,8]      [,9]      [,10]      [,11]      [,12]
[1,] 16.49097 -5.288718 -7.3329874 -13.351045 -16.646879
[2,] 22.70903  7.288718 -0.2670126   6.151045   6.246879

$out
[1]  8.6 22.7  5.1  6.9

$group
[1] 2 5 5 8

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

Warning message:
In bxp(list(stats = c(-6.6, -3.2, -2.10000000000001, 1.40000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6rfva1225742004.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,]  97.90  99.60  96.50 100.60 111.0 116.7
[2,] 103.45 104.90 107.20 110.70 112.9 116.7
[3,] 109.20 113.15 112.75 118.15 117.6 116.7
[4,] 116.30 122.35 118.80 125.90 125.0 116.7
[5,] 123.90 127.70 126.30 129.60 130.6 116.7

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]  [,6]
[1,] 103.3390 105.1909 107.4592 111.2172 112.0811 116.7
[2,] 115.0610 121.1091 118.0408 125.0828 123.1189 116.7

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(97.9, 103.45, 109.2, 116.3, 123.9, 99.6, 104.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/73axu1225742004.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,] 103.7200 103.90 103.550
[2,] 108.2900 107.95 108.625
[3,] 114.6383 114.30 115.275
[4,] 120.9800 121.15 122.800
[5,] 124.7400 126.30 124.950

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 108.8503 108.2794 108.8097
[2,] 120.4263 120.3206 121.7403

$out
numeric(0)

$group
numeric(0)

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

> dev.off()
null device 
          1 
> 
> system("convert tmp/1vdgr1225742004.ps tmp/1vdgr1225742004.png")
> system("convert tmp/21x371225742004.ps tmp/21x371225742004.png")
> system("convert tmp/31bmc1225742004.ps tmp/31bmc1225742004.png")
> system("convert tmp/4wg9w1225742004.ps tmp/4wg9w1225742004.png")
> system("convert tmp/578ss1225742004.ps tmp/578ss1225742004.png")
> system("convert tmp/6rfva1225742004.ps tmp/6rfva1225742004.png")
> system("convert tmp/73axu1225742004.ps tmp/73axu1225742004.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.449   0.987   1.751