R version 2.11.1 (2010-05-31)
Copyright (C) 2010 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(17.78,17.75,17.74,17.71,17.72,17.7,17.62,17.65,17.56,17.5,17.45,17.43,17.39,17.34,17.27,17.16,17.16,17.11,16.95,16.92,16.91,16.81,16.66,16.64,16.58,16.5,16.49,16.42,16.3,16.18,16.15,16.11,16.07,16.08,16.06,15.99,15.99,15.95,15.93,15.94,15.89,15.86,15.75,15.65,15.59,15.55,15.54,15.47,15.41,15.4,15.4,15.36,15.33,15.25,15.13,14.89,14.89,14.84,14.78,14.7,14.68,14.65,14.61,14.62,14.62,14.58,14.53,14.48,14.46,14.44,14.48,14.39,14.33,14.32,14.29)
> 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] 75
> (np <- floor(n / par1))
[1] 6
> 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] 7 7 7 6 6 6 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 17.78 17.39 16.58 15.99 15.41 14.68 14.33
 [2,] 17.75 17.34 16.50 15.95 15.40 14.65 14.32
 [3,] 17.74 17.27 16.49 15.93 15.40 14.61 14.29
 [4,] 17.71 17.16 16.42 15.94 15.36 14.62    NA
 [5,] 17.72 17.16 16.30 15.89 15.33 14.62    NA
 [6,] 17.70 17.11 16.18 15.86 15.25 14.58    NA
 [7,] 17.62 16.95 16.15 15.75 15.13 14.53    NA
 [8,] 17.65 16.92 16.11 15.65 14.89 14.48    NA
 [9,] 17.56 16.91 16.07 15.59 14.89 14.46    NA
[10,] 17.50 16.81 16.08 15.55 14.84 14.44    NA
[11,] 17.45 16.66 16.06 15.54 14.78 14.48    NA
[12,] 17.43 16.64 15.99 15.47 14.70 14.39    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] -0.03 -0.05 -0.08 -0.04 -0.01 -0.03 -0.01
 [2,] -0.01 -0.07 -0.01 -0.02  0.00 -0.04 -0.03
 [3,] -0.03 -0.11 -0.07  0.01 -0.04  0.01    NA
 [4,]  0.01  0.00 -0.12 -0.05 -0.03  0.00    NA
 [5,] -0.02 -0.05 -0.12 -0.03 -0.08 -0.04    NA
 [6,] -0.08 -0.16 -0.03 -0.11 -0.12 -0.05    NA
 [7,]  0.03 -0.03 -0.04 -0.10 -0.24 -0.05    NA
 [8,] -0.09 -0.01 -0.04 -0.06  0.00 -0.02    NA
 [9,] -0.06 -0.10  0.01 -0.04 -0.05 -0.02    NA
[10,] -0.05 -0.15 -0.02 -0.01 -0.06  0.04    NA
[11,] -0.02 -0.02 -0.07 -0.07 -0.08 -0.09    NA
[12,] -0.04 -0.06  0.00 -0.06 -0.02 -0.06    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/rcomp/tmp/1bf3q1290678963.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/rcomp/tmp/2bf3q1290678963.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/rcomp/tmp/3lokt1290678963.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/rcomp/tmp/4lokt1290678963.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,] 14.330 14.320 14.290 14.62 14.620 14.58 14.53 14.48 14.46 14.440 14.48
[2,] 15.045 15.025 15.005 15.36 15.330 15.25 15.13 14.89 14.89 14.840 14.78
[3,] 15.990 15.950 15.930 16.18 16.095 16.02 15.95 15.88 15.83 15.815 15.80
[4,] 16.985 16.920 16.880 17.16 17.160 17.11 16.95 16.92 16.91 16.810 16.66
[5,] 17.780 17.750 17.740 17.71 17.720 17.70 17.62 17.65 17.56 17.500 17.45
     [,12]
[1,] 14.39
[2,] 14.70
[3,] 15.73
[4,] 16.64
[5,] 17.43

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 14.83146 14.81834 14.81028 15.01894 14.91459 14.82024 14.77604 14.57058
[2,] 17.14854 17.08166 17.04972 17.34106 17.27541 17.21976 17.12396 17.18942
         [,9]    [,10]    [,11]    [,12]
[1,] 14.52703 14.54429 14.58734 14.47864
[2,] 17.13297 17.08571 17.01266 16.98136

$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(14.33, 15.045, 15.99, 16.985, 17.78, 14.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5lokt1290678963.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]
[1,] -0.080 -0.070 -0.110 -0.120 -0.120 -0.160 -0.100 -0.09 -0.100 -0.060 -0.09
[2,] -0.045 -0.035 -0.070 -0.050 -0.080 -0.120 -0.100 -0.06 -0.060 -0.060 -0.08
[3,] -0.030 -0.020 -0.035 -0.015 -0.045 -0.095 -0.045 -0.03 -0.045 -0.035 -0.07
[4,] -0.020 -0.010  0.010  0.000 -0.030 -0.050 -0.030 -0.01 -0.020 -0.010 -0.02
[5,] -0.010  0.000  0.010  0.010 -0.020 -0.030  0.030  0.00  0.010  0.040 -0.02
     [,12]
[1,] -0.06
[2,] -0.06
[3,] -0.05
[4,] -0.02
[5,]  0.00

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

$conf
            [,1]         [,2]        [,3]        [,4]        [,5]        [,6]
[1,] -0.04492960 -0.034929597 -0.08660258 -0.04725161 -0.07725161 -0.14015226
[2,] -0.01507040 -0.005070403  0.01660258  0.01725161 -0.01274839 -0.04984774
              [,7]         [,8]        [,9]        [,10]       [,11]
[1,] -0.0901522609 -0.062251615 -0.07080129 -0.067251615 -0.10870194
[2,]  0.0001522609  0.002251615 -0.01919871 -0.002748385 -0.03129806
           [,12]
[1,] -0.07580129
[2,] -0.02419871

$out
[1] -0.24 -0.15

$group
[1]  7 10

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

Warning message:
In bxp(list(stats = c(-0.0799999999999983, -0.0450000000000008,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6wfje1290678963.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]
[1,] 17.430 16.640 15.990 15.470 14.700 14.390 14.290
[2,] 17.530 16.860 16.075 15.570 14.865 14.470 14.305
[3,] 17.675 17.030 16.165 15.805 15.190 14.555 14.320
[4,] 17.730 17.215 16.455 15.935 15.380 14.620 14.325
[5,] 17.780 17.390 16.580 15.990 15.410 14.680 14.330

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 17.58378 16.86808 15.99168 15.63852 14.95511 14.48658 14.30176
[2,] 17.76622 17.19192 16.33832 15.97148 15.42489 14.62342 14.33824

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(17.43, 17.53, 17.675, 17.73, 17.78, 16.64,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7wfje1290678963.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.77000 15.7300 15.68500
[2,] 15.89167 15.8225 15.85250
[3,] 15.97429 15.9400 15.95750
[4,] 16.06810 16.0050 16.07875
[5,] 16.20167 16.1800 16.24000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 15.89382 15.85676 15.85431
[2,] 16.05476 16.02324 16.06069

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15.77, 15.8916666666667, 15.9742857142857,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1bf3q1290678963.ps tmp/1bf3q1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/2bf3q1290678963.ps tmp/2bf3q1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/3lokt1290678963.ps tmp/3lokt1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/4lokt1290678963.ps tmp/4lokt1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lokt1290678963.ps tmp/5lokt1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/6wfje1290678963.ps tmp/6wfje1290678963.png",intern=TRUE))
character(0)
> try(system("convert tmp/7wfje1290678963.ps tmp/7wfje1290678963.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.940   1.030   2.898