R version 2.11.1 (2010-05-31)
Copyright (C) 2010 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(102.89,102.64,103.33,103.56,103.60,104.24,105.31,105.40,105.89,105.89,105.54,106.15,106.14,105.85,106.27,106.51,106.82,106.53,107.14,107.39,107.33,107.53,107.42,108.25,108.26,108.93,109.43,109.61,109.74,110.12,110.16,110.44,111.23,112.86,112.77,113.04,112.79,113.87,114.28,115.51,116.76,116.91,116.47,116.94,117.24,116.82,117.48,117.11,117.31,117.77,118.37,117.91,118.12,118.02,117.77,117.85,118.68,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] 58
> (np <- floor(n / par1))
[1] 4
> 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 4 4
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]
 [1,] 102.89 106.14 108.26 112.79 117.31
 [2,] 102.64 105.85 108.93 113.87 117.77
 [3,] 103.33 106.27 109.43 114.28 118.37
 [4,] 103.56 106.51 109.61 115.51 117.91
 [5,] 103.60 106.82 109.74 116.76 118.12
 [6,] 104.24 106.53 110.12 116.91 118.02
 [7,] 105.31 107.14 110.16 116.47 117.77
 [8,] 105.40 107.39 110.44 116.94 117.85
 [9,] 105.89 107.33 111.23 117.24 118.68
[10,] 105.89 107.53 112.86 116.82 118.90
[11,] 105.54 107.42 112.77 117.48     NA
[12,] 106.15 108.25 113.04 117.11     NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] -0.25 -0.29  0.67  1.08  0.46
 [2,]  0.69  0.42  0.50  0.41  0.60
 [3,]  0.23  0.24  0.18  1.23 -0.46
 [4,]  0.04  0.31  0.13  1.25  0.21
 [5,]  0.64 -0.29  0.38  0.15 -0.10
 [6,]  1.07  0.61  0.04 -0.44 -0.25
 [7,]  0.09  0.25  0.28  0.47  0.08
 [8,]  0.49 -0.06  0.79  0.30  0.83
 [9,]  0.00  0.20  1.63 -0.42  0.22
[10,] -0.35 -0.11 -0.09  0.66    NA
[11,]  0.61  0.83  0.27 -0.37    NA
[12,] -0.01  0.01 -0.25  0.20    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/rcomp/tmp/12mo61289573154.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/rcomp/tmp/22mo61289573154.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/rcomp/tmp/32mo61289573154.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/rcomp/tmp/4ve6r1289573154.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,] 102.89 102.64 103.33 103.56 103.60 104.24 105.31 105.40 105.89 105.89
[2,] 106.14 105.85 106.27 106.51 106.82 106.53 107.14 107.39 107.33 107.53
[3,] 108.26 108.93 109.43 109.61 109.74 110.12 110.16 110.44 111.23 112.86
[4,] 112.79 113.87 114.28 115.51 116.76 116.91 116.47 116.94 117.24 116.82
[5,] 117.31 117.77 118.37 117.91 118.12 118.02 117.77 117.85 118.68 118.90
       [,11]   [,12]
[1,] 105.540 106.150
[2,] 106.480 107.200
[3,] 110.095 110.645
[4,] 115.125 115.075
[5,] 117.480 117.110

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 103.5611 103.2631 103.7702 103.2506 102.7164 102.7855 103.5674 103.692
[2,] 112.9589 114.5969 115.0898 115.9694 116.7636 117.4545 116.7526 117.188
         [,9]    [,10]    [,11]    [,12]
[1,] 104.2276 106.2957 103.2655 104.4238
[2,] 118.2324 119.4243 116.9245 116.8663

$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(102.89, 106.14, 108.26, 112.79, 117.31, 102.64,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5ve6r1289573154.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.29 0.41 0.18 0.04 -0.29 -0.44 0.08 -0.06 0.00 -0.350 -0.37 -0.250
[2,] -0.25 0.42 0.18 0.13 -0.10 -0.25 0.09  0.30 0.00 -0.230 -0.05 -0.130
[3,]  0.46 0.50 0.23 0.21  0.15  0.04 0.25  0.49 0.20 -0.100  0.44  0.000
[4,]  0.67 0.60 0.24 0.31  0.38  0.61 0.28  0.79 0.22  0.285  0.72  0.105
[5,]  1.08 0.69 0.24 0.31  0.64  1.07 0.47  0.83 0.22  0.660  0.83  0.200

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

$conf
           [,1]      [,2]      [,3]       [,4]       [,5]       [,6]      [,7]
[1,] -0.1900697 0.3728125 0.1876042 0.08281245 -0.1891668 -0.5676738 0.1157465
[2,]  1.1100697 0.6271875 0.2723958 0.33718755  0.4891668  0.6476738 0.3842535
          [,8]       [,9]    [,10]   [,11]    [,12]
[1,] 0.1437672 0.04454855 -0.50685 -0.1683 -0.18565
[2,] 0.8362328 0.35545145  0.30685  1.0483  0.18565

$out
[1]  1.23 -0.46  1.25  1.63 -0.42

$group
[1] 3 3 4 9 9

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

Warning message:
In bxp(list(stats = c(-0.290000000000006, -0.25, 0.459999999999994,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6655b1289573154.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]
[1,] 102.640 105.850 108.26 112.790 117.310
[2,] 103.445 106.390 109.52 114.895 117.770
[3,] 104.775 106.980 110.14 116.790 117.965
[4,] 105.715 107.405 112.00 117.025 118.370
[5,] 106.150 108.250 113.04 117.480 118.900

$n
[1] 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]
[1,] 103.7396 106.5171 109.0089 115.8185 117.6652
[2,] 105.8104 107.4429 111.2711 117.7615 118.2648

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(102.64, 103.445, 104.775, 105.715, 106.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7655b1289573154.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,] 109.4780 108.2600 109.4650
[2,] 110.4780 109.5200 110.3619
[3,] 111.0727 110.1075 111.3650
[4,] 111.4870 110.5425 111.9850
[5,] 112.4000 111.2300 112.2850

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 110.6125 109.6411 110.6247
[2,] 111.5330 110.5739 112.1053

$out
[1] 112.86

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(109.478, 110.478, 111.07275, 111.487, 112.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/12mo61289573154.ps tmp/12mo61289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/22mo61289573154.ps tmp/22mo61289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/32mo61289573154.ps tmp/32mo61289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/4ve6r1289573154.ps tmp/4ve6r1289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ve6r1289573154.ps tmp/5ve6r1289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/6655b1289573154.ps tmp/6655b1289573154.png",intern=TRUE))
character(0)
> try(system("convert tmp/7655b1289573154.ps tmp/7655b1289573154.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.040   0.620   2.666