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.7,19.7,18.4,17.4,17,18,23.8,25.5,25.6,23.7,22,21.3,20.7,20.4,20.3,20.4,19.8,19.5,23.1,23.5,23.5,22.9,21.9,21.5,20.5,20.2,19.4,19.2,18.8,18.8,22.6,23.3,23,21.4,19.9,18.8,18.6,18.4,18.6,19.9,19.2,18.4,21.1,20.5,19.1,18.1,17,17.1,17.4,16.8,15.3,14.3,13.4,15.3,22.1,23.7,22.2,19.5,16.6,17.3,19.8,21.2,21.5,20.6,19.1,19.6,23.5,24,23.2)
> par1 = '6'
> #'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] 69
> (np <- floor(n / par1))
[1] 11
> 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] 12 12 12 11 11 11
> arr
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
[1,] 20.7 23.8 20.7 23.1 20.5 22.6 18.6 21.1 17.4  22.1  19.8  23.5
[2,] 19.7 25.5 20.4 23.5 20.2 23.3 18.4 20.5 16.8  23.7  21.2  24.0
[3,] 18.4 25.6 20.3 23.5 19.4 23.0 18.6 19.1 15.3  22.2  21.5  23.2
[4,] 17.4 23.7 20.4 22.9 19.2 21.4 19.9 18.1 14.3  19.5  20.6    NA
[5,] 17.0 22.0 19.8 21.9 18.8 19.9 19.2 17.0 13.4  16.6  19.1    NA
[6,] 18.0 21.3 19.5 21.5 18.8 18.8 18.4 17.1 15.3  17.3  19.6    NA
> darr
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
[1,] -1.0  1.7 -0.3  0.4 -0.3  0.7 -0.2 -0.6 -0.6   1.6   1.4   0.5
[2,] -1.3  0.1 -0.1  0.0 -0.8 -0.3  0.2 -1.4 -1.5  -1.5   0.3  -0.8
[3,] -1.0 -1.9  0.1 -0.6 -0.2 -1.6  1.3 -1.0 -1.0  -2.7  -0.9    NA
[4,] -0.4 -1.7 -0.6 -1.0 -0.4 -1.5 -0.7 -1.1 -0.9  -2.9  -1.5    NA
[5,]  1.0 -0.7 -0.3 -0.4  0.0 -1.1 -0.8  0.1  1.9   0.7   0.5    NA
[6,]  5.8 -0.6  3.6 -1.0  3.8 -0.2  2.7  0.3  6.8   2.5   3.9    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/13bk41259247979.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/2x99y1259247979.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/33xg11259247979.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/494io1259247979.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]
[1,] 17.40 16.80 15.30 17.40 13.40 15.30
[2,] 20.15 19.95 18.85 18.65 17.00 17.65
[3,] 20.90 20.85 20.90 19.90 19.10 18.80
[4,] 22.85 23.60 23.10 21.00 19.85 19.55
[5,] 23.80 25.50 25.60 23.70 22.00 21.50

$n
[1] 12 12 12 11 11 11

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]
[1,] 19.66851 19.18521 18.96155 18.78049 17.74229 17.89486
[2,] 22.13149 22.51479 22.83845 21.01951 20.45771 19.70514

$out
[1] 14.3

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(17.4, 20.15, 20.9, 22.85, 23.8, 16.8, 19.95,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/578md1259247979.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]
[1,] -1.00 -1.50 -1.9 -1.70 -1.10 -1.00
[2,] -0.45 -1.35 -1.3 -1.50 -0.55  0.05
[3,]  0.10 -0.55 -1.0 -1.00  0.00  2.70
[4,]  1.05  0.05 -0.4 -0.65  0.60  3.85
[5,]  1.70  0.30  0.1 -0.40  1.90  6.80

$n
[1] 12 12 11 11 11 11

$conf
           [,1]       [,2]       [,3]       [,4]       [,5]      [,6]
[1,] -0.5841601 -1.1885494 -1.4287491 -1.4049297 -0.5478461 0.8897259
[2,]  0.7841601  0.0885494 -0.5712509 -0.5950703  0.5478461 4.5102741

$out
[1]  1.3 -2.7 -2.9

$group
[1] 3 3 4

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

Warning message:
In bxp(list(stats = c(-1, -0.449999999999999, 0.0999999999999979,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/61d0v1259247979.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] [,11] [,12]
[1,] 17.0 21.30 19.50 21.5 18.8 18.8 18.4 17.0 13.4  16.6  19.1 23.20
[2,] 17.4 22.00 19.80 21.9 18.8 19.9 18.4 17.1 14.3  17.3  19.6 23.35
[3,] 18.2 23.75 20.35 23.0 19.3 22.0 18.6 18.6 15.3  20.8  20.2 23.50
[4,] 19.7 25.50 20.40 23.5 20.2 23.0 19.2 20.5 16.8  22.2  21.2 23.75
[5,] 20.7 25.60 20.70 23.5 20.5 23.3 19.9 21.1 17.4  23.7  21.5 24.00

$n
 [1] 6 6 6 6 6 6 6 6 6 6 6 3

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]     [,7]     [,8]
[1,] 16.71643 21.49239 19.96298 21.96795 18.39695 20.0004 18.08397 16.40689
[2,] 19.68357 26.00761 20.73702 24.03205 20.20305 23.9996 19.11603 20.79311
         [,9]    [,10]    [,11]    [,12]
[1,] 13.68742 17.63934 19.16795 23.13511
[2,] 16.91258 23.96066 21.23205 23.86489

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(17, 17.4, 18.2, 19.7, 20.7, 21.3, 22, 23.75,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7lejt1259247979.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,] 18.60909 18.800 18.42500
[2,] 18.69091 19.100 18.60000
[3,] 20.30265 20.375 20.41875
[4,] 21.15833 20.900 21.52500
[5,] 21.43333 20.900 21.81250

$n
[1] 6 6 6

$conf
         [,1]     [,2]     [,3]
[1,] 18.71108 19.21394 18.53203
[2,] 21.89422 21.53606 22.30547

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(18.6090909090909, 18.6909090909091, 20.3026515151515,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/13bk41259247979.ps tmp/13bk41259247979.png")
> system("convert tmp/2x99y1259247979.ps tmp/2x99y1259247979.png")
> system("convert tmp/33xg11259247979.ps tmp/33xg11259247979.png")
> system("convert tmp/494io1259247979.ps tmp/494io1259247979.png")
> system("convert tmp/578md1259247979.ps tmp/578md1259247979.png")
> system("convert tmp/61d0v1259247979.ps tmp/61d0v1259247979.png")
> system("convert tmp/7lejt1259247979.ps tmp/7lejt1259247979.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.166   0.922   1.473