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(13.81,13.9,13.91,13.94,13.96,14.01,14.01,14.06,14.09,14.13,14.12,14.13,14.14,14.16,14.21,14.26,14.29,14.32,14.33,14.39,14.48,14.44,14.46,14.48,14.53,14.58,14.62,14.62,14.61,14.65,14.68,14.7,14.78,14.84,14.89,14.89,15.13,15.25,15.33,15.36,15.4,15.4,15.41,15.47,15.54,15.55,15.59,15.65,15.75,15.86,15.89,15.94,15.93,15.95,15.99,15.99,16.06,16.08,16.07,16.11,16.15,16.18,16.3,16.42,16.49,16.5,16.58,16.64,16.66,16.81,16.91,16.92,16.95,17.11,17.16,17.16,17.27,17.34,17.39,17.43,17.45,17.5,17.56,17.65,17.62)
> 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] 85
> (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 7 7 7 7 7 7 7 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 13.81 14.14 14.53 15.13 15.75 16.15 16.95 17.62
 [2,] 13.90 14.16 14.58 15.25 15.86 16.18 17.11    NA
 [3,] 13.91 14.21 14.62 15.33 15.89 16.30 17.16    NA
 [4,] 13.94 14.26 14.62 15.36 15.94 16.42 17.16    NA
 [5,] 13.96 14.29 14.61 15.40 15.93 16.49 17.27    NA
 [6,] 14.01 14.32 14.65 15.40 15.95 16.50 17.34    NA
 [7,] 14.01 14.33 14.68 15.41 15.99 16.58 17.39    NA
 [8,] 14.06 14.39 14.70 15.47 15.99 16.64 17.43    NA
 [9,] 14.09 14.48 14.78 15.54 16.06 16.66 17.45    NA
[10,] 14.13 14.44 14.84 15.55 16.08 16.81 17.50    NA
[11,] 14.12 14.46 14.89 15.59 16.07 16.91 17.56    NA
[12,] 14.13 14.48 14.89 15.65 16.11 16.92 17.65    NA
> darr
       [,1]  [,2]  [,3] [,4]  [,5] [,6]  [,7] [,8]
 [1,]  0.09  0.02  0.05 0.12  0.11 0.03  0.16   NA
 [2,]  0.01  0.05  0.04 0.08  0.03 0.12  0.05   NA
 [3,]  0.03  0.05  0.00 0.03  0.05 0.12  0.00   NA
 [4,]  0.02  0.03 -0.01 0.04 -0.01 0.07  0.11   NA
 [5,]  0.05  0.03  0.04 0.00  0.02 0.01  0.07   NA
 [6,]  0.00  0.01  0.03 0.01  0.04 0.08  0.05   NA
 [7,]  0.05  0.06  0.02 0.06  0.00 0.06  0.04   NA
 [8,]  0.03  0.09  0.08 0.07  0.07 0.02  0.02   NA
 [9,]  0.04 -0.04  0.06 0.01  0.02 0.15  0.05   NA
[10,] -0.01  0.02  0.05 0.04 -0.01 0.10  0.06   NA
[11,]  0.01  0.02  0.00 0.06  0.04 0.01  0.09   NA
[12,]  0.01  0.05  0.24 0.10  0.04 0.03 -0.03   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/1zqq81289671289.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/2zqq81289671289.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/3zqq81289671289.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/4sh7b1289671289.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]
[1,] 13.810 13.90 13.910 13.94 13.96 14.010 14.010 14.060 14.09 14.130 14.120
[2,] 14.335 14.37 14.415 14.44 14.45 14.485 14.505 14.545 14.63 14.640 14.675
[3,] 15.440 15.25 15.330 15.36 15.40 15.400 15.410 15.470 15.54 15.550 15.590
[4,] 16.550 16.02 16.095 16.18 16.21 16.225 16.285 16.315 16.36 16.445 16.490
[5,] 17.620 17.11 17.160 17.16 17.27 17.340 17.390 17.430 17.45 17.500 17.560
      [,12]
[1,] 14.130
[2,] 14.685
[3,] 15.650
[4,] 16.515
[5,] 17.650

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]    [,6]     [,7]     [,8]
[1,] 14.20267 14.26465 14.32673 14.3209 14.34896 14.3609 14.34701 14.41298
[2,] 16.67733 16.23535 16.33327 16.3991 16.45104 16.4391 16.47299 16.52702
         [,9]    [,10]    [,11]    [,12]
[1,] 14.50687 14.47208 14.50611 14.55715
[2,] 16.57313 16.62792 16.67389 16.74285

$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(13.81, 14.335, 15.44, 16.55, 17.62, 13.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5sh7b1289671289.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.020 0.010 0.000 -0.010 0.000 0.000 0.00 0.020 -0.040 -0.010  0.00 -0.030
[2,] 0.040 0.035 0.015  0.005 0.015 0.010 0.03 0.025  0.015  0.005  0.01  0.020
[3,] 0.090 0.050 0.030  0.030 0.030 0.030 0.05 0.070  0.040  0.040  0.02  0.040
[4,] 0.115 0.065 0.050  0.055 0.045 0.045 0.06 0.075  0.055  0.055  0.05  0.075
[5,] 0.160 0.080 0.050  0.110 0.070 0.080 0.06 0.090  0.060  0.100  0.09  0.100

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

$conf
           [,1]       [,2]        [,3]         [,4]       [,5]        [,6]
[1,] 0.04521121 0.03208448 0.009098565 0.0001408066 0.01208448 0.009098565
[2,] 0.13478879 0.06791552 0.050901435 0.0598591934 0.04791552 0.050901435
           [,7]       [,8]       [,9]      [,10]        [,11]       [,12]
[1,] 0.03208448 0.04014081 0.01611265 0.01014081 -0.003887355 0.007154887
[2,] 0.06791552 0.09985919 0.06388735 0.06985919  0.043887355 0.072845113

$out
[1] 0.12 0.12 0.15 0.24

$group
[1]  2  3  9 12

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

Warning message:
In bxp(list(stats = c(0.0199999999999996, 0.0400000000000009, 0.0899999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6sh7b1289671289.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]
[1,] 13.810 14.140 14.530 15.130 15.750 16.150 16.950 17.62
[2,] 13.925 14.235 14.615 15.345 15.910 16.360 17.160 17.62
[3,] 14.010 14.325 14.665 15.405 15.970 16.540 17.365 17.62
[4,] 14.105 14.450 14.810 15.545 16.065 16.735 17.475 17.62
[5,] 14.130 14.480 14.890 15.650 16.110 16.920 17.650 17.62

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]  [,8]
[1,] 13.9279 14.22694 14.57606 15.31378 15.89930 16.36896 17.22133 17.62
[2,] 14.0921 14.42306 14.75394 15.49622 16.04070 16.71104 17.50867 17.62

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(13.81, 13.925, 14.01, 14.105, 14.13, 14.14,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7286e1289671289.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,] 15.29143 15.250 15.19500
[2,] 15.40357 15.380 15.32000
[3,] 15.49714 15.425 15.39312
[4,] 15.60071 15.545 15.51875
[5,] 15.69000 15.650 15.60000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 15.40722 15.34974 15.30247
[2,] 15.58706 15.50026 15.48378

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15.2914285714286, 15.4035714285714, 15.4971428571429,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1zqq81289671289.ps tmp/1zqq81289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/2zqq81289671289.ps tmp/2zqq81289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/3zqq81289671289.ps tmp/3zqq81289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/4sh7b1289671289.ps tmp/4sh7b1289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/5sh7b1289671289.ps tmp/5sh7b1289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sh7b1289671289.ps tmp/6sh7b1289671289.png",intern=TRUE))
character(0)
> try(system("convert tmp/7286e1289671289.ps tmp/7286e1289671289.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.917   1.459   2.178