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(0.163174,0.201921,0.232752,0.228971,0.245583,0.254378,0.224034,0.230488,0.145256,0.151829,0.165472,0.196739,0.187848,0.204079,0.207498,0.240080,0.239256,0.228852,0.245555,0.222047,0.266542,0.267053,0.231662,0.206036,0.257087,0.216368,0.225624,0.231318,0.253926,0.227463,0.263592,0.260886,0.232007,0.233059,0.299560,0.270942,0.283427,0.301381,0.316548,0.292079,0.312963,0.361612,0.358814,0.358174,0.220870,0.288982,0.291079,0.327455,0.335642,0.284180,0.325792,0.326291,0.356246,0.335943,0.414564,0.374787,0.361182,0.393189,0.344907,0.352865)
> 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,] 0.163174 0.187848 0.257087 0.283427 0.335642   NA
 [2,] 0.201921 0.204079 0.216368 0.301381 0.284180   NA
 [3,] 0.232752 0.207498 0.225624 0.316548 0.325792   NA
 [4,] 0.228971 0.240080 0.231318 0.292079 0.326291   NA
 [5,] 0.245583 0.239256 0.253926 0.312963 0.356246   NA
 [6,] 0.254378 0.228852 0.227463 0.361612 0.335943   NA
 [7,] 0.224034 0.245555 0.263592 0.358814 0.414564   NA
 [8,] 0.230488 0.222047 0.260886 0.358174 0.374787   NA
 [9,] 0.145256 0.266542 0.232007 0.220870 0.361182   NA
[10,] 0.151829 0.267053 0.233059 0.288982 0.393189   NA
[11,] 0.165472 0.231662 0.299560 0.291079 0.344907   NA
[12,] 0.196739 0.206036 0.270942 0.327455 0.352865   NA
> darr
           [,1]      [,2]      [,3]      [,4]      [,5] [,6]
 [1,]  0.038747  0.016231 -0.040719  0.017954 -0.051462   NA
 [2,]  0.030831  0.003419  0.009256  0.015167  0.041612   NA
 [3,] -0.003781  0.032582  0.005694 -0.024469  0.000499   NA
 [4,]  0.016612 -0.000824  0.022608  0.020884  0.029955   NA
 [5,]  0.008795 -0.010404 -0.026463  0.048649 -0.020303   NA
 [6,] -0.030344  0.016703  0.036129 -0.002798  0.078621   NA
 [7,]  0.006454 -0.023508 -0.002706 -0.000640 -0.039777   NA
 [8,] -0.085232  0.044495 -0.028879 -0.137304 -0.013605   NA
 [9,]  0.006573  0.000511  0.001052  0.068112  0.032007   NA
[10,]  0.013643 -0.035391  0.066501  0.002097 -0.048282   NA
[11,]  0.031267 -0.025626 -0.028618  0.036376  0.007958   NA
[12,] -0.008891  0.051051  0.012485  0.008187        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/rcomp/tmp/1s0og1289583416.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/2krn11289583416.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/3krn11289583416.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/4v1441289583416.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]
[1,] 0.163174 0.201921 0.207498 0.228971 0.239256 0.227463 0.224034 0.222047
[2,] 0.187848 0.204079 0.225624 0.231318 0.245583 0.228852 0.245555 0.230488
[3,] 0.257087 0.216368 0.232752 0.240080 0.253926 0.254378 0.263592 0.260886
[4,] 0.283427 0.284180 0.316548 0.292079 0.312963 0.335943 0.358814 0.358174
[5,] 0.335642 0.301381 0.325792 0.326291 0.356246 0.361612 0.414564 0.374787
         [,9]    [,10]    [,11]    [,12]
[1,] 0.220870 0.151829 0.165472 0.196739
[2,] 0.220870 0.233059 0.231662 0.206036
[3,] 0.232007 0.267053 0.291079 0.270942
[4,] 0.266542 0.288982 0.299560 0.327455
[5,] 0.266542 0.288982 0.344907 0.352865

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.1895511 0.1597688 0.1685053 0.1971464 0.2063155 0.1787078 0.1835635
[2,] 0.3246229 0.2729672 0.2969987 0.2830136 0.3015365 0.3300482 0.3436205
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.1706634 0.1997353 0.2275379 0.2431024 0.1851476
[2,] 0.3511086 0.2642787 0.3065681 0.3390556 0.3567364

$out
[1] 0.145256 0.361182 0.393189

$group
[1]  9  9 10

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

Warning message:
In bxp(list(stats = c(0.163174, 0.187848, 0.257087, 0.283427, 0.335642,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5v1441289583416.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]
[1,] -0.051462 0.003419 -0.003781 0.016612 -0.026463 -0.030344 -0.039777
[2,] -0.040719 0.009256 -0.003781 0.016612 -0.020303 -0.002798 -0.023508
[3,]  0.016231 0.015167  0.000499 0.020884 -0.010404  0.016703 -0.002706
[4,]  0.017954 0.030831  0.005694 0.022608  0.008795  0.036129 -0.000640
[5,]  0.038747 0.041612  0.005694 0.029955  0.048649  0.078621  0.006454
          [,8]     [,9]     [,10]     [,11]     [,12]
[1,] -0.137304 0.000511 -0.048282 -0.028618 -0.008891
[2,] -0.085232 0.001052 -0.035391 -0.025626 -0.000352
[3,] -0.028879 0.006573  0.002097  0.007958  0.010336
[4,] -0.013605 0.032007  0.013643  0.031267  0.031768
[5,]  0.044495 0.068112  0.066501  0.036376  0.051051

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

$conf
            [,1]          [,2]         [,3]       [,4]        [,5]        [,6]
[1,] -0.02522719 -7.784065e-05 -0.006196011 0.01664724 -0.03096457 -0.01080272
[2,]  0.05768919  3.041184e-02  0.007194011 0.02512076  0.01015657  0.04420872
            [,7]        [,8]        [,9]      [,10]       [,11]      [,12]
[1,] -0.01886447 -0.07949046 -0.01529973 -0.0325503 -0.03224245 -0.0150388
[2,]  0.01345247  0.02173246  0.02844573  0.0367443  0.04815845  0.0357108

$out
[1]  0.032582 -0.024469 -0.000824

$group
[1] 3 3 4

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

Warning message:
In bxp(list(stats = c(-0.051462, -0.040719, 0.0162310000000000,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6v1441289583416.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,] 0.1452560 0.1878480 0.2163680 0.2208700 0.2841800   NA
[2,] 0.1643230 0.2067670 0.2293905 0.2900305 0.3309665   NA
[3,] 0.2129775 0.2302570 0.2434925 0.3071720 0.3488860   NA
[4,] 0.2316200 0.2428175 0.2622390 0.3428145 0.3679845   NA
[5,] 0.2543780 0.2670530 0.2995600 0.3616120 0.4145640   NA

$n
[1] 12 12 12 12 12  0

$conf
          [,1]      [,2]      [,3]      [,4]      [,5] [,6]
[1,] 0.1822829 0.2138141 0.2285101 0.2830969 0.3320018   NA
[2,] 0.2436721 0.2466999 0.2584749 0.3312471 0.3657702   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.145256, 0.164323, 0.2129775, 0.23162, 0.254378,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/76smp1289583416.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,] 0.2415858 0.2163680 0.2356375
[2,] 0.2535392 0.2364160 0.2525750
[3,] 0.2666792 0.2557325 0.2661782
[4,] 0.2816222 0.2653225 0.2808352
[5,] 0.3013118 0.2910790 0.3021845

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.2538704 0.2425481 0.2532886
[2,] 0.2794880 0.2689169 0.2790679

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.2415858, 0.2535392, 0.2666792, 0.2816222,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1s0og1289583416.ps tmp/1s0og1289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/2krn11289583416.ps tmp/2krn11289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/3krn11289583416.ps tmp/3krn11289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/4v1441289583416.ps tmp/4v1441289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/5v1441289583416.ps tmp/5v1441289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/6v1441289583416.ps tmp/6v1441289583416.png",intern=TRUE))
character(0)
> try(system("convert tmp/76smp1289583416.ps tmp/76smp1289583416.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.980   0.600   2.555