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(13,14,8,10,16,26,32,27,18,32,36,24,22,23,22,18,25,21,21,14,8,11,14,23,18,17,19,20,22,19,13,26,13,14,22,24,21,22,26,21,23,24,27,41,31,27,35,26,28,36,39,21,17,22,17,19,15,34,10,15,22,18,15,20,15,22,19,16,30,27,29,23,20,16,21,21,25,16,18,15,18,14,10,15,8,15,6,11,8,7,13,10,23,16,15,25,22,20,16)
> par1 = '8'
> #'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] 99
> (np <- floor(n / par1))
[1] 12
> 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] 13 13 13 12 12 12 12 12
> arr
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,]   13   18   25   18   13   23   28   15   15    20    18     8    22
[2,]   14   32   21   17   14   24   36   34   22    16    14     7    20
[3,]    8   36   21   19   22   27   39   10   19    21    10    13    16
[4,]   10   24   14   20   24   41   21   15   16    21    15    10    NA
[5,]   16   22    8   22   21   31   17   22   30    25     8    23    NA
[6,]   26   23   11   19   22   27   22   18   27    16    15    16    NA
[7,]   32   22   14   13   26   35   17   15   29    18     6    15    NA
[8,]   27   18   23   26   21   26   19   20   23    15    11    25    NA
> darr
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,]    1   14   -4   -1    1    1    8   19    7    -4    -4    -1    -2
[2,]   -6    4    0    2    8    3    3  -24   -3     5    -4     6    -4
[3,]    2  -12   -7    1    2   14  -18    5   -3     0     5    -3    NA
[4,]    6   -2   -6    2   -3  -10   -4    7   14     4    -7    13    NA
[5,]   10    1    3   -3    1   -4    5   -4   -3    -9     7    -7    NA
[6,]    6   -1    3   -6    4    8   -5   -3    2     2    -9    -1    NA
[7,]   -5   -4    9   13   -5   -9    2    5   -6    -3     5    10    NA
[8,]   -9    7   -5  -13    2    2   -4   -5   -3     3    -3    -3    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/freestat/rcomp/tmp/1hbgy1272309347.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/freestat/rcomp/tmp/2hbgy1272309347.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/freestat/rcomp/tmp/39kx11272309347.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/freestat/rcomp/tmp/49kx11272309347.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]
[1,]    8    7    8 10.0  8.0 11.0  6.0 11.0
[2,]   15   14   13 14.5 16.5 16.0 14.5 18.5
[3,]   18   20   19 18.0 22.0 20.5 17.5 22.0
[4,]   22   24   22 22.5 24.0 24.5 27.5 25.5
[5,]   28   36   27 24.0 31.0 27.0 35.0 27.0

$n
[1] 13 13 13 12 12 12 12 12

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 14.93251 15.61787 15.05608 14.35115 18.5792 16.62309 11.57061 18.80725
[2,] 21.06749 24.38213 22.94392 21.64885 25.4208 24.37691 23.42939 25.19275

$out
[1] 36 39 41

$group
[1] 3 3 4

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

Warning message:
In bxp(list(stats = c(8, 15, 18, 22, 28, 7, 14, 20, 24, 36, 8, 13,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/59kx11272309347.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]
[1,]   -4   -6 -12.0 -10.0   -9 -9.0 -9.0  -13
[2,]   -2   -4  -5.0  -5.0   -4 -4.0 -5.0   -5
[3,]    1    2   0.5   0.0   -1  0.5 -0.5   -3
[4,]    7    4   3.5   6.5    4  3.5  7.0    2
[5,]   19    8  14.0  14.0   10  8.0 13.0    7

$n
[1] 13 13 12 12 12 12 12 12

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -2.943918 -1.505705 -3.376907 -5.245227 -4.648854 -2.920800 -5.973281
[2,]  4.943918  5.505705  4.376907  5.245227  2.648854  3.920800  4.973281
          [,8]
[1,] -6.192747
[2,]  0.192747

$out
[1] -24 -18

$group
[1] 2 3

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

Warning message:
In bxp(list(stats = c(-4, -2, 1, 7, 19, -6, -4, 2, 4, 8, -12, -5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/62uxm1272309347.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] [,13]
[1,]  8.0 18.0  8.0 13.0 13.0   23 17.0 10.0 15.0    15   6.0   7.0    16
[2,] 11.5 20.0 12.5 17.5 17.5   25 18.0 15.0 17.5    16   9.0   9.0    18
[3,] 15.0 22.5 17.5 19.0 21.5   27 21.5 16.5 22.5    19  12.5  14.0    20
[4,] 26.5 28.0 22.0 21.0 23.0   33 32.0 21.0 28.0    21  15.0  19.5    21
[5,] 32.0 36.0 25.0 26.0 26.0   41 39.0 22.0 30.0    25  18.0  25.0    22

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 8 3

$conf
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,]  6.620785 18.03109 12.19316 17.04485 18.42762 22.53109 13.6794 13.14831
[2,] 23.379215 26.96891 22.80684 20.95515 24.57238 31.46891 29.3206 19.85169
         [,9]    [,10]     [,11]    [,12]    [,13]
[1,] 16.63455 16.20693  9.148314  8.13455 17.26336
[2,] 28.36545 21.79307 15.851686 19.86545 22.73664

$out
[1] 34

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(8, 11.5, 15, 26.5, 32, 18, 20, 22.5, 28, 36,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/72uxm1272309347.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,] 19.25000 17.50 17.5000
[2,] 19.66346 18.00 18.3750
[3,] 20.16667 19.50 19.4375
[4,] 20.63141 21.25 20.4375
[5,] 21.16667 22.00 22.0000

$n
[1] 8 8 8

$conf
         [,1]     [,2]     [,3]
[1,] 19.62596 17.68450 18.28536
[2,] 20.70738 21.31550 20.58964

$out
[1] 18.15385

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(19.25, 19.6634615384615, 20.1666666666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1hbgy1272309347.ps tmp/1hbgy1272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/2hbgy1272309347.ps tmp/2hbgy1272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/39kx11272309347.ps tmp/39kx11272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/49kx11272309347.ps tmp/49kx11272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/59kx11272309347.ps tmp/59kx11272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/62uxm1272309347.ps tmp/62uxm1272309347.png",intern=TRUE))
character(0)
> try(system("convert tmp/72uxm1272309347.ps tmp/72uxm1272309347.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.897   1.491   2.212