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(20.3,20,19.2,21.8,21.3,21.5,19.5,19.5,19.7,18.7,19.7,20,19.7,19.2,19.7,22,21.8,22.8,21,25,23.3,25,26.8,25.3,26.5,27.8,22,22.3,28,25,27.3,25.8,27.3,23.5,24.5,18,21.3,21.8,20.5,22.3,18.7,22.3,17.7,19.7,20.5,18.5,10,14.2,15.5,16.5,20.5,15.7,11.7,7.5,3.5,4.5,2.2,5,2.3,6.1,3.3)
> 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,] 20.3 19.7 26.5 21.3 15.5  3.3
 [2,] 20.0 19.2 27.8 21.8 16.5   NA
 [3,] 19.2 19.7 22.0 20.5 20.5   NA
 [4,] 21.8 22.0 22.3 22.3 15.7   NA
 [5,] 21.3 21.8 28.0 18.7 11.7   NA
 [6,] 21.5 22.8 25.0 22.3  7.5   NA
 [7,] 19.5 21.0 27.3 17.7  3.5   NA
 [8,] 19.5 25.0 25.8 19.7  4.5   NA
 [9,] 19.7 23.3 27.3 20.5  2.2   NA
[10,] 18.7 25.0 23.5 18.5  5.0   NA
[11,] 19.7 26.8 24.5 10.0  2.3   NA
[12,] 20.0 25.3 18.0 14.2  6.1   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,] -0.3 -0.5  1.3  0.5  1.0   NA
 [2,] -0.8  0.5 -5.8 -1.3  4.0   NA
 [3,]  2.6  2.3  0.3  1.8 -4.8   NA
 [4,] -0.5 -0.2  5.7 -3.6 -4.0   NA
 [5,]  0.2  1.0 -3.0  3.6 -4.2   NA
 [6,] -2.0 -1.8  2.3 -4.6 -4.0   NA
 [7,]  0.0  4.0 -1.5  2.0  1.0   NA
 [8,]  0.2 -1.7  1.5  0.8 -2.3   NA
 [9,] -1.0  1.7 -3.8 -2.0  2.8   NA
[10,]  1.0  1.8  1.0 -8.5 -2.7   NA
[11,]  0.3 -1.5 -6.5  4.2  3.8   NA
[12,] -0.3  1.2  3.3  1.3 -2.8   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/115vr1258124333.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/22e3f1258124333.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/3vnxi1258124333.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/4ef571258124333.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,] 15.5 16.5 19.2 21.8 18.7 21.5 17.7 19.5 19.7  18.5   2.3   6.1
[2,] 15.5 19.2 19.7 21.8 18.7 21.5 17.7 19.5 19.7  18.5  10.0  14.2
[3,] 20.0 20.0 20.5 22.0 21.3 22.3 19.5 19.7 20.5  18.7  19.7  18.0
[4,] 21.3 21.8 20.5 22.3 21.8 22.8 21.0 25.0 23.3  23.5  24.5  20.0
[5,] 26.5 21.8 20.5 22.3 21.8 22.8 21.0 25.8 27.3  25.0  26.8  25.3

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 16.25881 18.16285 19.93472 21.6467 19.10955 21.38142 17.16823 15.81371
[2,] 23.74119 21.83715 21.06528 22.3533 23.49045 23.21858 21.83177 23.58629
         [,9]    [,10]     [,11]    [,12]
[1,] 17.95625 15.16701  9.454337 13.90173
[2,] 23.04375 22.23299 29.945663 22.09827

$out
 [1]  3.3 27.8 22.0 15.7 28.0 11.7 25.0  7.5 27.3  3.5  4.5  2.2  5.0

$group
 [1]  1  2  3  4  5  5  6  6  7  7  8  9 10

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

Warning message:
In bxp(list(stats = c(15.5, 15.5, 20, 21.3, 26.5, 16.5, 19.2, 20,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5mew31258124333.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,] -0.5 -1.3  0.3 -4.0 -4.2 -4.6 -1.5 -2.3 -3.8  -2.7  -6.5  -0.3
[2,] -0.3 -1.3  0.3 -3.6 -3.0 -4.0  0.0 -1.7 -2.0  -2.7  -1.5  -0.3
[3,]  0.5 -0.8  1.8 -0.5  0.2 -2.0  1.0  0.2 -1.0   1.0   0.3   1.2
[4,]  1.0  0.5  2.3 -0.2  1.0 -1.8  2.0  0.8  1.7   1.0   3.8   1.3
[5,]  1.3  0.5  2.6 -0.2  3.6 -1.8  4.0  1.5  2.8   1.8   4.2   3.3

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

$conf
           [,1]       [,2]     [,3]      [,4]     [,5]       [,6]      [,7]
[1,] -0.4185767 -2.0718755 0.386805 -2.902431 -2.62639 -3.5545145 -0.413195
[2,]  1.4185767  0.4718755 3.213195  1.902431  3.02639 -0.4454855  2.413195
          [,8]      [,9]     [,10]     [,11]      [,12]
[1,] -1.566494 -3.614411 -1.614411 -3.444967 0.06944403
[2,]  1.966494  1.614411  3.614411  4.044967 2.33055597

$out
[1] -5.8  4.0 -4.8  5.7  2.3 -8.5 -2.8

$group
[1]  2  2  3  4  6 10 12

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

Warning message:
In bxp(list(stats = c(-0.5, -0.300000000000001, 0.5, 1, 1.3, -1.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6l32j1258124333.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,] 18.70 19.20 18.0 14.20  2.2  3.3
[2,] 19.50 20.35 22.9 18.10  4.0  3.3
[3,] 19.85 22.40 25.4 20.10  6.8  3.3
[4,] 20.80 25.00 27.3 21.55 15.6  3.3
[5,] 21.80 26.80 28.0 22.30 20.5  3.3

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]      [,5] [,6]
[1,] 19.25706 20.27910 23.39313 18.52643  1.509162  3.3
[2,] 20.44294 24.52090 27.40687 21.67357 12.090838  3.3

$out
[1] 10

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(18.7, 19.5, 19.85, 20.8, 21.8, 19.2, 20.35,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7akxf1258124333.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,] 16.66000 18.0 17.100
[2,] 17.78333 19.6 19.075
[3,] 18.75000 20.0 20.375
[4,] 20.34000 20.9 21.775
[5,] 21.06000 22.3 22.250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 17.58389 19.40706 19.14351
[2,] 19.91611 20.59294 21.60649

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(16.66, 17.7833333333333, 18.75, 20.34, 21.06,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/115vr1258124333.ps tmp/115vr1258124333.png")
> system("convert tmp/22e3f1258124333.ps tmp/22e3f1258124333.png")
> system("convert tmp/3vnxi1258124333.ps tmp/3vnxi1258124333.png")
> system("convert tmp/4ef571258124333.ps tmp/4ef571258124333.png")
> system("convert tmp/5mew31258124333.ps tmp/5mew31258124333.png")
> system("convert tmp/6l32j1258124333.ps tmp/6l32j1258124333.png")
> system("convert tmp/7akxf1258124333.ps tmp/7akxf1258124333.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.191   0.954   1.458