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(599,588,566,557,561,549,532,526,511,499,555,565,542,527,510,514,517,508,493,490,469,478,528,534,518,506,502,516,528,533,536,537,524,536,587,597,581,564,558,575,580,575,563,552,537,545,601,604,586)
> 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] 49
> (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 4 4 4 4 4 4 4 4 4 4 4
> arr
      [,1] [,2] [,3] [,4] [,5]
 [1,]  599  542  518  581  586
 [2,]  588  527  506  564   NA
 [3,]  566  510  502  558   NA
 [4,]  557  514  516  575   NA
 [5,]  561  517  528  580   NA
 [6,]  549  508  533  575   NA
 [7,]  532  493  536  563   NA
 [8,]  526  490  537  552   NA
 [9,]  511  469  524  537   NA
[10,]  499  478  536  545   NA
[11,]  555  528  587  601   NA
[12,]  565  534  597  604   NA
> darr
      [,1] [,2] [,3] [,4] [,5]
 [1,]  -11  -15  -12  -17   NA
 [2,]  -22  -17   -4   -6   NA
 [3,]   -9    4   14   17   NA
 [4,]    4    3   12    5   NA
 [5,]  -12   -9    5   -5   NA
 [6,]  -17  -15    3  -12   NA
 [7,]   -6   -3    1  -11   NA
 [8,]  -15  -21  -13  -15   NA
 [9,]  -12    9   12    8   NA
[10,]   56   50   51   56   NA
[11,]   10    6   10    3   NA
[12,]  -23  -16  -16  -18   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/126z11289754041.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/226z11289754041.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/326z11289754041.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/4dfgm1289754041.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] [,11] [,12]
[1,]  518 506.0  502 514.0 517.0 508.0 493.0 490.0 469.0 478.0 528.0 534.0
[2,]  542 516.5  506 515.0 522.5 520.5 512.5 508.0 490.0 488.5 541.5 549.5
[3,]  581 545.5  534 536.5 544.5 541.0 534.0 531.5 517.5 517.5 571.0 581.0
[4,]  586 576.0  562 566.0 570.5 562.0 549.5 544.5 530.5 540.5 594.0 600.5
[5,]  599 588.0  566 575.0 580.0 575.0 563.0 552.0 537.0 545.0 601.0 604.0

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

$conf
         [,1]    [,2]   [,3]   [,4]   [,5]    [,6]   [,7]    [,8]    [,9]
[1,] 549.9097 498.495 489.76 496.21 506.58 508.215 504.77 502.665 485.505
[2,] 612.0903 592.505 578.24 576.79 582.42 573.785 563.23 560.335 549.495
      [,10]   [,11]  [,12]
[1,] 476.42 529.525 540.71
[2,] 558.58 612.475 621.29

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(518, 542, 581, 586, 599, 506, 516.5, 545.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5dfgm1289754041.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] [,10] [,11] [,12]
[1,] -17.0 -22.0 -9.0  3.0 -12.0 -17.0 -11.0  -21 -12.0  50.0   3.0 -23.0
[2,] -16.0 -19.5 -2.5  3.5 -10.5 -16.0  -8.5  -18  -2.0  50.5   4.5 -20.5
[3,] -13.5 -11.5  9.0  4.5  -7.0 -13.5  -4.5  -15   8.5  53.5   8.0 -17.0
[4,] -11.5  -5.0 15.5  8.5   0.0  -4.5  -1.0  -14  10.5  56.0  10.0 -16.0
[5,] -11.0  -4.0 17.0 12.0   5.0   3.0   1.0  -13  12.0  56.0  10.0 -16.0

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

$conf
        [,1]    [,2]  [,3] [,4]    [,5]    [,6]    [,7]   [,8]   [,9]  [,10]
[1,] -17.055 -22.955 -5.22 0.55 -15.295 -22.585 -10.425 -18.16 -1.375 49.155
[2,]  -9.945  -0.045 23.22 8.45   1.295  -4.415   1.425 -11.84 18.375 57.845
      [,11]   [,12]
[1,]  3.655 -20.555
[2,] 12.345 -13.445

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(-17, -16, -13.5, -11.5, -11, -22, -19.5, -11.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/65oxp1289754041.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,] 499.0 469.0 502.0 537.0  586
[2,] 529.0 491.5 517.0 555.0  586
[3,] 556.0 512.0 530.5 569.5  586
[4,] 565.5 527.5 536.5 580.5  586
[5,] 599.0 542.0 537.0 604.0  586

$n
[1] 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4] [,5]
[1,] 539.3521 495.5802 521.6059 557.8693  586
[2,] 572.6479 528.4198 539.3941 581.1307  586

$out
[1] 587 597

$group
[1] 3 3

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

Warning message:
In bxp(list(stats = c(499, 529, 556, 565.5, 599, 469, 491.5, 512,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/75oxp1289754041.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,] 510.250 517.50 513.8750
[2,] 528.625 532.75 530.6875
[3,] 540.875 538.75 539.8125
[4,] 555.850 558.25 554.9375
[5,] 575.000 581.00 578.0000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 528.4575 527.1193 528.7519
[2,] 553.2925 550.3807 550.8731

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(510.25, 528.625, 540.875, 555.85, 575, 517.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> try(system("convert tmp/126z11289754041.ps tmp/126z11289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/226z11289754041.ps tmp/226z11289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/326z11289754041.ps tmp/326z11289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/4dfgm1289754041.ps tmp/4dfgm1289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/5dfgm1289754041.ps tmp/5dfgm1289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/65oxp1289754041.ps tmp/65oxp1289754041.png",intern=TRUE))
character(0)
> try(system("convert tmp/75oxp1289754041.ps tmp/75oxp1289754041.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.179   0.935   1.573