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.892201,1.007095,0.966709,1.003381,0.950054,0.942173,0.938196,0.731490,0.815207,0.841577,0.896191,0.918043,0.819296,0.907075,0.892324,0.927952,0.897811,0.917274,0.909615,0.876498,0.918068,0.888221,0.928924,0.921158,0.832321,0.918037,0.902913,0.934450,0.886786,0.903352,0.883465,0.849764,0.906730,0.872937,0.909101,0.916899,0.861511,0.926093,0.890331,0.917292,0.887005,0.922855,0.879873,0.683812,0.839532,0.874177,0.909158,0.926634,0.854157,0.939830,0.917559,0.967147,0.918967,0.971321,0.974110,0.965109,0.989029,0.944283,0.980111,0.976797)
> 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.892201 0.819296 0.832321 0.861511 0.854157   NA
 [2,] 1.007095 0.907075 0.918037 0.926093 0.939830   NA
 [3,] 0.966709 0.892324 0.902913 0.890331 0.917559   NA
 [4,] 1.003381 0.927952 0.934450 0.917292 0.967147   NA
 [5,] 0.950054 0.897811 0.886786 0.887005 0.918967   NA
 [6,] 0.942173 0.917274 0.903352 0.922855 0.971321   NA
 [7,] 0.938196 0.909615 0.883465 0.879873 0.974110   NA
 [8,] 0.731490 0.876498 0.849764 0.683812 0.965109   NA
 [9,] 0.815207 0.918068 0.906730 0.839532 0.989029   NA
[10,] 0.841577 0.888221 0.872937 0.874177 0.944283   NA
[11,] 0.896191 0.928924 0.909101 0.909158 0.980111   NA
[12,] 0.918043 0.921158 0.916899 0.926634 0.976797   NA
> darr
           [,1]      [,2]      [,3]      [,4]      [,5] [,6]
 [1,]  0.114894  0.087779  0.085716  0.064582  0.085673   NA
 [2,] -0.040386 -0.014751 -0.015124 -0.035762 -0.022271   NA
 [3,]  0.036672  0.035628  0.031537  0.026961  0.049588   NA
 [4,] -0.053327 -0.030141 -0.047664 -0.030287 -0.048180   NA
 [5,] -0.007881  0.019463  0.016566  0.035850  0.052354   NA
 [6,] -0.003977 -0.007659 -0.019887 -0.042982  0.002789   NA
 [7,] -0.206706 -0.033117 -0.033701 -0.196061 -0.009001   NA
 [8,]  0.083717  0.041570  0.056966  0.155720  0.023920   NA
 [9,]  0.026370 -0.029847 -0.033793  0.034645 -0.044746   NA
[10,]  0.054614  0.040703  0.036164  0.034981  0.035828   NA
[11,]  0.021852 -0.007766  0.007798  0.017476 -0.003314   NA
[12,] -0.098747 -0.088837 -0.055388 -0.072477        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/1kt421289581503.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/2kt421289581503.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/3kt421289581503.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/45c6i1289581504.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.819296 0.907075 0.890331 0.917292 0.886786 0.903352 0.879873 0.683812
[2,] 0.832321 0.918037 0.892324 0.927952 0.887005 0.917274 0.883465 0.731490
[3,] 0.854157 0.926093 0.902913 0.934450 0.897811 0.922855 0.909615 0.849764
[4,] 0.861511 0.939830 0.917559 0.967147 0.918967 0.942173 0.938196 0.876498
[5,] 0.892201 0.939830 0.917559 1.003381 0.950054 0.971321 0.974110 0.965109
         [,9]    [,10]    [,11]    [,12]
[1,] 0.815207 0.872937 0.896191 0.916899
[2,] 0.839532 0.872937 0.909101 0.918043
[3,] 0.906730 0.874177 0.909158 0.921158
[4,] 0.918068 0.888221 0.928924 0.926634
[5,] 0.989029 0.888221 0.928924 0.926634

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

$conf
          [,1]      [,2]     [,3]     [,4]      [,5]      [,6]      [,7]
[1,] 0.8335314 0.9106941 0.885082 0.906755 0.8752267 0.9052614 0.8709422
[2,] 0.8747826 0.9414919 0.920744 0.962145 0.9203953 0.9404486 0.9482878
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.7473017 0.8512367 0.8633774 0.8951511 0.9150876
[2,] 0.9522263 0.9622233 0.8849766 0.9231649 0.9272284

$out
[1] 1.007095 0.966709 0.841577 0.944283 0.980111 0.976797

$group
[1]  2  3 10 10 11 12

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

Warning message:
In bxp(list(stats = c(0.819296, 0.832321, 0.854157, 0.861511, 0.892201,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/55c6i1289581504.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.085673 -0.040386 0.026961 -0.053327 -0.007881 -0.042982 -0.206706
[2,] 0.085673 -0.035762 0.031537 -0.048180  0.016566 -0.019887 -0.196061
[3,] 0.085716 -0.022271 0.035628 -0.047664  0.019463 -0.007659 -0.033701
[4,] 0.087779 -0.015124 0.036672 -0.030287  0.035850 -0.003977 -0.033117
[5,] 0.087779 -0.014751 0.036672 -0.030141  0.052354  0.002789 -0.009001
         [,8]      [,9]    [,10]     [,11]      [,12]
[1,] 0.023920 -0.044746 0.034981 -0.007766 -0.0987470
[2,] 0.041570 -0.033793 0.035828 -0.003314 -0.0937920
[3,] 0.056966 -0.029847 0.036164  0.007798 -0.0806570
[4,] 0.083717  0.026370 0.040703  0.017476 -0.0639325
[5,] 0.083717  0.034645 0.040703  0.021852 -0.0553880

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

$conf
          [,1]         [,2]       [,3]        [,4]        [,5]         [,6]
[1,] 0.0842279 -0.036853759 0.03199962 -0.06030715 0.005836974 -0.018900966
[2,] 0.0872041 -0.007688241 0.03925638 -0.03502085 0.033089026  0.003582966
            [,7]       [,8]        [,9]      [,10]        [,11]     [,12]
[1,] -0.14883682 0.02718504 -0.07235802 0.03271934 -0.006892162 -0.104246
[2,]  0.08143482 0.08674696  0.01266402 0.03960866  0.022488162 -0.057068

$out
[1] 0.114894 0.064582 0.049588 0.155720 0.054614

$group
[1]  1  1  3  8 10

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

Warning message:
In bxp(list(stats = c(0.085673, 0.085673, 0.085716, 0.087779, 0.087779,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/65c6i1289581504.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.7314900 0.8764980 0.8323210 0.8395320 0.9175590   NA
[2,] 0.8668890 0.8902725 0.8782010 0.8678440 0.9293985   NA
[3,] 0.9281195 0.9083450 0.9031325 0.8886680 0.9661280   NA
[4,] 0.9583815 0.9196130 0.9130000 0.9200735 0.9754535   NA
[5,] 1.0070950 0.9289240 0.9344500 0.9266340 0.9890290   NA

$n
[1] 12 12 12 12 12  0

$conf
          [,1]      [,2]      [,3]      [,4]     [,5] [,6]
[1,] 0.8863892 0.8949626 0.8872604 0.8648458 0.945122   NA
[2,] 0.9698498 0.9217274 0.9190046 0.9124902 0.987134   NA

$out
[1] 0.819296 0.683812 0.854157

$group
[1] 2 4 5

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

Warning message:
In bxp(list(stats = c(0.73149, 0.866889, 0.9281195, 0.9583815, 1.007095,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7rd561289581504.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.8518972 0.8497640 0.8469160
[2,] 0.8889761 0.8859940 0.8796895
[3,] 0.9155095 0.9079440 0.9078860
[4,] 0.9316506 0.9220065 0.9256360
[5,] 0.9500444 0.9344500 0.9475495

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.8960454 0.8915185 0.8869295
[2,] 0.9349736 0.9243695 0.9288425

$out
[1] 0.8213346 0.8039940

$group
[1] 1 3

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

Warning message:
In bxp(list(stats = c(0.8518972, 0.8889761, 0.9155095, 0.9316506,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1kt421289581503.ps tmp/1kt421289581503.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kt421289581503.ps tmp/2kt421289581503.png",intern=TRUE))
character(0)
> try(system("convert tmp/3kt421289581503.ps tmp/3kt421289581503.png",intern=TRUE))
character(0)
> try(system("convert tmp/45c6i1289581504.ps tmp/45c6i1289581504.png",intern=TRUE))
character(0)
> try(system("convert tmp/55c6i1289581504.ps tmp/55c6i1289581504.png",intern=TRUE))
character(0)
> try(system("convert tmp/65c6i1289581504.ps tmp/65c6i1289581504.png",intern=TRUE))
character(0)
> try(system("convert tmp/7rd561289581504.ps tmp/7rd561289581504.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.219   0.954   1.622