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.

  Natural language support but running in an English locale

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/freestat/rcomp/tmp/1nt9w1289835454.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/freestat/rcomp/tmp/2nt9w1289835454.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/freestat/rcomp/tmp/3f28h1289835454.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/freestat/rcomp/tmp/4f28h1289835454.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/freestat/rcomp/tmp/5f28h1289835454.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/freestat/rcomp/tmp/6qtp11289835454.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/freestat/rcomp/tmp/7qtp11289835454.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/1nt9w1289835454.ps tmp/1nt9w1289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/2nt9w1289835454.ps tmp/2nt9w1289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/3f28h1289835454.ps tmp/3f28h1289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/4f28h1289835454.ps tmp/4f28h1289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/5f28h1289835454.ps tmp/5f28h1289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/6qtp11289835454.ps tmp/6qtp11289835454.png",intern=TRUE))
character(0)
> try(system("convert tmp/7qtp11289835454.ps tmp/7qtp11289835454.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.942   1.463   2.205