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(100,100,100,100,101.1,101.35,101.45,101.49,101.68,101.92,102.04,102.55,104.02,105.41,105.48,105.54,105.16,105.16,105.16,105.16,105.16,105.17,105.17,105.54,106.9,107.27,107.31,107.39,107.41,107.46,113.14,117,119.28,119.39,119.5,119.67,119.67,119.73,119.77,119.77,119.78,119.78,119.78,121.28,122.44,122.72,122.75,122.8,122.81,122.83,122.83,122.83,122.84,122.85,123.61,124.74,125.1,125.29,125.45,125.51,125.55)
> 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] 61
> (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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]
 [1,] 100.00 104.02 106.90 119.67 122.81 125.55
 [2,] 100.00 105.41 107.27 119.73 122.83     NA
 [3,] 100.00 105.48 107.31 119.77 122.83     NA
 [4,] 100.00 105.54 107.39 119.77 122.83     NA
 [5,] 101.10 105.16 107.41 119.78 122.84     NA
 [6,] 101.35 105.16 107.46 119.78 122.85     NA
 [7,] 101.45 105.16 113.14 119.78 123.61     NA
 [8,] 101.49 105.16 117.00 121.28 124.74     NA
 [9,] 101.68 105.16 119.28 122.44 125.10     NA
[10,] 101.92 105.17 119.39 122.72 125.29     NA
[11,] 102.04 105.17 119.50 122.75 125.45     NA
[12,] 102.55 105.54 119.67 122.80 125.51     NA
> darr
      [,1]  [,2] [,3] [,4] [,5] [,6]
 [1,] 0.00  1.39 0.37 0.06 0.02   NA
 [2,] 0.00  0.07 0.04 0.04 0.00   NA
 [3,] 0.00  0.06 0.08 0.00 0.00   NA
 [4,] 1.10 -0.38 0.02 0.01 0.01   NA
 [5,] 0.25  0.00 0.05 0.00 0.01   NA
 [6,] 0.10  0.00 5.68 0.00 0.76   NA
 [7,] 0.04  0.00 3.86 1.50 1.13   NA
 [8,] 0.19  0.00 2.28 1.16 0.36   NA
 [9,] 0.24  0.01 0.11 0.28 0.19   NA
[10,] 0.12  0.00 0.11 0.03 0.16   NA
[11,] 0.51  0.37 0.17 0.05 0.06   NA
[12,] 1.47  1.36 0.00 0.01 0.04   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/17c1q1258138507.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/2brvm1258138507.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/3cbo91258138507.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/4wqoa1258138507.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,] 100.000 100.00 100.00 100.00 101.10 101.35 101.45 101.49 101.68 101.92
[2,] 104.020 105.41 105.48 105.54 105.16 105.16 105.16 105.16 105.16 105.17
[3,] 113.285 107.27 107.31 107.39 107.41 107.46 113.14 117.00 119.28 119.39
[4,] 122.810 119.73 119.77 119.77 119.78 119.78 119.78 121.28 122.44 122.72
[5,] 125.550 122.83 122.83 122.83 122.84 122.85 123.61 124.74 125.10 125.29
      [,11]  [,12]
[1,] 102.04 102.55
[2,] 105.17 105.54
[3,] 119.50 119.67
[4,] 122.75 122.80
[5,] 125.45 125.51

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

$conf
         [,1]      [,2]      [,3]      [,4]      [,5]      [,6]     [,7]
[1,] 101.1648  97.15152  97.21272  97.33512  97.07954  97.12954 102.8095
[2,] 125.4052 117.38848 117.40728 117.44488 117.74046 117.79046 123.4705
         [,8]   [,9]    [,10]    [,11]    [,12]
[1,] 105.6096 107.07 106.9892 107.0780 107.4741
[2,] 128.3904 131.49 131.7908 131.9220 131.8659

$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(100, 104.02, 113.285, 122.81, 125.55, 100,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/56qdx1258138507.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.00 0.00 0.00 0.01 0.00 0.00 0.00 0.00 0.01  0.00  0.05  0.00
[2,] 0.02 0.00 0.00 0.01 0.00 0.00 0.04 0.19 0.11  0.03  0.06  0.01
[3,] 0.06 0.04 0.00 0.01 0.01 0.10 1.13 0.36 0.19  0.11  0.17  0.04
[4,] 0.37 0.04 0.06 0.02 0.05 0.76 1.50 1.16 0.24  0.12  0.37  1.36
[5,] 0.37 0.07 0.08 0.02 0.05 0.76 1.50 2.28 0.28  0.16  0.51  1.47

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

$conf
           [,1]      [,2]        [,3]        [,4]        [,5]       [,6]
[1,] -0.1873091 0.0117361 -0.04239585 0.002934025 -0.02532987 -0.4370141
[2,]  0.3073091 0.0682639  0.04239585 0.017065975  0.04532987  0.6370141
           [,7]       [,8]       [,9]      [,10]       [,11]      [,12]
[1,] 0.09836768 -0.3253996 0.09814233 0.04640623 -0.04904522 -0.9139066
[2,] 2.16163232  1.0453996 0.28185767 0.17359377  0.38904522  0.9939066

$out
[1]  1.39  1.10 -0.38  0.25  5.68  3.86

$group
[1] 1 4 4 5 6 7

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

Warning message:
In bxp(list(stats = c(0, 0.019999999999996, 0.0600000000000023,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6vaqj1258138507.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,] 100.00 105.160 106.900 119.67 122.810 125.55
[2,] 100.00 105.160 107.350 119.77 122.830 125.55
[3,] 101.40 105.165 110.300 119.78 123.230 125.55
[4,] 101.80 105.445 119.335 122.58 125.195 125.55
[5,] 102.55 105.540 119.670 122.80 125.510 125.55

$n
[1] 12 12 12 12 12  1

$conf
        [,1]    [,2]     [,3]     [,4]     [,5]   [,6]
[1,] 100.579 105.035 104.8336 118.4983 122.1513 125.55
[2,] 102.221 105.295 115.7664 121.0617 124.3087 125.55

$out
[1] 104.02

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(100, 100, 101.4, 101.8, 102.55, 105.16, 105.16,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/78rpw1258138507.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.0480 107.2700 112.4700
[2,] 111.1820 107.4000 112.5200
[3,] 112.8932 113.2125 112.9375
[4,] 114.8150 119.3350 113.8725
[5,] 115.2140 119.6700 114.1700

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 111.2361 107.7689 112.3206
[2,] 114.5502 118.6561 113.5544

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(111.048, 111.182, 112.893166666667, 114.815,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/17c1q1258138507.ps tmp/17c1q1258138507.png")
> system("convert tmp/2brvm1258138507.ps tmp/2brvm1258138507.png")
> system("convert tmp/3cbo91258138507.ps tmp/3cbo91258138507.png")
> system("convert tmp/4wqoa1258138507.ps tmp/4wqoa1258138507.png")
> system("convert tmp/56qdx1258138507.ps tmp/56qdx1258138507.png")
> system("convert tmp/6vaqj1258138507.ps tmp/6vaqj1258138507.png")
> system("convert tmp/78rpw1258138507.ps tmp/78rpw1258138507.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.253   0.992   1.939