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(288.6,269.1,268.7,264.3,264.3,267.6,298.1,279.8,263.2,272.5,263.7,273.7,261.4,241.1,253.4,228.6,244.9,206.1,177,165.1,148.1,152.9,146.5,188,252,351.6,403,468.8,464,435.4,382.2,360.6,329.5,320.2,315,322.7,289.7,270.3,247.8,259.6,241,230,230.3,214,202.9,188.5,215.6,205.6,203.7,218.2,253,255.4,240.7,242.2,240.2,215.2,211.1,219.3)
> 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] 58
> (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 5 5 5 5 5 5 5 5 5 4 4
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] 288.6 261.4 252.0 289.7 203.7
 [2,] 269.1 241.1 351.6 270.3 218.2
 [3,] 268.7 253.4 403.0 247.8 253.0
 [4,] 264.3 228.6 468.8 259.6 255.4
 [5,] 264.3 244.9 464.0 241.0 240.7
 [6,] 267.6 206.1 435.4 230.0 242.2
 [7,] 298.1 177.0 382.2 230.3 240.2
 [8,] 279.8 165.1 360.6 214.0 215.2
 [9,] 263.2 148.1 329.5 202.9 211.1
[10,] 272.5 152.9 320.2 188.5 219.3
[11,] 263.7 146.5 315.0 215.6    NA
[12,] 273.7 188.0 322.7 205.6    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]
 [1,] -19.5 -20.3  99.6 -19.4  14.5
 [2,]  -0.4  12.3  51.4 -22.5  34.8
 [3,]  -4.4 -24.8  65.8  11.8   2.4
 [4,]   0.0  16.3  -4.8 -18.6 -14.7
 [5,]   3.3 -38.8 -28.6 -11.0   1.5
 [6,]  30.5 -29.1 -53.2   0.3  -2.0
 [7,] -18.3 -11.9 -21.6 -16.3 -25.0
 [8,] -16.6 -17.0 -31.1 -11.1  -4.1
 [9,]   9.3   4.8  -9.3 -14.4   8.2
[10,]  -8.8  -6.4  -5.2  27.1    NA
[11,]  10.0  41.5   7.7 -10.0    NA
[12,] -12.3  64.0 -33.0  -1.9    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/1mr2c1289851994.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/2mr2c1289851994.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/37s4s1289851995.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/47s4s1289851995.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,] 203.7 218.2 247.8 255.4 240.7 206.1 177.0 165.1 148.1 152.9 146.50 188.00
[2,] 252.0 241.1 253.0 255.4 241.0 230.0 230.3 214.0 202.9 188.5 181.05 196.80
[3,] 261.4 269.1 253.4 259.6 244.9 242.2 240.2 215.2 211.1 219.3 239.65 239.65
[4,] 288.6 270.3 268.7 264.3 264.3 267.6 298.1 279.8 263.2 272.5 289.35 298.20
[5,] 289.7 270.3 268.7 264.3 264.3 267.6 382.2 360.6 329.5 320.2 315.00 322.70

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 235.5385 248.4674 242.3064 253.3113 228.4363 215.6319 192.2927 168.7059
[2,] 287.2615 289.7326 264.4936 265.8887 261.3637 268.7681 288.1073 261.6941
         [,9]    [,10]   [,11]   [,12]
[1,] 168.4922 159.9458 154.093 159.544
[2,] 253.7078 278.6542 325.207 319.756

$out
[1] 351.6 403.0 228.6 468.8 464.0 435.4

$group
[1] 2 3 4 4 5 6

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

Warning message:
In bxp(list(stats = c(203.7, 252, 261.4, 288.6, 289.7, 218.2, 241.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/57s4s1289851995.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,] -20.3 -22.5 -24.8 -18.6 -38.8 -53.2 -25.0 -17.0 -14.4 -8.80 -10.00 -33.00
[2,] -19.5  -0.4  -4.4 -14.7 -28.6 -29.1 -21.6 -17.0  -9.3 -7.60  -1.15 -22.65
[3,] -19.4  12.3   2.4  -4.8 -11.0  -2.0 -18.3 -16.6   4.8 -5.80   8.85  -7.10
[4,]  14.5  34.8  11.8   0.0   1.5   0.3 -16.3 -11.1   8.2 10.95  25.75  31.05
[5,]  14.5  51.4  11.8  16.3   3.3  30.5 -11.9  -4.1   9.3 27.10  41.50  64.00

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

$conf
           [,1]      [,2]     [,3]       [,4]      [,5]      [,6]      [,7]
[1,] -43.424314 -12.57223 -9.04688 -15.186983 -32.26858 -22.77397 -22.04497
[2,]   4.624314  37.17223 13.84688   5.586983  10.26858  18.77397 -14.55503
          [,8]      [,9]    [,10]   [,11]   [,12]
[1,] -20.76893 -7.565456 -20.4545 -12.401 -49.523
[2,] -12.43107 17.165456   8.8545  30.101  35.323

$out
[1]  99.6  65.8 -31.1

$group
[1] 1 3 8

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

Warning message:
In bxp(list(stats = c(-20.3, -19.5, -19.4000000000000, 14.5, 14.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6zk3v1289851995.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,] 263.20 146.50 252.00 188.50 203.70
[2,] 264.30 159.00 321.45 209.80 215.20
[3,] 268.90 197.05 356.10 230.15 229.75
[4,] 276.75 243.00 419.20 253.70 242.20
[5,] 288.60 261.40 468.80 289.70 255.40

$n
[1] 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]
[1,] 263.2215 158.7370 311.5156 210.1269 216.2597
[2,] 274.5785 235.3630 400.6844 250.1731 243.2403

$out
[1] 298.1

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(263.2, 264.3, 268.9, 276.75, 288.6, 146.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7zk3v1289851995.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,] 230.68 211.100 230.500
[2,] 241.07 229.475 240.500
[3,] 262.32 241.200 250.725
[4,] 280.72 256.500 260.350
[5,] 295.34 269.100 270.300

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 244.2354 228.8737 241.6713
[2,] 280.4046 253.5263 259.7787

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(230.68, 241.07, 262.32, 280.72, 295.34, 211.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1mr2c1289851994.ps tmp/1mr2c1289851994.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mr2c1289851994.ps tmp/2mr2c1289851994.png",intern=TRUE))
character(0)
> try(system("convert tmp/37s4s1289851995.ps tmp/37s4s1289851995.png",intern=TRUE))
character(0)
> try(system("convert tmp/47s4s1289851995.ps tmp/47s4s1289851995.png",intern=TRUE))
character(0)
> try(system("convert tmp/57s4s1289851995.ps tmp/57s4s1289851995.png",intern=TRUE))
character(0)
> try(system("convert tmp/6zk3v1289851995.ps tmp/6zk3v1289851995.png",intern=TRUE))
character(0)
> try(system("convert tmp/7zk3v1289851995.ps tmp/7zk3v1289851995.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.912   1.469   2.256