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(104.3,119.8,116.8,118.2,107.4,110.8,94.8,96.5,113.4,109.8,118.7,117.2,110.3,111.6,128.1,121.3,107.3,120.5,98.5,97.7,113.2,114.6,118.3,123.9,113.6,117.5,130.1,124.7,114.2,127.3,105.9,101.5,120.2,117.1,131.1,130,120.6,123.1,135.3,134.1,123.7,134.6,108.3,110.4,127.8,126.6,131.4,141.1,127,127.3,143.6,149.4,126.6,136.5,116,118,131.4,140.7,144.9,143.9,127.1)
> 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] 61
> (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 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]
 [1,] 104.3 110.3 113.6 120.6 127.0 127.1
 [2,] 119.8 111.6 117.5 123.1 127.3    NA
 [3,] 116.8 128.1 130.1 135.3 143.6    NA
 [4,] 118.2 121.3 124.7 134.1 149.4    NA
 [5,] 107.4 107.3 114.2 123.7 126.6    NA
 [6,] 110.8 120.5 127.3 134.6 136.5    NA
 [7,]  94.8  98.5 105.9 108.3 116.0    NA
 [8,]  96.5  97.7 101.5 110.4 118.0    NA
 [9,] 113.4 113.2 120.2 127.8 131.4    NA
[10,] 109.8 114.6 117.1 126.6 140.7    NA
[11,] 118.7 118.3 131.1 131.4 144.9    NA
[12,] 117.2 123.9 130.0 141.1 143.9    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6]
 [1,]  15.5   1.3   3.9   2.5   0.3   NA
 [2,]  -3.0  16.5  12.6  12.2  16.3   NA
 [3,]   1.4  -6.8  -5.4  -1.2   5.8   NA
 [4,] -10.8 -14.0 -10.5 -10.4 -22.8   NA
 [5,]   3.4  13.2  13.1  10.9   9.9   NA
 [6,] -16.0 -22.0 -21.4 -26.3 -20.5   NA
 [7,]   1.7  -0.8  -4.4   2.1   2.0   NA
 [8,]  16.9  15.5  18.7  17.4  13.4   NA
 [9,]  -3.6   1.4  -3.1  -1.2   9.3   NA
[10,]   8.9   3.7  14.0   4.8   4.2   NA
[11,]  -1.5   5.6  -1.1   9.7  -1.0   NA
[12,]  -6.9 -10.3  -9.4 -14.1 -16.8   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/17nln1229600067.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/2fmmz1229600067.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/3on1q1229600067.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/4ajpy1229600067.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,] 104.3 111.6 128.1 118.2 107.3 110.8  94.8  96.5 113.2 109.8 118.3 117.2
[2,] 110.3 117.5 128.1 121.3 107.4 120.5  98.5  97.7 113.4 114.6 118.7 123.9
[3,] 117.1 119.8 130.1 124.7 114.2 127.3 105.9 101.5 120.2 117.1 131.1 130.0
[4,] 127.0 123.1 135.3 134.1 123.7 134.6 108.3 110.4 127.8 126.6 131.4 141.1
[5,] 127.1 127.3 143.6 149.4 126.6 136.5 116.0 118.0 131.4 140.7 144.9 143.9

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]      [,7]      [,8]
[1,] 106.3280 115.8431 125.0125 115.6556 102.6825 117.3370  98.97534  92.52621
[2,] 127.8720 123.7569 135.1875 133.7444 125.7175 137.2630 112.82466 110.47379
        [,9]    [,10]    [,11]    [,12]
[1,] 110.025 108.6208 122.1262 117.8465
[2,] 130.375 125.5792 140.0738 142.1535

$out
[1] 116.8

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(104.3, 110.3, 117.1, 127, 127.1, 111.6, 117.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5x70z1229600067.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.3 12.2 -6.8 -14.0  9.9 -22.0 -4.4 13.4 -3.6   3.7  -1.5 -16.8
[2,]  1.3 12.2 -5.4 -14.0  9.9 -22.0 -0.8 15.5 -3.1   4.2  -1.1 -14.1
[3,]  2.5 12.6 -1.2 -10.8 10.9 -21.4  1.7 16.9 -1.2   4.8  -1.0 -10.3
[4,]  3.9 16.3  1.4 -10.5 13.1 -20.5  2.0 17.4  1.4   8.9   5.6  -9.4
[5,]  3.9 16.5  5.8 -10.4 13.2 -20.5  2.1 18.7  1.4  14.0   9.7  -6.9

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

$conf
          [,1]     [,2]      [,3]       [,4]      [,5]      [,6]       [,7]
[1,] 0.6628465  9.70295 -6.004863 -13.273091  8.638888 -22.45990 -0.2784729
[2,] 4.3371535 15.49705  3.604863  -8.326909 13.161112 -20.34010  3.6784729
         [,8]      [,9]    [,10]     [,11]      [,12]
[1,] 15.55746 -4.379689 1.478992 -5.734203 -13.621008
[2,] 18.24254  1.979689 8.121008  3.734203  -6.978992

$out
[1]  15.5  -3.0 -22.8   3.4 -16.0 -26.3   9.3

$group
[1] 1 2 4 5 6 6 9

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

Warning message:
In bxp(list(stats = c(0.299999999999997, 1.30000000000000, 2.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/67n9d1229600067.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,]  94.80  97.7 101.50 108.30 116.00 127.1
[2,] 105.85 108.8 113.90 121.85 126.80 127.1
[3,] 112.10 113.9 118.85 127.20 133.95 127.1
[4,] 117.70 120.9 128.65 134.35 143.75 127.1
[5,] 119.80 128.1 131.10 141.10 149.40 127.1

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]  [,6]
[1,] 106.6951 108.3811 112.1224 121.4987 126.219 127.1
[2,] 117.5049 119.4189 125.5776 132.9013 141.681 127.1

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(94.8, 105.85, 112.1, 117.7, 119.8, 97.7, 108.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/75qox1229600067.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,] 104.700 101.50 103.4000
[2,] 116.495 115.65 116.9062
[3,] 121.480 120.00 120.6000
[4,] 129.210 128.65 127.6250
[5,] 131.220 131.10 132.5000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 115.6806 114.0706 115.7111
[2,] 127.2794 125.9294 125.4889

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(104.7, 116.495, 121.48, 129.21, 131.22, 101.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/17nln1229600067.ps tmp/17nln1229600067.png")
> system("convert tmp/2fmmz1229600067.ps tmp/2fmmz1229600067.png")
> system("convert tmp/3on1q1229600067.ps tmp/3on1q1229600067.png")
> system("convert tmp/4ajpy1229600067.ps tmp/4ajpy1229600067.png")
> system("convert tmp/5x70z1229600067.ps tmp/5x70z1229600067.png")
> system("convert tmp/67n9d1229600067.ps tmp/67n9d1229600067.png")
> system("convert tmp/75qox1229600067.ps tmp/75qox1229600067.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.228   0.940   1.494