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(1.01135,1.01747,1.03066,1.04183,1.05315,1.06308,1.07429,1.08626,1.10025,1.12714,1.14231,1.14965,1.16536,1.18178,1.19467,1.21240,1.22326,1.23535,1.25060,1.26165,1.27666,1.29593,1.30998,1.32252,1.33510,1.34013,1.35081,1.36653,1.37428,1.38213,1.39470,1.40414,1.41745,1.43398,1.39402,1.40632,1.41982,1.43532,1.44687,1.46064,1.47084,1.48212,1.49316,1.50550,1.51772,1.31869,1.33392,1.34300,1.35538,1.21333,1.18188,1.10824,1.04660,1.05462,1.06373,0.95098,0.95915,0.97592,0.99204,1.00234)
> par1 = '4'
> #'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] 15
> 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] 15 15 15 15
> arr
        [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
[1,] 1.01135 1.05315 1.10025 1.16536 1.22326 1.27666 1.33510 1.37428 1.41745
[2,] 1.01747 1.06308 1.12714 1.18178 1.23535 1.29593 1.34013 1.38213 1.43398
[3,] 1.03066 1.07429 1.14231 1.19467 1.25060 1.30998 1.35081 1.39470 1.39402
[4,] 1.04183 1.08626 1.14965 1.21240 1.26165 1.32252 1.36653 1.40414 1.40632
       [,10]   [,11]   [,12]   [,13]   [,14]   [,15] [,16]
[1,] 1.41982 1.47084 1.51772 1.35538 1.04660 0.95915    NA
[2,] 1.43532 1.48212 1.31869 1.21333 1.05462 0.97592    NA
[3,] 1.44687 1.49316 1.33392 1.18188 1.06373 0.99204    NA
[4,] 1.46064 1.50550 1.34300 1.10824 0.95098 1.00234    NA
> darr
        [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]     [,9]
[1,] 0.00612 0.00993 0.02689 0.01642 0.01209 0.01927 0.00503 0.00785  0.01653
[2,] 0.01319 0.01121 0.01517 0.01289 0.01525 0.01405 0.01068 0.01257 -0.03996
[3,] 0.01117 0.01197 0.00734 0.01773 0.01105 0.01254 0.01572 0.00944  0.01230
[4,] 0.01132 0.01399 0.01571 0.01086 0.01501 0.01258 0.00775 0.01331  0.01350
       [,10]   [,11]    [,12]    [,13]    [,14]   [,15] [,16]
[1,] 0.01550 0.01128 -0.19903 -0.14205  0.00802 0.01677    NA
[2,] 0.01155 0.01104  0.01523 -0.03145  0.00911 0.01612    NA
[3,] 0.01377 0.01234  0.00908 -0.07364 -0.11275 0.01030    NA
[4,] 0.01020 0.01222  0.01238 -0.06164  0.00817      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/1hsmn1258108136.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/2lart1258108136.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/36mip1258108136.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/4sh4z1258108136.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]
[1,] 0.959150 0.97592 0.992040 0.950980
[2,] 1.076700 1.09511 1.108300 1.097250
[3,] 1.276660 1.23535 1.250600 1.261650
[4,] 1.395865 1.36113 1.372415 1.385335
[5,] 1.517720 1.48212 1.493160 1.505500

$n
[1] 15 15 15 15

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 1.146455 1.126826 1.142853 1.144125
[2,] 1.406865 1.343874 1.358347 1.379175

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.95915, 1.0767, 1.27666, 1.395865, 1.51772,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5n8821258108136.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]
[1,] 0.005030 0.00911 0.00734 0.00775
[2,] 0.006985 0.01086 0.00926 0.01020
[3,] 0.011280 0.01257 0.01117 0.01230
[4,] 0.016475 0.01461 0.01244 0.01350
[5,] 0.026890 0.01612 0.01572 0.01571

$n
[1] 15 15 15 14

$conf
            [,1]       [,2]        [,3]      [,4]
[1,] 0.007408514 0.01104017 0.009872705 0.0109065
[2,] 0.015151486 0.01409983 0.012467295 0.0136935

$out
[1] -0.19903 -0.14205 -0.03996 -0.03145  0.01773 -0.07364 -0.11275 -0.06164

$group
[1] 1 1 2 2 3 3 3 4

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

Warning message:
In bxp(list(stats = c(0.00503000000000009, 0.00698500000000013,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6adkb1258108136.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]     [,8]
[1,] 1.011350 1.053150 1.100250 1.165360 1.223260 1.276660 1.335100 1.374280
[2,] 1.014410 1.058115 1.113695 1.173570 1.229305 1.286295 1.337615 1.378205
[3,] 1.024065 1.068685 1.134725 1.188225 1.242975 1.302955 1.345470 1.388415
[4,] 1.036245 1.080275 1.145980 1.203535 1.256125 1.316250 1.358670 1.399420
[5,] 1.041830 1.086260 1.149650 1.212400 1.261650 1.322520 1.366530 1.404140
         [,9]    [,10]   [,11]    [,12]    [,13]    [,14]    [,15] [,16]
[1,] 1.394020 1.419820 1.47084 1.318690 1.108240 0.950980 0.959150    NA
[2,] 1.400170 1.427570 1.47648 1.326305 1.145060 0.998790 0.967535    NA
[3,] 1.411885 1.441095 1.48764 1.338460 1.197605 1.050610 0.983980    NA
[4,] 1.425715 1.453755 1.49933 1.430360 1.284355 1.059175 0.997190    NA
[5,] 1.433980 1.460640 1.50550 1.517720 1.355380 1.063730 1.002340    NA

$n
 [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1.006815 1.051179 1.109220 1.164553 1.221787 1.279291 1.328837 1.371655
[2,] 1.041315 1.086191 1.160230 1.211897 1.264163 1.326619 1.362103 1.405175
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]     [,15] [,16]
[1,] 1.391704 1.420409 1.469588 1.256257 1.087562 1.002906 0.9605526    NA
[2,] 1.432066 1.461781 1.505692 1.420663 1.307648 1.098314 1.0074075    NA

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] NA  

Warning message:
In bxp(list(stats = c(1.01135, 1.01441, 1.024065, 1.036245, 1.04183,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7tedx1258108136.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,] 1.237133 1.235350 1.228120
[2,] 1.239300 1.242975 1.232201
[3,] 1.242521 1.256125 1.238320
[4,] 1.246000 1.269155 1.240825
[5,] 1.248425 1.276660 1.241292

$n
[1] 4 4 4

$conf
         [,1]     [,2]     [,3]
[1,] 1.237228 1.235443 1.231507
[2,] 1.247815 1.276807 1.245133

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.23713266666667, 1.23929966666667, 1.24252133333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1hsmn1258108136.ps tmp/1hsmn1258108136.png")
> system("convert tmp/2lart1258108136.ps tmp/2lart1258108136.png")
> system("convert tmp/36mip1258108136.ps tmp/36mip1258108136.png")
> system("convert tmp/4sh4z1258108136.ps tmp/4sh4z1258108136.png")
> system("convert tmp/5n8821258108136.ps tmp/5n8821258108136.png")
> system("convert tmp/6adkb1258108136.ps tmp/6adkb1258108136.png")
> system("convert tmp/7tedx1258108136.ps tmp/7tedx1258108136.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.191   0.920   1.566