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.

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(74.8,93.1,103.9,83.9,77.7,141.5,58.9,75.3,108.4,91,84.6,179.8,85.6,76.4,109.7,99.1,86.7,111.4,78.4,76.7,114.2,99.7,94.2,173.5,83.1,88.9,132,122.1,105.1,133.7,63.6,112.7,120.5,112,126.2,209.2,91,116.7,137.6,108.1,136.6,152.3,114.3,120.7,131.8,129.4,187.5,189.5,109.2,158.1,176.2,125.5,155,170.3,99.4,139.2,169.6,136.1,168.2,318.6,154.1,161.4,183.4,167.2,205.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] 65
> (np <- floor(n / par1))
[1] 5
> 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] 6 6 6 6 6 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]
 [1,]  74.8  85.6  83.1  91.0 109.2 154.1
 [2,]  93.1  76.4  88.9 116.7 158.1 161.4
 [3,] 103.9 109.7 132.0 137.6 176.2 183.4
 [4,]  83.9  99.1 122.1 108.1 125.5 167.2
 [5,]  77.7  86.7 105.1 136.6 155.0 205.3
 [6,] 141.5 111.4 133.7 152.3 170.3    NA
 [7,]  58.9  78.4  63.6 114.3  99.4    NA
 [8,]  75.3  76.7 112.7 120.7 139.2    NA
 [9,] 108.4 114.2 120.5 131.8 169.6    NA
[10,]  91.0  99.7 112.0 129.4 136.1    NA
[11,]  84.6  94.2 126.2 187.5 168.2    NA
[12,] 179.8 173.5 209.2 189.5 318.6    NA
> darr
       [,1]  [,2]   [,3]  [,4]   [,5]  [,6]
 [1,]  18.3  -9.2    5.8  25.7   48.9   7.3
 [2,]  10.8  33.3   43.1  20.9   18.1  22.0
 [3,] -20.0 -10.6   -9.9 -29.5  -50.7 -16.2
 [4,]  -6.2 -12.4  -17.0  28.5   29.5  38.1
 [5,]  63.8  24.7   28.6  15.7   15.3    NA
 [6,] -82.6 -33.0  -70.1 -38.0  -70.9    NA
 [7,]  16.4  -1.7   49.1   6.4   39.8    NA
 [8,]  33.1  37.5    7.8  11.1   30.4    NA
 [9,] -17.4 -14.5   -8.5  -2.4  -33.5    NA
[10,]  -6.4  -5.5   14.2  58.1   32.1    NA
[11,]  95.2  79.3   83.0   2.0  150.4    NA
[12,] -94.2 -90.4 -118.2 -80.3 -164.5    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/1im9c1228353275.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/2fc2a1228353275.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/3zgrg1228353275.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/4hqxy1228353275.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,]  74.8  76.4 103.9  83.9  77.70 111.4  58.9  75.3 108.4  91.0  84.6 173.5
[2,]  83.1  88.9 109.7  99.1  86.70 133.7  63.6  76.7 114.2  99.7  94.2 179.8
[3,]  88.3 104.9 134.8 115.1 120.85 141.5  78.4 112.7 120.5 112.0 126.2 189.5
[4,] 109.2 158.1 176.2 125.5 155.00 152.3  99.4 120.7 131.8 129.4 168.2 209.2
[5,] 109.2 161.4 183.4 125.5 205.30 170.3 114.3 139.2 131.8 136.1 187.5 209.2

$n
 [1] 6 6 6 6 6 5 5 5 5 5 5 5

$conf
          [,1]      [,2]      [,3]      [,4]     [,5]     [,6]      [,7]
[1,]  71.46466  60.26376  91.90535  98.07115  76.7943 128.3573  53.10381
[2,] 105.13534 149.53624 177.69465 132.12885 164.9057 154.6427 103.69619
          [,8]     [,9]     [,10]     [,11]    [,12]
[1,]  81.60971 108.0639  91.01405  73.91179 168.7260
[2,] 143.79029 132.9361 132.98595 178.48821 210.2740

$out
[1] 154.1 167.2 169.6 318.6

$group
[1]  1  4  9 12

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

Warning message:
In bxp(list(stats = c(74.8, 83.1, 88.3, 109.2, 109.2, 76.4, 88.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/53foy1228353275.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,] -9.2 10.80 -50.7 -17.00 15.3 -82.6 -1.7  7.8 -17.4  -6.4  79.3 -118.2
[2,]  5.8 18.10 -29.5 -12.40 15.7 -70.9  6.4 11.1 -17.4  -5.5  79.3 -118.2
[3,] 12.8 21.45 -18.1  11.15 24.7 -70.1 16.4 30.4 -14.5  14.2  83.0  -94.2
[4,] 25.7 33.30 -10.6  29.50 28.6 -38.0 39.8 33.1  -8.5  32.1  95.2  -90.4
[5,] 48.9 43.10  -9.9  38.10 28.6 -33.0 49.1 37.5  -2.4  58.1  95.2  -80.3

$n
 [1] 6 6 6 6 5 5 5 5 5 5 5 5

$conf
            [,1]     [,2]      [,3]      [,4]     [,5]      [,6]      [,7]
[1,] -0.03614275 11.64551 -30.29111 -15.87685 15.58489 -93.34706 -7.200356
[2,] 25.63614275 31.25449  -5.90889  38.17685 33.81511 -46.85294 40.000356
         [,8]       [,9]     [,10]   [,11]      [,12]
[1,] 14.85486 -20.788718 -12.36807 71.7651 -113.84341
[2,] 45.94514  -8.211282  40.76807 94.2349  -74.55659

$out
[1]   63.8  -33.5    2.0  150.4 -164.5

$group
[1]  5  9 11 11 12

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

Warning message:
In bxp(list(stats = c(-9.19999999999999, 5.80000000000001, 12.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/64lz41228353275.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]
[1,]  58.90  76.40  63.6  91.00  99.40 154.1
[2,]  76.50  82.00  97.0 115.50 130.80 161.4
[3,]  87.80  96.65 116.6 130.60 156.55 167.2
[4,] 106.15 110.55 129.1 144.95 169.95 183.4
[5,] 141.50 114.20 133.7 187.50 176.20 205.3

$n
[1] 12 12 12 12 12  5

$conf
          [,1]      [,2]     [,3]     [,4]     [,5]     [,6]
[1,]  74.27644  83.62815 101.9590 117.1677 138.6934 151.6549
[2,] 101.32356 109.67185 131.2410 144.0323 174.4066 182.7451

$out
[1] 179.8 173.5 209.2 189.5 318.6

$group
[1] 1 2 3 4 5

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

Warning message:
In bxp(list(stats = c(58.9, 76.5, 87.8, 106.15, 141.5, 76.4, 82,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7mnln1228353275.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,]  82.9200  78.40  81.5000
[2,] 109.2800 108.45 105.8500
[3,] 122.6917 117.80 119.8500
[4,] 136.3033 130.50 136.0562
[5,] 141.8400 141.50 143.0000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 110.3661 107.7428 106.0727
[2,] 135.0172 127.8572 133.6273

$out
[1] 214.12 189.50 194.50

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(82.92, 109.28, 122.691666666667, 136.303333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1im9c1228353275.ps tmp/1im9c1228353275.png")
> system("convert tmp/2fc2a1228353275.ps tmp/2fc2a1228353275.png")
> system("convert tmp/3zgrg1228353275.ps tmp/3zgrg1228353275.png")
> system("convert tmp/4hqxy1228353275.ps tmp/4hqxy1228353275.png")
> system("convert tmp/53foy1228353275.ps tmp/53foy1228353275.png")
> system("convert tmp/64lz41228353275.ps tmp/64lz41228353275.png")
> system("convert tmp/7mnln1228353275.ps tmp/7mnln1228353275.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.240   0.990   1.536