R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(299.9,339.2,374.2,393.5,389.2,381.7,375.2,369,357.4,352.1,346.5,342.9,340.3,328.3,322.9,314.3,308.9,294,285.6,281.2,280.3,278.8,274.5,270.4,263.4,259.9,258,262.7,284.7,311.3,322.1,327,331.3,333.3,321.4,327,320,314.7,316.7,314.4,321.3,318.2,307.2,301.3,287.5,277.7,274.4,258.8,253.3,251,248.4,249.5,246.1,244.5,243.6,244,240.8,249.8,248,259.4,260.5,260.8,261.3,259.5,256.6,257.9,256.5,254.2,253.3,253.8,255.5,257.1,257.3,253.2,252.8,252,250.7,252.2,250,251,253.4,251.2,255.6,261.1,258.9,259.9,261.2,264.7,267.1,266.4,267.7,268.6,267.5,268.5,268.5)
> 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] 95
> (np <- floor(n / par1))
[1] 7
> 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] 8 8 8 8 8 8 8 8 8 8 8 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 299.9 340.3 263.4 320.0 253.3 260.5 257.3 258.9
 [2,] 339.2 328.3 259.9 314.7 251.0 260.8 253.2 259.9
 [3,] 374.2 322.9 258.0 316.7 248.4 261.3 252.8 261.2
 [4,] 393.5 314.3 262.7 314.4 249.5 259.5 252.0 264.7
 [5,] 389.2 308.9 284.7 321.3 246.1 256.6 250.7 267.1
 [6,] 381.7 294.0 311.3 318.2 244.5 257.9 252.2 266.4
 [7,] 375.2 285.6 322.1 307.2 243.6 256.5 250.0 267.7
 [8,] 369.0 281.2 327.0 301.3 244.0 254.2 251.0 268.6
 [9,] 357.4 280.3 331.3 287.5 240.8 253.3 253.4 267.5
[10,] 352.1 278.8 333.3 277.7 249.8 253.8 251.2 268.5
[11,] 346.5 274.5 321.4 274.4 248.0 255.5 255.6 268.5
[12,] 342.9 270.4 327.0 258.8 259.4 257.1 261.1    NA
> darr
       [,1]  [,2]  [,3]  [,4] [,5] [,6] [,7] [,8]
 [1,]  39.3 -12.0  -3.5  -5.3 -2.3  0.3 -4.1  1.0
 [2,]  35.0  -5.4  -1.9   2.0 -2.6  0.5 -0.4  1.3
 [3,]  19.3  -8.6   4.7  -2.3  1.1 -1.8 -0.8  3.5
 [4,]  -4.3  -5.4  22.0   6.9 -3.4 -2.9 -1.3  2.4
 [5,]  -7.5 -14.9  26.6  -3.1 -1.6  1.3  1.5 -0.7
 [6,]  -6.5  -8.4  10.8 -11.0 -0.9 -1.4 -2.2  1.3
 [7,]  -6.2  -4.4   4.9  -5.9  0.4 -2.3  1.0  0.9
 [8,] -11.6  -0.9   4.3 -13.8 -3.2 -0.9  2.4 -1.1
 [9,]  -5.3  -1.5   2.0  -9.8  9.0  0.5 -2.2  1.0
[10,]  -5.6  -4.3 -11.9  -3.3 -1.8  1.7  4.4  0.0
[11,]  -3.6  -4.1   5.6 -15.6 11.4  1.6  5.5   NA
[12,]  -2.6  -7.0  -7.0  -5.5  1.1  0.2 -2.2   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/wessaorg/rcomp/tmp/1b3vo1321540435.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/wessaorg/rcomp/tmp/2fbfw1321540435.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/wessaorg/rcomp/tmp/33nmk1321540435.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/wessaorg/rcomp/tmp/43fw61321540435.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]   [,9]  [,10]
[1,] 253.30 251.00 248.40 249.50 246.10 244.50 243.60 244.00 240.80 249.80
[2,] 258.10 256.55 255.40 255.75 253.65 255.05 253.25 252.60 253.35 252.50
[3,] 261.95 260.35 261.25 263.70 275.90 280.20 276.65 274.90 273.90 273.10
[4,] 309.95 321.50 319.80 314.35 315.10 314.75 314.65 314.15 309.40 306.05
[5,] 340.30 339.20 374.20 393.50 389.20 381.70 375.20 369.00 357.40 352.10
      [,11] [,12]
[1,] 248.00 257.1
[2,] 255.55 259.1
[3,] 271.45 261.1
[4,] 297.95 298.7
[5,] 346.50 342.9

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 7

$conf
         [,1]    [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 232.9858 224.068 225.2752 230.9652 241.5731 246.8507 242.3511 240.5173
[2,] 290.9142 296.632 297.2248 296.4348 310.2269 313.5493 310.9489 309.2827
         [,9]    [,10]    [,11]    [,12]
[1,] 242.5897 243.1862 247.7648 237.4515
[2,] 305.2103 303.0138 295.1352 284.7485

$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(253.3, 258.1, 261.95, 309.95, 340.3, 251,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5mzjg1321540435.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]   [,8]  [,9]  [,10]  [,11]
[1,] -12.00 -5.40 -8.60 -5.40 -14.90 -11.00 -6.20 -13.80 -9.80 -11.90 -15.60
[2,]  -4.70 -2.25 -2.05 -3.85  -5.30  -7.45 -5.15  -7.40 -3.75  -4.95  -3.85
[3,]  -2.90  0.05  0.15 -2.10  -1.15  -1.80 -0.95  -1.00 -0.50  -2.55   1.60
[4,]   0.65  1.65  4.10  4.65   1.40   0.20  0.95   0.75  1.50   0.85   5.55
[5,]   1.00  2.00  4.70  6.90   1.50  10.80  4.90   4.30  9.00   4.40  11.40
     [,12]
[1,] -7.00
[2,] -6.25
[3,] -2.60
[4,] -1.00
[5,]  1.10

$n
 [1] 8 8 8 8 8 8 8 8 8 8 7 7

$conf
            [,1]      [,2]      [,3]      [,4]      [,5]    [,6]      [,7]
[1,] -5.88858681 -2.128596 -3.285478 -6.848222 -4.892716 -6.0734 -4.357548
[2,]  0.08858681  2.228596  3.585478  2.648222  2.592716  2.4734  2.457548
          [,8]      [,9]      [,10]     [,11]      [,12]
[1,] -5.552707 -3.432725 -5.7899633 -4.013528 -5.7352153
[2,]  3.552707  2.432725  0.6899633  7.213528  0.5352153

$out
[1] 39.3 35.0 19.3 22.0 26.6

$group
[1] 1 2 3 4 5

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

Warning message:
In bxp(list(stats = c(-12, -4.70000000000002, -2.90000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6ksw71321540435.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]  [,7]   [,8]
[1,] 299.90 270.40 258.00 258.80 240.80 253.30 250.0 258.90
[2,] 344.70 279.55 263.05 282.60 244.25 254.85 251.1 262.95
[3,] 363.20 289.80 316.35 310.80 248.20 256.85 252.5 267.10
[4,] 378.45 318.60 327.00 317.45 250.40 260.00 254.5 268.10
[5,] 393.50 340.30 333.30 321.30 259.40 261.30 257.3 268.60

$n
[1] 12 12 12 12 12 12 12 11

$conf
         [,1]    [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 347.8064 271.989 287.182 294.9047 245.3949 254.5011 250.9492 264.6466
[2,] 378.5936 307.611 345.518 326.6953 251.0051 259.1989 254.0508 269.5534

$out
[1] 261.1

$group
[1] 7

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

Warning message:
In bxp(list(stats = c(299.9, 344.7, 363.2, 378.45, 393.5, 270.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/7esqd1321540435.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,] 280.5500 260.350 270.9000
[2,] 282.7679 261.600 277.4062
[3,] 285.4375 272.275 282.3062
[4,] 288.6562 275.400 285.3625
[5,] 290.7750 280.200 288.1625

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 282.7518 265.9807 278.6774
[2,] 288.1232 278.5693 285.9351

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(280.55, 282.767857142857, 285.4375, 288.65625,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1b3vo1321540435.ps tmp/1b3vo1321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/2fbfw1321540435.ps tmp/2fbfw1321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/33nmk1321540435.ps tmp/33nmk1321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/43fw61321540435.ps tmp/43fw61321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/5mzjg1321540435.ps tmp/5mzjg1321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ksw71321540435.ps tmp/6ksw71321540435.png",intern=TRUE))
character(0)
> try(system("convert tmp/7esqd1321540435.ps tmp/7esqd1321540435.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.568   0.339   1.940