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(107.11,107.57,107.81,108.75,109.43,109.62,109.54,109.53,109.84,109.67,109.79,109.56,110.22,110.40,110.69,110.72,110.89,110.58,110.94,110.91,111.22,111.09,111.00,111.06,111.55,112.32,112.64,112.36,112.04,112.37,112.59,112.89,113.22,112.85,113.06,112.99,113.32,113.74,113.91,114.52,114.96,114.91,115.30,115.44,115.52,116.08,115.94,115.56,115.88,116.66,117.41,117.68,117.85,118.21,118.92,119.03,119.17,118.95,118.92,118.90)
> 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,] 107.11 110.22 111.55 113.32 115.88   NA
 [2,] 107.57 110.40 112.32 113.74 116.66   NA
 [3,] 107.81 110.69 112.64 113.91 117.41   NA
 [4,] 108.75 110.72 112.36 114.52 117.68   NA
 [5,] 109.43 110.89 112.04 114.96 117.85   NA
 [6,] 109.62 110.58 112.37 114.91 118.21   NA
 [7,] 109.54 110.94 112.59 115.30 118.92   NA
 [8,] 109.53 110.91 112.89 115.44 119.03   NA
 [9,] 109.84 111.22 113.22 115.52 119.17   NA
[10,] 109.67 111.09 112.85 116.08 118.95   NA
[11,] 109.79 111.00 113.06 115.94 118.92   NA
[12,] 109.56 111.06 112.99 115.56 118.90   NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,]  0.46  0.18  0.77  0.42  0.78   NA
 [2,]  0.24  0.29  0.32  0.17  0.75   NA
 [3,]  0.94  0.03 -0.28  0.61  0.27   NA
 [4,]  0.68  0.17 -0.32  0.44  0.17   NA
 [5,]  0.19 -0.31  0.33 -0.05  0.36   NA
 [6,] -0.08  0.36  0.22  0.39  0.71   NA
 [7,] -0.01 -0.03  0.30  0.14  0.11   NA
 [8,]  0.31  0.31  0.33  0.08  0.14   NA
 [9,] -0.17 -0.13 -0.37  0.56 -0.22   NA
[10,]  0.12 -0.09  0.21 -0.14 -0.03   NA
[11,] -0.23  0.06 -0.07 -0.38 -0.02   NA
[12,]  0.66  0.49  0.33  0.32    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/1ulkw1258113027.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/23e8i1258113027.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/3smue1258113027.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/4twmc1258113027.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]
[1,] 107.11 107.57 107.81 108.75 109.43 109.62 109.54 109.53 109.84 109.67
[2,] 110.22 110.40 110.69 110.72 110.89 110.58 110.94 110.91 111.22 111.09
[3,] 111.55 112.32 112.64 112.36 112.04 112.37 112.59 112.89 113.22 112.85
[4,] 113.32 113.74 113.91 114.52 114.96 114.91 115.30 115.44 115.52 116.08
[5,] 115.88 116.66 117.41 117.68 117.85 118.21 118.92 119.03 119.17 118.95
      [,11]  [,12]
[1,] 109.79 109.56
[2,] 111.00 111.06
[3,] 113.06 112.99
[4,] 115.94 115.56
[5,] 118.92 118.90

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 109.3595 109.9600 110.3648 109.6749 109.1641 109.3104 109.5092 109.6891
[2,] 113.7405 114.6800 114.9152 115.0451 114.9159 115.4296 115.6708 116.0909
         [,9]    [,10]    [,11]    [,12]
[1,] 110.1816 109.3241 109.5694 109.8103
[2,] 116.2584 116.3759 116.5506 116.1697

$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(107.11, 110.22, 111.55, 113.32, 115.88, 107.57,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5bmd51258113027.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.18 0.17 -0.28 0.17 -0.31 0.22 -0.03 0.08 -0.22 -0.14 -0.38 0.320
[2,] 0.42 0.24  0.03 0.17 -0.05 0.22 -0.01 0.14 -0.22 -0.09 -0.23 0.325
[3,] 0.46 0.29  0.27 0.17  0.19 0.36  0.11 0.31 -0.17 -0.03 -0.07 0.410
[4,] 0.77 0.32  0.61 0.44  0.33 0.39  0.14 0.31 -0.13  0.12 -0.02 0.575
[5,] 0.78 0.32  0.94 0.68  0.36 0.39  0.30 0.33 -0.13  0.21  0.06 0.660

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

$conf
          [,1]      [,2]       [,3]        [,4]        [,5]      [,6]
[1,] 0.2126909 0.2334722 -0.1398265 -0.02078132 -0.07850704 0.2398784
[2,] 0.7073091 0.3465278  0.6798265  0.36078132  0.45850704 0.4801216
            [,7]      [,8]       [,9]      [,10]       [,11]  [,12]
[1,] 0.004010378 0.1898784 -0.2335938 -0.1783855 -0.21838547 0.2125
[2,] 0.215989622 0.4301216 -0.1064062  0.1183855  0.07838547 0.6075

$out
[1]  0.75 -0.32 -0.08  0.71 -0.37  0.56

$group
[1] 2 4 6 6 9 9

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

Warning message:
In bxp(list(stats = c(0.180000000000007, 0.420000000000002, 0.459999999999994,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6f3fx1258113027.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,] 107.110 110.220 111.550 113.320 115.880   NA
[2,] 108.280 110.635 112.340 114.215 117.545   NA
[3,] 109.535 110.900 112.615 115.130 118.555   NA
[4,] 109.645 111.030 112.940 115.540 118.935   NA
[5,] 109.840 111.220 113.220 116.080 119.170   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]    [,5] [,6]
[1,] 108.9124 110.7198 112.3413 114.5257 117.921   NA
[2,] 110.1576 111.0802 112.8887 115.7343 119.189   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(107.11, 108.28, 109.535, 109.645, 109.84,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7oami1258113027.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,] 111.616 111.550 111.7700
[2,] 112.649 112.340 112.4600
[3,] 113.298 112.615 113.0225
[4,] 113.671 112.940 113.3400
[5,] 113.794 113.220 113.5850

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 112.8319 112.3413 112.6211
[2,] 113.7641 112.8887 113.4239

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(111.616, 112.649, 113.298, 113.671, 113.794,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ulkw1258113027.ps tmp/1ulkw1258113027.png")
> system("convert tmp/23e8i1258113027.ps tmp/23e8i1258113027.png")
> system("convert tmp/3smue1258113027.ps tmp/3smue1258113027.png")
> system("convert tmp/4twmc1258113027.ps tmp/4twmc1258113027.png")
> system("convert tmp/5bmd51258113027.ps tmp/5bmd51258113027.png")
> system("convert tmp/6f3fx1258113027.ps tmp/6f3fx1258113027.png")
> system("convert tmp/7oami1258113027.ps tmp/7oami1258113027.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.215   0.948   3.041