R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.60773,0.58933,0.60039,0.61342,0.6348,0.634,0.62915,0.62168,0.61328,0.6089,0.60857,0.62672,0.62291,0.62393,0.61838,0.62012,0.61659,0.6116,0.61573,0.61407,0.62823,0.64405,0.6387,0.63633,0.63059,0.62994,0.63709,0.64217,0.65711,0.66977,0.68255,0.68902,0.71322,0.70224,0.70045,0.69919,0.69693,0.69763,0.69278,0.70196,0.69215,0.6769,0.67124,0.66532,0.67157,0.66428,0.66576,0.66942,0.6813,0.69144,0.69862,0.695,0.69867,0.68968,0.69233,0.68293,0.68399,0.66895,0.68756,0.68527,0.6776,0.68137,0.67933,0.67922,0.68598,0.68297,0.68935,0.69463,0.6833,0.68666,0.68782,0.67669,0.67511,0.67254,0.67397,0.67286,0.66341,0.668,0.68021,0.67934,0.68136,0.67562,0.6744,0.67766,0.68887,0.69614,0.70896,0.72064)
> 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] 88
> (np <- floor(n / par1))
[1] 7
> 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] 8 8 8 8 7 7 7 7 7 7 7 7
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
 [1,] 0.60773 0.62291 0.63059 0.69693 0.68130 0.67760 0.67511 0.68887
 [2,] 0.58933 0.62393 0.62994 0.69763 0.69144 0.68137 0.67254 0.69614
 [3,] 0.60039 0.61838 0.63709 0.69278 0.69862 0.67933 0.67397 0.70896
 [4,] 0.61342 0.62012 0.64217 0.70196 0.69500 0.67922 0.67286 0.72064
 [5,] 0.63480 0.61659 0.65711 0.69215 0.69867 0.68598 0.66341      NA
 [6,] 0.63400 0.61160 0.66977 0.67690 0.68968 0.68297 0.66800      NA
 [7,] 0.62915 0.61573 0.68255 0.67124 0.69233 0.68935 0.68021      NA
 [8,] 0.62168 0.61407 0.68902 0.66532 0.68293 0.69463 0.67934      NA
 [9,] 0.61328 0.62823 0.71322 0.67157 0.68399 0.68330 0.68136      NA
[10,] 0.60890 0.64405 0.70224 0.66428 0.66895 0.68666 0.67562      NA
[11,] 0.60857 0.63870 0.70045 0.66576 0.68756 0.68782 0.67440      NA
[12,] 0.62672 0.63633 0.69919 0.66942 0.68527 0.67669 0.67766      NA
> darr
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
 [1,] -0.01840  0.00102 -0.00065  0.00070  0.01014  0.00377 -0.00257 0.00727
 [2,]  0.01106 -0.00555  0.00715 -0.00485  0.00718 -0.00204  0.00143 0.01282
 [3,]  0.01303  0.00174  0.00508  0.00918 -0.00362 -0.00011 -0.00111 0.01168
 [4,]  0.02138 -0.00353  0.01494 -0.00981  0.00367  0.00676 -0.00945      NA
 [5,] -0.00080 -0.00499  0.01266 -0.01525 -0.00899 -0.00301  0.00459      NA
 [6,] -0.00485  0.00413  0.01278 -0.00566  0.00265  0.00638  0.01221      NA
 [7,] -0.00747 -0.00166  0.00647 -0.00592 -0.00940  0.00528 -0.00087      NA
 [8,] -0.00840  0.01416  0.02420  0.00625  0.00106 -0.01133  0.00202      NA
 [9,] -0.00438  0.01582 -0.01098 -0.00729 -0.01504  0.00336 -0.00574      NA
[10,] -0.00033 -0.00535 -0.00179  0.00148  0.01861  0.00116 -0.00122      NA
[11,]  0.01815 -0.00237 -0.00126  0.00366 -0.00229 -0.01113  0.00326      NA
[12,] -0.00381 -0.00574 -0.00226  0.01188 -0.00767 -0.00158  0.01121      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/1kemv1227967420.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/2uz061227967420.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/3f9i81227967420.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/4iwhb1227967420.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.607730 0.589330 0.600390 0.613420 0.616590 0.611600 0.615730 0.614070
[2,] 0.626750 0.626935 0.627735 0.631145 0.645955 0.651000 0.650195 0.643500
[3,] 0.676355 0.676955 0.676650 0.676040 0.663410 0.669770 0.680210 0.679340
[4,] 0.685085 0.693790 0.695700 0.698480 0.689065 0.679935 0.685950 0.685975
[5,] 0.696930 0.697630 0.708960 0.720640 0.698670 0.689680 0.692330 0.694630
         [,9]    [,10]   [,11]    [,12]
[1,] 0.613280 0.644050 0.60857 0.626720
[2,] 0.649900 0.654165 0.65223 0.652875
[3,] 0.681360 0.668950 0.67440 0.676690
[4,] 0.683645 0.681140 0.68769 0.681465
[5,] 0.713220 0.702240 0.70045 0.699190

$n
 [1] 8 8 8 8 7 7 7 7 7 7 7 7

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.6437682 0.6396088 0.6386838 0.6384257 0.6376654 0.6524905 0.6588577
[2,] 0.7089418 0.7143012 0.7146162 0.7136543 0.6891546 0.6870495 0.7015623
          [,8]     [,9]    [,10]     [,11]     [,12]
[1,] 0.6539746 0.661208 0.652841 0.6532239 0.6596165
[2,] 0.7047054 0.701512 0.685059 0.6955761 0.6937635

$out
[1] 0.6089

$group
[1] 10

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

Warning message:
In bxp(list(stats = c(0.60773, 0.62675, 0.676355, 0.685085, 0.69693,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/51xqg1227967421.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.00257 -0.005550 -0.00362 -0.00981 -0.015250 -0.005660 -0.009400
[2,] -0.00161 -0.003445 -0.00061 -0.00649 -0.006990 -0.001100 -0.006695
[3,]  0.00086  0.004290  0.00341  0.00367 -0.003010  0.004130 -0.001660
[4,]  0.00552  0.009120  0.01043  0.01085  0.001895  0.009295  0.002205
[5,]  0.01014  0.012820  0.01303  0.02138  0.012660  0.012780  0.006470
          [,8]      [,9]     [,10]    [,11]     [,12]
[1,] -0.011330 -0.015040 -0.005350 -0.00237 -0.007670
[2,] -0.003670 -0.009135 -0.001505 -0.00233 -0.004775
[3,]  0.002020 -0.005740 -0.000330 -0.00126 -0.002260
[4,]  0.010205 -0.000510  0.001320  0.00346  0.004815
[5,]  0.024200  0.003360  0.001480  0.00366  0.011880

$n
 [1] 8 8 8 7 7 7 7 7 7 7 7 7

$conf
             [,1]         [,2]         [,3]         [,4]         [,5]
[1,] -0.003122920 -0.002728989 -0.002757103 -0.006685168 -0.008315979
[2,]  0.004842920  0.011308989  0.009577103  0.014025168  0.002295979
             [,6]         [,7]         [,8]          [,9]        [,10]
[1,] -0.002077726 -0.006974936 -0.006265926 -0.0108907109 -0.002017044
[2,]  0.010337726  0.003654936  0.010305926 -0.0005892891  0.001357044
            [,11]        [,12]
[1,] -0.004717695 -0.007986993
[2,]  0.002197695  0.003466993

$out
[1] -0.01840  0.01582  0.01861  0.01815 -0.01113

$group
[1]  1  9 10 11 11

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

Warning message:
In bxp(list(stats = c(-0.00256999999999996, -0.00160999999999994,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6dfmw1227967421.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,] 0.589330 0.611600 0.62994 0.664280 0.668950 0.676690 0.668000 0.688870
[2,] 0.608150 0.616160 0.63963 0.667590 0.683460 0.679275 0.672700 0.692505
[3,] 0.613350 0.621515 0.67616 0.674235 0.688620 0.683135 0.674755 0.702550
[4,] 0.627935 0.632280 0.69982 0.694855 0.693665 0.687240 0.678500 0.714800
[5,] 0.634800 0.644050 0.71322 0.701960 0.698670 0.694630 0.681360 0.720640

$n
[1] 12 12 12 12 12 12 12  4

$conf
          [,1]      [,2]     [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.6043259 0.6141626 0.648707 0.6617993 0.6839654 0.6795021 0.6721096
[2,] 0.6223741 0.6288674 0.703613 0.6866707 0.6932746 0.6867679 0.6774004
         [,8]
[1,] 0.684937
[2,] 0.720163

$out
[1] 0.66341

$group
[1] 7

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

Warning message:
In bxp(list(stats = c(0.58933, 0.60815, 0.61335, 0.627935, 0.6348,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7lj891227967421.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.6601300 0.6634100 0.6605262
[2,] 0.6627679 0.6720850 0.6640319
[3,] 0.6642436 0.6765025 0.6667356
[4,] 0.6667529 0.6781475 0.6675813
[5,] 0.6681738 0.6813600 0.6699600

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.6624260 0.6737374 0.6651167
[2,] 0.6660612 0.6792676 0.6683545

$out
[1] 0.6559312

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(0.66013, 0.662767857142857, 0.664243571428571,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1kemv1227967420.ps tmp/1kemv1227967420.png")
> system("convert tmp/2uz061227967420.ps tmp/2uz061227967420.png")
> system("convert tmp/3f9i81227967420.ps tmp/3f9i81227967420.png")
> system("convert tmp/4iwhb1227967420.ps tmp/4iwhb1227967420.png")
> system("convert tmp/51xqg1227967421.ps tmp/51xqg1227967421.png")
> system("convert tmp/6dfmw1227967421.ps tmp/6dfmw1227967421.png")
> system("convert tmp/7lj891227967421.ps tmp/7lj891227967421.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.184   0.928   1.741