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(395.3,395.1,403.5,403.3,405.7,406.7,407.2,412.4,415.9,414,411.8,409.9,412.4,415.9,416.3,417.2,421.8,421.4,415.1,412.4,411.8,408.8,404.5,402.5,409.4,410.7,413.4,415.2,417.7,417.8,417.9,418.4,418.2,416.6,418.9,421,423.5,432.3,432.3,428.6,426.7,427.3,428.5,437,442,444.9,441.4,440.3,447.1,455.3,478.6,486.5,487.8,485.9,483.8,488.4,494,493.6,487.3,482.1,484.2,496.8,501.1,499.8,495.5,498.1,503.8,516.2,526.1,527.1,525.1,528.9,540.1,549,556,568.9,589.1,590.3,603.3,638.8,643,656.7,656.1,654.1,659.9,662.1,669.2,673.1,678.3,677.4,678.5,672.4,665.3,667.9,672.1,662.5,682.3,692.1,702.7,721.4,733.2,747.7,737.6,729.3,706.1,674.3,659,645.7)
> 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] 108
> (np <- floor(n / par1))
[1] 9
> 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] 9 9 9 9 9 9 9 9 9 9 9 9
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10]
 [1,] 395.3 412.4 409.4 423.5 447.1 484.2 540.1 659.9 682.3    NA
 [2,] 395.1 415.9 410.7 432.3 455.3 496.8 549.0 662.1 692.1    NA
 [3,] 403.5 416.3 413.4 432.3 478.6 501.1 556.0 669.2 702.7    NA
 [4,] 403.3 417.2 415.2 428.6 486.5 499.8 568.9 673.1 721.4    NA
 [5,] 405.7 421.8 417.7 426.7 487.8 495.5 589.1 678.3 733.2    NA
 [6,] 406.7 421.4 417.8 427.3 485.9 498.1 590.3 677.4 747.7    NA
 [7,] 407.2 415.1 417.9 428.5 483.8 503.8 603.3 678.5 737.6    NA
 [8,] 412.4 412.4 418.4 437.0 488.4 516.2 638.8 672.4 729.3    NA
 [9,] 415.9 411.8 418.2 442.0 494.0 526.1 643.0 665.3 706.1    NA
[10,] 414.0 408.8 416.6 444.9 493.6 527.1 656.7 667.9 674.3    NA
[11,] 411.8 404.5 418.9 441.4 487.3 525.1 656.1 672.1 659.0    NA
[12,] 409.9 402.5 421.0 440.3 482.1 528.9 654.1 662.5 645.7    NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]  [,9] [,10]
 [1,] -0.2  3.5  1.3  8.8  8.2 12.6  8.9  2.2   9.8    NA
 [2,]  8.4  0.4  2.7  0.0 23.3  4.3  7.0  7.1  10.6    NA
 [3,] -0.2  0.9  1.8 -3.7  7.9 -1.3 12.9  3.9  18.7    NA
 [4,]  2.4  4.6  2.5 -1.9  1.3 -4.3 20.2  5.2  11.8    NA
 [5,]  1.0 -0.4  0.1  0.6 -1.9  2.6  1.2 -0.9  14.5    NA
 [6,]  0.5 -6.3  0.1  1.2 -2.1  5.7 13.0  1.1 -10.1    NA
 [7,]  5.2 -2.7  0.5  8.5  4.6 12.4 35.5 -6.1  -8.3    NA
 [8,]  3.5 -0.6 -0.2  5.0  5.6  9.9  4.2 -7.1 -23.2    NA
 [9,] -1.9 -3.0 -1.6  2.9 -0.4  1.0 13.7  2.6 -31.8    NA
[10,] -2.2 -4.3  2.3 -3.5 -6.3 -2.0 -0.6  4.2 -15.3    NA
[11,] -1.9 -2.0  2.1 -1.1 -5.2  3.8 -2.0 -9.6 -13.3    NA
[12,]  2.5  6.9  2.5  6.8  2.1 11.2  5.8 19.8    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/1f98h1257850272.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/2oyfk1257850272.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/3fyqp1257850273.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/4ujgq1257850273.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,] 395.3 395.1 403.5 403.3 405.7 406.7 407.2 412.4 411.8 408.8 404.5 402.5
[2,] 412.4 415.9 416.3 417.2 421.8 421.4 417.9 418.4 418.2 416.6 418.9 421.0
[3,] 447.1 455.3 478.6 486.5 487.8 485.9 483.8 488.4 494.0 493.6 487.3 482.1
[4,] 540.1 549.0 556.0 568.9 589.1 590.3 603.3 638.8 643.0 656.7 656.1 645.7
[5,] 682.3 692.1 702.7 721.4 733.2 747.7 737.6 729.3 706.1 674.3 672.1 662.5

$n
 [1] 9 9 9 9 9 9 9 9 9 9 9 9

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]    [,7]     [,8]
[1,] 379.8447 385.2007 405.0247 406.6047 399.6887 396.946 386.156 372.3227
[2,] 514.3553 525.3993 552.1753 566.3953 575.9113 574.854 581.444 604.4773
         [,9]    [,10]    [,11]   [,12]
[1,] 375.6053 367.1473 362.3747 363.758
[2,] 612.3947 620.0527 612.2253 600.442

$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(395.3, 412.4, 447.1, 540.1, 682.3, 395.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/543ah1257850273.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,] -0.2  0.0 -3.7 -4.3 -1.9 -6.3 -8.3 -7.1 -3.0  -6.3  -9.6  2.10
[2,]  2.2  2.7 -0.2  1.3 -0.4 -2.1 -2.7 -0.6 -1.9  -4.3  -5.2  2.50
[3,]  8.2  7.0  1.8  2.5  0.6  0.5  4.6  3.5 -0.4  -2.2  -2.0  6.30
[4,]  8.9  8.4  7.9  5.2  1.2  1.2  8.5  5.0  2.6  -0.6  -1.1  9.05
[5,] 12.6 10.6 18.7  5.2  2.6  5.7 12.4  9.9  2.9   4.2   3.8 11.20

$n
 [1] 9 9 9 9 9 9 9 9 9 9 9 8

$conf
          [,1]   [,2]   [,3]  [,4]       [,5]   [,6]      [,7]      [,8]  [,9]
[1,]  4.671333  3.998 -2.466 0.446 -0.2426667 -1.238 -1.298667 0.5506667 -2.77
[2,] 11.728667 10.002  6.066 4.554  1.4426667  2.238 10.498667 6.4493333  1.97
          [,10]      [,11]    [,12]
[1,] -4.1486667 -4.1593333 2.641076
[2,] -0.2513333  0.1593333 9.958924

$out
 [1]  23.3  20.2  11.8  14.5  13.0 -10.1  35.5 -23.2  13.7 -31.8 -15.3 -13.3
[13]  19.8

$group
 [1]  2  4  4  5  6  6  7  8  9  9 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.199999999999989, 2.20000000000005, 8.19999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6pif51257850273.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]   [,7]   [,8]   [,9] [,10]
[1,] 395.10 402.50 409.40 423.50 478.60 484.20 540.10 659.90 645.70    NA
[2,] 403.40 410.30 414.30 427.90 480.35 497.45 562.45 663.90 678.30    NA
[3,] 406.95 413.75 417.75 432.30 486.20 502.45 596.80 670.65 704.40    NA
[4,] 412.10 416.75 418.30 440.85 488.10 525.60 648.55 675.25 731.25    NA
[5,] 415.90 421.80 421.00 444.90 494.00 528.90 656.70 678.50 747.70    NA

$n
 [1] 12 12 12 12 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 402.9819 410.8081 415.9256 426.3934 482.6652 489.6106 557.5292 665.4732
[2,] 410.9181 416.6919 419.5744 438.2066 489.7348 515.2894 636.0708 675.8268
         [,9] [,10]
[1,] 680.2491    NA
[2,] 728.5509    NA

$out
[1] 447.1 455.3

$group
[1] 5 5

$names
 [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA 

Warning message:
In bxp(list(stats = c(395.1, 403.4, 406.95, 412.1, 415.9, 402.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/79a8m1257850273.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,] 494.9111 478.60 476.250
[2,] 510.3944 480.35 489.600
[3,] 518.2444 486.20 508.225
[4,] 521.1167 488.10 531.975
[5,] 525.0333 494.00 537.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 513.3540 482.6652 488.8975
[2,] 523.1349 489.7348 527.5525

$out
[1] 447.1 455.3

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(494.911111111111, 510.394444444444, 518.244444444445,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1f98h1257850272.ps tmp/1f98h1257850272.png")
> system("convert tmp/2oyfk1257850272.ps tmp/2oyfk1257850272.png")
> system("convert tmp/3fyqp1257850273.ps tmp/3fyqp1257850273.png")
> system("convert tmp/4ujgq1257850273.ps tmp/4ujgq1257850273.png")
> system("convert tmp/543ah1257850273.ps tmp/543ah1257850273.png")
> system("convert tmp/6pif51257850273.ps tmp/6pif51257850273.png")
> system("convert tmp/79a8m1257850273.ps tmp/79a8m1257850273.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.205   0.955   1.511