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(13.2,13.8,16.2,14.7,13.9,16.0,14.4,12.3,15.9,15.9,15.5,15.1,14.5,15.1,17.4,16.2,15.6,17.2,14.9,13.8,17.5,16.2,17.5,16.6,16.2,16.6,19.6,15.9,18.0,18.3,16.3,14.9,18.2,18.4,18.5,16.0,17.4,17.2,19.6,17.2,18.3,19.3,18.1,16.2,18.4,20.5,19.0,16.5,18.7,19.0,19.2,20.5,19.3,20.6,20.1,16.1,20.4,19.7,15.6,14.4,13.9,14.3,15.3,14.4,13.8,15.7,14.7,12.5,16.2,16.1,16,15.8,15.2)
> 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] 73
> (np <- floor(n / par1))
[1] 6
> 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] 7 6 6 6 6 6 6 6 6 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,] 13.2 14.5 16.2 17.4 18.7 13.9 15.2
 [2,] 13.8 15.1 16.6 17.2 19.0 14.3   NA
 [3,] 16.2 17.4 19.6 19.6 19.2 15.3   NA
 [4,] 14.7 16.2 15.9 17.2 20.5 14.4   NA
 [5,] 13.9 15.6 18.0 18.3 19.3 13.8   NA
 [6,] 16.0 17.2 18.3 19.3 20.6 15.7   NA
 [7,] 14.4 14.9 16.3 18.1 20.1 14.7   NA
 [8,] 12.3 13.8 14.9 16.2 16.1 12.5   NA
 [9,] 15.9 17.5 18.2 18.4 20.4 16.2   NA
[10,] 15.9 16.2 18.4 20.5 19.7 16.1   NA
[11,] 15.5 17.5 18.5 19.0 15.6 16.0   NA
[12,] 15.1 16.6 16.0 16.5 14.4 15.8   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]  0.6  0.6  0.4 -0.2  0.3  0.4   NA
 [2,]  2.4  2.3  3.0  2.4  0.2  1.0   NA
 [3,] -1.5 -1.2 -3.7 -2.4  1.3 -0.9   NA
 [4,] -0.8 -0.6  2.1  1.1 -1.2 -0.6   NA
 [5,]  2.1  1.6  0.3  1.0  1.3  1.9   NA
 [6,] -1.6 -2.3 -2.0 -1.2 -0.5 -1.0   NA
 [7,] -2.1 -1.1 -1.4 -1.9 -4.0 -2.2   NA
 [8,]  3.6  3.7  3.3  2.2  4.3  3.7   NA
 [9,]  0.0 -1.3  0.2  2.1 -0.7 -0.1   NA
[10,] -0.4  1.3  0.1 -1.5 -4.1 -0.1   NA
[11,] -0.4 -0.9 -2.5 -2.5 -1.2 -0.2   NA
[12,] -0.6 -0.4  1.4  2.2 -0.5 -0.6   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/13oe31289994301.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/23oe31289994301.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/3wxdn1289994301.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/4wxdn1289994301.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,] 13.2 13.80 15.3 14.40 13.8 15.70 14.4 12.30 15.90  15.9 15.50  14.4
[2,] 14.2 14.30 16.2 14.70 13.9 16.00 14.7 12.50 16.20  16.1 15.60  15.1
[3,] 15.2 15.85 18.3 16.05 16.8 17.75 15.6 14.35 17.85  17.3 16.75  15.9
[4,] 16.8 17.20 19.6 17.20 18.3 19.30 18.1 16.10 18.40  19.7 18.50  16.5
[5,] 18.7 19.00 19.6 20.50 19.3 20.60 20.1 16.20 20.40  20.5 19.00  16.6

$n
 [1] 7 6 6 6 6 6 6 6 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 13.64732 13.97941 16.10689 14.43742 13.96186 15.62139 13.40689 12.02788
[2,] 16.75268 17.72059 20.49311 17.66258 19.63814 19.87861 17.79311 16.67212
         [,9]    [,10]    [,11]    [,12]
[1,] 16.43093 14.97788 14.87941 14.99695
[2,] 19.26907 19.62212 18.62059 16.80305

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(13.2, 14.2, 15.2, 16.8, 18.7, 13.8, 14.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5podr1289994301.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.3 0.20 -3.70 -1.2 0.30 -2.3 -2.2 3.30 -1.30 -1.50 -2.50 -0.60
[2,]  0.3 1.00 -2.40 -0.8 1.00 -2.0 -2.2 3.30 -0.70 -1.50 -2.50 -0.60
[3,]  0.4 2.35 -1.35 -0.6 1.45 -1.4 -2.0 3.65 -0.05 -0.25 -1.05 -0.45
[4,]  0.6 2.40 -0.90  1.1 1.90 -1.0 -1.4 3.70  0.20  0.10 -0.40  1.40
[5,]  0.6 3.00  1.30  2.1 2.10 -0.5 -1.1 4.30  0.20  1.30 -0.20  2.20

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

$conf
          [,1]     [,2]       [,3]       [,4]     [,5]       [,6]      [,7]
[1,] 0.2064903 1.446955 -2.3175484 -1.8255614 0.869471 -2.0450323 -2.516026
[2,] 0.5935097 3.253045 -0.3824516  0.6255614 2.030529 -0.7549677 -1.483974
         [,8]       [,9]      [,10]      [,11]      [,12]
[1,] 3.391987 -0.6305291 -1.2820517 -2.4045678 -1.7400646
[2,] 3.908013  0.5305291  0.7820517  0.3045678  0.8400646

$out
[1] -0.2 -4.0  2.2  2.1 -4.1

$group
[1]  1  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(0.300000000000001, 0.300000000000001, 0.400000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6podr1289994301.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]
[1,] 12.30 13.8 14.90 16.20 14.40 12.5 15.2
[2,] 13.85 15.0 16.10 17.20 17.40 14.1 15.2
[3,] 14.90 16.2 17.30 18.20 19.25 15.0 15.2
[4,] 15.90 17.3 18.35 19.15 20.25 15.9 15.2
[5,] 16.20 17.5 19.60 20.50 20.60 16.2 15.2

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 13.96498 15.15095 16.27376 17.31059 17.95010 14.17901 15.2
[2,] 15.83502 17.24905 18.32624 19.08941 20.54990 15.82099 15.2

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(12.3, 13.85, 14.9, 15.9, 16.2, 13.8, 15, 16.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7podr1289994301.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,] 14.30000 14.350 14.31250
[2,] 15.86667 15.725 15.80000
[3,] 16.48333 16.400 16.23750
[4,] 17.78333 17.525 17.55625
[5,] 17.88333 18.300 18.00000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 15.60913 15.57901 15.43646
[2,] 17.35754 17.22099 17.03854

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(14.3, 15.8666666666667, 16.4833333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/13oe31289994301.ps tmp/13oe31289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/23oe31289994301.ps tmp/23oe31289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/3wxdn1289994301.ps tmp/3wxdn1289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/4wxdn1289994301.ps tmp/4wxdn1289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/5podr1289994301.ps tmp/5podr1289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/6podr1289994301.ps tmp/6podr1289994301.png",intern=TRUE))
character(0)
> try(system("convert tmp/7podr1289994301.ps tmp/7podr1289994301.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.199   0.954   2.915