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(104.9,110.9,104.8,94.1,95.8,99.3,101.1,104.0,99.0,105.4,107.1,110.7,117.1,118.7,126.5,127.5,134.6,131.8,135.9,142.7,141.7,153.4,145.0,137.7,148.3,152.2,169.4,168.6,161.1,174.1,179.0,190.6,190.0,181.6,174.8,180.5,196.8,193.8,197.0,216.3,221.4,217.9,229.7,227.4,204.2,196.6,198.8,207.5,190.7,201.6,210.5,223.5,223.8,231.2,244.0,234.7,250.2,265.7,287.6,283.3,295.4,312.3,333.8,347.7,383.2,407.1,413.6,362.7,321.9,239.4,191.0,159.7,166.7)
> par1 = '4'
> #'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] 73
> (np <- floor(n / par1))
[1] 18
> 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] 19 18 18 18
> arr
      [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
[1,] 104.9  95.8  99.0 117.1 134.6 141.7 148.3 161.1 190.0 196.8 221.4 204.2
[2,] 110.9  99.3 105.4 118.7 131.8 153.4 152.2 174.1 181.6 193.8 217.9 196.6
[3,] 104.8 101.1 107.1 126.5 135.9 145.0 169.4 179.0 174.8 197.0 229.7 198.8
[4,]  94.1 104.0 110.7 127.5 142.7 137.7 168.6 190.6 180.5 216.3 227.4 207.5
     [,13] [,14] [,15] [,16] [,17] [,18] [,19]
[1,] 190.7 223.8 250.2 295.4 383.2 321.9 166.7
[2,] 201.6 231.2 265.7 312.3 407.1 239.4    NA
[3,] 210.5 244.0 287.6 333.8 413.6 191.0    NA
[4,] 223.5 234.7 283.3 347.7 362.7 159.7    NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
[1,]   6.0  3.5  6.4  1.6 -2.8 11.7  3.9 13.0 -8.4  -3.0  -3.5  -7.6  10.9
[2,]  -6.1  1.8  1.7  7.8  4.1 -8.4 17.2  4.9 -6.8   3.2  11.8   2.2   8.9
[3,] -10.7  2.9  3.6  1.0  6.8 -7.3 -0.8 11.6  5.7  19.3  -2.3   8.7  13.0
[4,]   1.7 -5.0  6.4  7.1 -1.0 10.6 -7.5 -0.6 16.3   5.1 -23.2 -16.8   0.3
     [,14] [,15] [,16] [,17] [,18] [,19]
[1,]   7.4  15.5  16.9  23.9 -82.5    NA
[2,]  12.8  21.9  21.5   6.5 -48.4    NA
[3,]  -9.3  -4.3  13.9 -50.9 -31.3    NA
[4,]  15.5  12.1  35.5 -40.8   7.0    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/1w1gv1242290905.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/2tl8q1242290905.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/3n1ut1242290905.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/4oo6z1242290905.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]
[1,]  95.80  99.3 101.1  94.10
[2,] 138.15 131.8 135.9 137.70
[3,] 190.00 187.7 185.0 185.55
[4,] 222.60 231.2 229.7 227.40
[5,] 321.90 312.3 333.8 347.70

$n
[1] 19 18 18 18

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 159.3888 150.6825 150.0680 152.1449
[2,] 220.6112 224.7175 219.9320 218.9551

$out
[1] 383.2 407.1 413.6 362.7

$group
[1] 1 2 3 4

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

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5eeol1242290905.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]
[1,] -8.40 -8.4 -31.30 -23.2
[2,] -3.00  1.7  -7.30  -5.0
[3,]  4.95  4.5   1.95   3.4
[4,] 11.70 11.8   8.70  10.6
[5,] 23.90 21.9  19.30  16.3

$n
[1] 18 18 18 18

$conf
           [,1]      [,2]      [,3]      [,4]
[1,] -0.5244207 0.7386633 -4.008553 -2.409589
[2,] 10.4244207 8.2613367  7.908553  9.209589

$out
[1] -82.5 -48.4 -50.9  35.5 -40.8

$group
[1] 1 2 3 4 4

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

Warning message:
In bxp(list(stats = c(-8.4, -3, 4.94999999999999, 11.7000000000000,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/674ta1242290905.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,]  94.10  95.80  99.00 117.1 131.80 137.70 148.30 161.10 174.80 193.80
[2,]  99.45  97.55 102.20 117.9 133.20 139.70 150.25 167.60 177.65 195.30
[3,] 104.85 100.20 106.25 122.6 135.25 143.35 160.40 176.55 181.05 196.90
[4,] 107.90 102.55 108.90 127.0 139.30 149.20 169.00 184.80 185.80 206.65
[5,] 110.90 104.00 110.70 127.5 142.70 153.40 169.40 190.60 190.00 216.30
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18] [,19]
[1,] 217.90 196.60 190.70 223.80 250.20 295.40 362.70 159.70 166.7
[2,] 219.65 197.70 196.15 227.50 257.95 303.85 372.95 175.35 166.7
[3,] 224.40 201.50 206.05 232.95 274.50 323.05 395.15 215.20 166.7
[4,] 228.55 205.85 217.00 239.35 285.45 340.75 410.35 280.65 166.7
[5,] 229.70 207.50 223.50 244.00 287.60 347.70 413.60 321.90 166.7

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

$conf
         [,1]   [,2]    [,3]    [,4]    [,5]    [,6]     [,7]    [,8]     [,9]
[1,]  98.1745  96.25 100.957 115.411 130.431 135.845 145.5875 162.962 174.6115
[2,] 111.5255 104.15 111.543 129.789 140.069 150.855 175.2125 190.138 187.4885
        [,10]   [,11]    [,12]    [,13]    [,14]   [,15]   [,16]   [,17]
[1,] 187.9335 217.369 195.0615 189.5785 223.5885 252.775 293.899 365.604
[2,] 205.8665 231.431 207.9385 222.5215 242.3115 296.225 352.201 424.696
       [,18] [,19]
[1,] 132.013 166.7
[2,] 298.387 166.7

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" NA  

Warning message:
In bxp(list(stats = c(94.1, 99.45, 104.85, 107.9, 110.9, 95.8, 97.55,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7ds8t1242290905.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,] 191.9368 185.000 180.3750
[2,] 192.9962 185.275 180.9563
[3,] 194.7833 186.625 181.9625
[4,] 196.3556 188.850 182.5375
[5,] 197.2000 190.000 182.6875

$n
[1] 4 4 4

$conf
         [,1]     [,2]     [,3]
[1,] 192.1294 183.8008 180.7133
[2,] 197.4372 189.4492 183.2117

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(191.936842105263, 192.996198830409, 194.783333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1w1gv1242290905.ps tmp/1w1gv1242290905.png")
> system("convert tmp/2tl8q1242290905.ps tmp/2tl8q1242290905.png")
> system("convert tmp/3n1ut1242290905.ps tmp/3n1ut1242290905.png")
> system("convert tmp/4oo6z1242290905.ps tmp/4oo6z1242290905.png")
> system("convert tmp/5eeol1242290905.ps tmp/5eeol1242290905.png")
> system("convert tmp/674ta1242290905.ps tmp/674ta1242290905.png")
> system("convert tmp/7ds8t1242290905.ps tmp/7ds8t1242290905.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.207   0.949   1.735