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(0.7461,0.7775,0.7790,0.7744,0.7905,0.7719,0.7811,0.7557,0.7637,0.7595,0.7471,0.7615,0.7487,0.7389,0.7337,0.7510,0.7382,0.7159,0.7542,0.7636,0.7433,0.7658,0.7627,0.7480,0.7692,0.7850,0.7913,0.7720,0.7880,0.8070,0.8268,0.8244,0.8487,0.8572,0.8214,0.8827,0.9216,0.8865,0.8816,0.8884,0.9466,0.9180,0.9337,0.9559,0.9626,0.9434,0.8639,0.7996,0.6680,0.6572,0.6928,0.6438,0.6454,0.6873,0.7265,0.7912,0.8114,0.8281,0.8393)
> 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] 59
> (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 5 4
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]
 [1,] 0.7461 0.7487 0.7692 0.9216 0.6680
 [2,] 0.7775 0.7389 0.7850 0.8865 0.6572
 [3,] 0.7790 0.7337 0.7913 0.8816 0.6928
 [4,] 0.7744 0.7510 0.7720 0.8884 0.6438
 [5,] 0.7905 0.7382 0.7880 0.9466 0.6454
 [6,] 0.7719 0.7159 0.8070 0.9180 0.6873
 [7,] 0.7811 0.7542 0.8268 0.9337 0.7265
 [8,] 0.7557 0.7636 0.8244 0.9559 0.7912
 [9,] 0.7637 0.7433 0.8487 0.9626 0.8114
[10,] 0.7595 0.7658 0.8572 0.9434 0.8281
[11,] 0.7471 0.7627 0.8214 0.8639 0.8393
[12,] 0.7615 0.7480 0.8827 0.7996     NA
> darr
         [,1]    [,2]    [,3]    [,4]    [,5]
 [1,]  0.0314 -0.0098  0.0158 -0.0351 -0.0108
 [2,]  0.0015 -0.0052  0.0063 -0.0049  0.0356
 [3,] -0.0046  0.0173 -0.0193  0.0068 -0.0490
 [4,]  0.0161 -0.0128  0.0160  0.0582  0.0016
 [5,] -0.0186 -0.0223  0.0190 -0.0286  0.0419
 [6,]  0.0092  0.0383  0.0198  0.0157  0.0392
 [7,] -0.0254  0.0094 -0.0024  0.0222  0.0647
 [8,]  0.0080 -0.0203  0.0243  0.0067  0.0202
 [9,] -0.0042  0.0225  0.0085 -0.0192  0.0167
[10,] -0.0124 -0.0031 -0.0358 -0.0795  0.0112
[11,]  0.0144 -0.0147  0.0613 -0.0643      NA
[12,] -0.0128  0.0212  0.0389 -0.1316      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/1rrvi1257450589.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/21uce1257450589.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/371zt1257450589.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/4zjtc1257450589.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,] 0.7461 0.7389 0.6928 0.7510 0.7382 0.6873 0.7265 0.7557 0.7433 0.7595
[2,] 0.7461 0.7389 0.7337 0.7510 0.7382 0.7159 0.7542 0.7636 0.7637 0.7658
[3,] 0.7487 0.7775 0.7790 0.7720 0.7880 0.7719 0.7811 0.7912 0.8114 0.8281
[4,] 0.7692 0.7850 0.7913 0.7744 0.7905 0.8070 0.8268 0.8244 0.8487 0.8572
[5,] 0.7692 0.7850 0.7913 0.7744 0.7905 0.9180 0.9337 0.8244 0.9626 0.9434
      [,11]   [,12]
[1,] 0.7471 0.74800
[2,] 0.7627 0.75475
[3,] 0.8214 0.78055
[4,] 0.8393 0.84115
[5,] 0.8639 0.88270

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

$conf
          [,1]      [,2]   [,3]      [,4]     [,5]     [,6]     [,7]      [,8]
[1,] 0.7323776 0.7449259 0.7383 0.7554656 0.751045 0.707529 0.729801 0.7482389
[2,] 0.7650224 0.8100741 0.8197 0.7885344 0.824955 0.836271 0.832399 0.8341611
          [,9]    [,10]     [,11]    [,12]
[1,] 0.7513392 0.763517 0.7672746 0.712294
[2,] 0.8714608 0.892683 0.8755254 0.848806

$out
 [1] 0.9216 0.6680 0.8865 0.6572 0.8816 0.8884 0.6438 0.9466 0.6454 0.9559

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

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

Warning message:
In bxp(list(stats = c(0.7461, 0.7461, 0.7487, 0.7692, 0.7692, 0.7389,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5tfiq1257450589.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]
[1,] -0.0351 -0.0052 -0.0490 -0.0128 -0.0286 0.0092 -0.0254 0.0067 -0.0192
[2,] -0.0108 -0.0049 -0.0193  0.0016 -0.0223 0.0157 -0.0024 0.0067 -0.0042
[3,] -0.0098  0.0015 -0.0046  0.0160 -0.0186 0.0198  0.0094 0.0080  0.0085
[4,]  0.0158  0.0063  0.0068  0.0161  0.0190 0.0383  0.0222 0.0202  0.0167
[5,]  0.0314  0.0063  0.0173  0.0161  0.0419 0.0392  0.0222 0.0243  0.0225
       [,10]    [,11]    [,12]
[1,] -0.0795 -0.06430 -0.13160
[2,] -0.0358 -0.03950 -0.07220
[3,] -0.0124 -0.00015  0.00420
[4,] -0.0031  0.03785  0.03005
[5,]  0.0112  0.06130  0.03890

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

$conf
             [,1]         [,2]        [,3]        [,4]        [,5]        [,6]
[1,] -0.028595493 -0.006413892 -0.02304219 0.005754337 -0.04778248 0.003830897
[2,]  0.008995493  0.009413892  0.01384219 0.026245663  0.01058248 0.035769103
             [,7]         [,8]         [,9]       [,10]      [,11]      [,12]
[1,] -0.007982298 -0.001539066 -0.006267887 -0.03550574 -0.0612565 -0.0765775
[2,]  0.026782298  0.017539066  0.023267887  0.01070574  0.0609565  0.0849775

$out
[1]  0.0356  0.0582  0.0647 -0.0203

$group
[1] 2 4 7 8

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

Warning message:
In bxp(list(stats = c(-0.0351, -0.0108000000000000, -0.00980000000000003,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6znyk1257450589.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,] 0.74610 0.71590 0.76920 0.79960 0.6438
[2,] 0.75760 0.73855 0.78650 0.88405 0.6626
[3,] 0.76780 0.74835 0.81420 0.91980 0.6928
[4,] 0.77825 0.75845 0.83775 0.94500 0.8013
[5,] 0.79050 0.76580 0.88270 0.96260 0.8393

$n
[1] 12 12 12 12 11

$conf
          [,1]      [,2]      [,3]      [,4]     [,5]
[1,] 0.7583814 0.7392735 0.7908245 0.8920003 0.626725
[2,] 0.7772186 0.7574265 0.8375755 0.9475997 0.758875

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.7461, 0.7576, 0.7678, 0.77825, 0.7905, 0.7159,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7c6k21257450589.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.765920 0.748700 0.757650
[2,] 0.773200 0.774750 0.762225
[3,] 0.789845 0.780825 0.776800
[4,] 0.812520 0.801300 0.797500
[5,] 0.830800 0.828100 0.811500

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.7719109 0.7687154 0.7607108
[2,] 0.8077791 0.7929346 0.7928892

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.76592, 0.7732, 0.789845, 0.81252, 0.8308,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1rrvi1257450589.ps tmp/1rrvi1257450589.png")
> system("convert tmp/21uce1257450589.ps tmp/21uce1257450589.png")
> system("convert tmp/371zt1257450589.ps tmp/371zt1257450589.png")
> system("convert tmp/4zjtc1257450589.ps tmp/4zjtc1257450589.png")
> system("convert tmp/5tfiq1257450589.ps tmp/5tfiq1257450589.png")
> system("convert tmp/6znyk1257450589.ps tmp/6znyk1257450589.png")
> system("convert tmp/7c6k21257450589.ps tmp/7c6k21257450589.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.205   0.959   1.823