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.762253,0.768403,0.757518,0.772917,0.787774,0.82203,0.830772,0.813537,0.815927,0.832293,0.848464,0.843455,0.826241,0.837661,0.831947,0.81493,0.783085,0.790514,0.788395,0.780579,0.785731,0.792959,0.776337,0.75683,0.76929,0.764877,0.755173,0.739864,0.740138,0.745212,0.729076,0.734107,0.719632,0.702889,0.681013,0.686342,0.67944,0.678058,0.644039,0.63488,0.642797,0.642963,0.634115,0.66778,0.695894,0.750638,0.785423,0.74355,0.755344,0.782167,0.766284,0.75815,0.732601,0.71347,0.709824,0.700869,0.686719,0.674946,0.670511,0.684275)
> 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.762253 0.826241 0.769290 0.679440 0.755344   NA
 [2,] 0.768403 0.837661 0.764877 0.678058 0.782167   NA
 [3,] 0.757518 0.831947 0.755173 0.644039 0.766284   NA
 [4,] 0.772917 0.814930 0.739864 0.634880 0.758150   NA
 [5,] 0.787774 0.783085 0.740138 0.642797 0.732601   NA
 [6,] 0.822030 0.790514 0.745212 0.642963 0.713470   NA
 [7,] 0.830772 0.788395 0.729076 0.634115 0.709824   NA
 [8,] 0.813537 0.780579 0.734107 0.667780 0.700869   NA
 [9,] 0.815927 0.785731 0.719632 0.695894 0.686719   NA
[10,] 0.832293 0.792959 0.702889 0.750638 0.674946   NA
[11,] 0.848464 0.776337 0.681013 0.785423 0.670511   NA
[12,] 0.843455 0.756830 0.686342 0.743550 0.684275   NA
> darr
           [,1]      [,2]      [,3]      [,4]      [,5] [,6]
 [1,]  0.006150  0.011420 -0.004413 -0.001382  0.026823   NA
 [2,] -0.010885 -0.005714 -0.009704 -0.034019 -0.015883   NA
 [3,]  0.015399 -0.017017 -0.015309 -0.009159 -0.008134   NA
 [4,]  0.014857 -0.031845  0.000274  0.007917 -0.025549   NA
 [5,]  0.034256  0.007429  0.005074  0.000166 -0.019131   NA
 [6,]  0.008742 -0.002119 -0.016136 -0.008848 -0.003646   NA
 [7,] -0.017235 -0.007816  0.005031  0.033665 -0.008955   NA
 [8,]  0.002390  0.005152 -0.014475  0.028114 -0.014150   NA
 [9,]  0.016366  0.007228 -0.016743  0.054744 -0.011773   NA
[10,]  0.016171 -0.016622 -0.021876  0.034785 -0.004435   NA
[11,] -0.005009 -0.019507  0.005329 -0.041873  0.013764   NA
[12,] -0.017214  0.012460 -0.006902  0.011794        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/1s1b81292794144.ps",horizontal=F,onefile=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/23sbt1292794144.ps",horizontal=F,onefile=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/33sbt1292794144.ps",horizontal=F,onefile=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/43sbt1292794144.ps",horizontal=F,onefile=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.755344 0.764877 0.755173 0.739864 0.732601 0.642963 0.634115 0.667780
[2,] 0.755344 0.764877 0.755173 0.739864 0.732601 0.713470 0.709824 0.700869
[3,] 0.762253 0.768403 0.757518 0.758150 0.740138 0.745212 0.729076 0.734107
[4,] 0.769290 0.782167 0.766284 0.772917 0.783085 0.790514 0.788395 0.780579
[5,] 0.769290 0.782167 0.766284 0.814930 0.787774 0.822030 0.830772 0.813537
         [,9]    [,10]    [,11]    [,12]
[1,] 0.686719 0.674946 0.670511 0.684275
[2,] 0.695894 0.702889 0.681013 0.686342
[3,] 0.719632 0.750638 0.776337 0.743550
[4,] 0.785731 0.792959 0.785423 0.756830
[5,] 0.815927 0.832293 0.848464 0.843455

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

$conf
          [,1]     [,2]     [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.7523988 0.756186 0.749667 0.7347948 0.7044661 0.6907729 0.6735579
[2,] 0.7721072 0.780620 0.765369 0.7815052 0.7758099 0.7996511 0.7845941
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.6777841 0.6561534 0.6869948 0.7025612 0.6937434
[2,] 0.7904299 0.7831106 0.8142812 0.8501128 0.7933566

$out
[1] 0.826241 0.679440 0.837661 0.678058 0.831947 0.644039 0.634880 0.642797

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

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

Warning message:
In bxp(list(stats = c(0.755344, 0.755344, 0.762253, 0.76929, 0.76929,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5ekse1292794144.ps",horizontal=F,onefile=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.004413 -0.015883 -0.017017 -0.031845 0.000166 -0.016136 -0.017235
[2,] -0.001382 -0.015883 -0.015309 -0.025549 0.000166 -0.008848 -0.008955
[3,]  0.006150 -0.010885 -0.009159  0.000274 0.005074 -0.003646 -0.007816
[4,]  0.011420 -0.009704 -0.008134  0.007917 0.007429 -0.002119  0.005031
[5,]  0.026823 -0.005714 -0.008134  0.014857 0.007429 -0.002119  0.005031
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -0.014475 -0.016743 -0.021876 -0.041873 -0.017214
[2,] -0.014150 -0.011773 -0.016622 -0.019507 -0.012058
[3,]  0.002390  0.007228 -0.004435 -0.005009  0.002446
[4,]  0.005152  0.016366  0.016171  0.005329  0.012127
[5,]  0.028114  0.054744  0.034785  0.013764  0.012460

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

$conf
             [,1]         [,2]         [,3]        [,4]          [,5]
[1,] -0.002895861 -0.015251066 -0.014228837 -0.02337299 -0.0000580175
[2,]  0.015195861 -0.006518934 -0.004089163  0.02392099  0.0102060175
             [,6]         [,7]        [,8]        [,9]       [,10]       [,11]
[1,] -0.008400694 -0.017698472 -0.01124874 -0.01265495 -0.02760645 -0.02255806
[2,]  0.001108694  0.002066472  0.01602874  0.02711095  0.01873645  0.01254006
           [,12]
[1,] -0.01666015
[2,]  0.02155215

$out
[1] -0.034019  0.015399  0.034256 -0.019131  0.008742  0.033665

$group
[1] 2 3 5 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.004413, -0.00138199999999999, 0.00614999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ekse1292794144.ps",horizontal=F,onefile=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.7575180 0.7568300 0.6810130 0.634115 0.670511   NA
[2,] 0.7706600 0.7818320 0.7112605 0.642880 0.685497   NA
[3,] 0.8147320 0.7894545 0.7369855 0.672919 0.711647   NA
[4,] 0.8315325 0.8205855 0.7501925 0.719722 0.756747   NA
[5,] 0.8484640 0.8376610 0.7692900 0.785423 0.782167   NA

$n
[1] 12 12 12 12 12  0

$conf
          [,1]      [,2]      [,3]      [,4]      [,5] [,6]
[1,] 0.7869676 0.7717788 0.7192284 0.6378708 0.6791494   NA
[2,] 0.8424964 0.8071302 0.7547426 0.7079672 0.7441446   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.757518, 0.77066, 0.814732, 0.8315325, 0.848464,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/77b9z1292794144.ps",horizontal=F,onefile=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.7372790 0.7196320 0.7215860
[2,] 0.7400775 0.7371225 0.7407682
[3,] 0.7435193 0.7479250 0.7505507
[4,] 0.7516709 0.7602015 0.7592857
[5,] 0.7662332 0.7763370 0.7735220

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.7382315 0.7373985 0.7421048
[2,] 0.7488071 0.7584515 0.7589967

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.737279, 0.7400775, 0.7435193, 0.7516709,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1s1b81292794144.ps tmp/1s1b81292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/23sbt1292794144.ps tmp/23sbt1292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/33sbt1292794144.ps tmp/33sbt1292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/43sbt1292794144.ps tmp/43sbt1292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ekse1292794144.ps tmp/5ekse1292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ekse1292794144.ps tmp/6ekse1292794144.png",intern=TRUE))
character(0)
> try(system("convert tmp/77b9z1292794144.ps tmp/77b9z1292794144.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.189   0.942   3.159