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(2.08,2.12,2.14,2.13,2.1,2.09,2.1,2.09,2.08,2.07,2.08,2.09,2.11,2.2,2.42,2.46,2.5,2.59,2.75,2.78,2.9,3.03,3.1,3.23,3.36,3.51,3.61,3.67,3.74,3.82,3.89,3.98,4.08,4.14,4.33,4.57,4.63,4.57,4.71,4.54,4.3,4.36,4.61,4.71,4.68,4.91,4.75,4.77,5.18,3.42,2.71,2.29,2,1.64,1.3,1.08,1,1,1,1)
> 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] 60
> (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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,] 2.08 2.11 3.36 4.63 5.18   NA
 [2,] 2.12 2.20 3.51 4.57 3.42   NA
 [3,] 2.14 2.42 3.61 4.71 2.71   NA
 [4,] 2.13 2.46 3.67 4.54 2.29   NA
 [5,] 2.10 2.50 3.74 4.30 2.00   NA
 [6,] 2.09 2.59 3.82 4.36 1.64   NA
 [7,] 2.10 2.75 3.89 4.61 1.30   NA
 [8,] 2.09 2.78 3.98 4.71 1.08   NA
 [9,] 2.08 2.90 4.08 4.68 1.00   NA
[10,] 2.07 3.03 4.14 4.91 1.00   NA
[11,] 2.08 3.10 4.33 4.75 1.00   NA
[12,] 2.09 3.23 4.57 4.77 1.00   NA
> darr
       [,1] [,2] [,3]  [,4]  [,5] [,6]
 [1,]  0.04 0.09 0.15 -0.06 -1.76   NA
 [2,]  0.02 0.22 0.10  0.14 -0.71   NA
 [3,] -0.01 0.04 0.06 -0.17 -0.42   NA
 [4,] -0.03 0.04 0.07 -0.24 -0.29   NA
 [5,] -0.01 0.09 0.08  0.06 -0.36   NA
 [6,]  0.01 0.16 0.07  0.25 -0.34   NA
 [7,] -0.01 0.03 0.09  0.10 -0.22   NA
 [8,] -0.01 0.12 0.10 -0.03 -0.08   NA
 [9,] -0.01 0.13 0.06  0.23  0.00   NA
[10,]  0.01 0.07 0.19 -0.16  0.00   NA
[11,]  0.01 0.13 0.24  0.02  0.00   NA
[12,]  0.02 0.13 0.06  0.41    NA   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/1nmy51258151968.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/2kole1258151968.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/31ljd1258151968.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/4mtzs1258151968.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,] 2.08 2.12 2.14 2.13 2.00 1.64 1.30 1.08 1.00  1.00  1.00  1.00
[2,] 2.11 2.20 2.42 2.29 2.10 2.09 2.10 2.09 2.08  2.07  2.08  2.09
[3,] 3.36 3.42 2.71 2.46 2.50 2.59 2.75 2.78 2.90  3.03  3.10  3.23
[4,] 4.63 3.51 3.61 3.67 3.74 3.82 3.89 3.98 4.08  4.14  4.33  4.57
[5,] 5.18 4.57 4.71 4.54 4.30 4.36 4.61 4.71 4.68  4.91  4.75  4.77

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 1.579374 2.494357 1.869149 1.484895 1.34118 1.367586 1.485191 1.444531
[2,] 5.140626 4.345643 3.550851 3.435105 3.65882 3.812414 4.014809 4.115469
         [,9]    [,10]    [,11]    [,12]
[1,] 1.486805 1.567343 1.510156 1.477638
[2,] 4.313195 4.492657 4.689844 4.982362

$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(2.08, 2.11, 3.36, 4.63, 5.18, 2.12, 2.2, 3.42,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/52au11258151968.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.06 0.02 -0.42 -0.29 -0.01 0.01 -0.01 -0.08 -0.01  0.00  0.00 0.020
[2,] -0.06 0.02 -0.17 -0.24 -0.01 0.01 -0.01 -0.03  0.00  0.00  0.01 0.040
[3,]  0.04 0.10 -0.01 -0.03  0.06 0.07  0.03 -0.01  0.06  0.01  0.02 0.095
[4,]  0.09 0.14  0.04  0.04  0.08 0.16  0.09  0.10  0.13  0.07  0.13 0.270
[5,]  0.15 0.22  0.06  0.07  0.09 0.25  0.10  0.12  0.23  0.07  0.24 0.410

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

$conf
            [,1]       [,2]       [,3]       [,4]         [,5]        [,6]
[1,] -0.06598962 0.01520830 -0.1583855 -0.2278473 -0.003593773 -0.03598962
[2,]  0.14598962 0.18479170  0.1383855  0.1678473  0.123593773  0.17598962
            [,7]        [,8]        [,9]       [,10]      [,11]   [,12]
[1,] -0.04065975 -0.10185767 -0.03185767 -0.03946182 -0.0647917 -0.0867
[2,]  0.10065975  0.08185767  0.15185767  0.05946182  0.1047917  0.2767

$out
[1] -1.76 -0.71 -0.36 -0.34 -0.22  0.19 -0.16

$group
[1]  1  2  5  6  7 10 10

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

Warning message:
In bxp(list(stats = c(-0.0599999999999996, -0.0599999999999996,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6bz4e1258151968.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,] 2.07 2.110 3.360 4.300 1.00   NA
[2,] 2.08 2.440 3.640 4.555 1.00   NA
[3,] 2.09 2.670 3.855 4.655 1.47   NA
[4,] 2.11 2.965 4.110 4.730 2.50   NA
[5,] 2.14 3.230 4.570 4.910 3.42   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]    [,3]     [,4]    [,5] [,6]
[1,] 2.076317 2.430544 3.64063 4.575181 0.78584   NA
[2,] 2.103683 2.909456 4.06937 4.734819 2.15416   NA

$out
[1] 5.18

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(2.07, 2.08, 2.09, 2.11, 2.14, 2.11, 2.44,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7xhtb1258151968.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,] 2.900 2.460 2.8550
[2,] 2.929 2.650 2.9675
[3,] 3.024 2.840 3.0250
[4,] 3.125 3.165 3.1550
[5,] 3.164 3.420 3.3700

$n
[1] 12 12 12

$conf
         [,1]     [,2]    [,3]
[1,] 2.934603 2.605105 2.93948
[2,] 3.113397 3.074895 3.11052

$out
[1] 3.472

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(2.9, 2.929, 3.024, 3.125, 3.164, 2.46, 2.65,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1nmy51258151968.ps tmp/1nmy51258151968.png")
> system("convert tmp/2kole1258151968.ps tmp/2kole1258151968.png")
> system("convert tmp/31ljd1258151968.ps tmp/31ljd1258151968.png")
> system("convert tmp/4mtzs1258151968.ps tmp/4mtzs1258151968.png")
> system("convert tmp/52au11258151968.ps tmp/52au11258151968.png")
> system("convert tmp/6bz4e1258151968.ps tmp/6bz4e1258151968.png")
> system("convert tmp/7xhtb1258151968.ps tmp/7xhtb1258151968.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.196   0.935   1.589