R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(31.514,27.071,29.462,26.105,22.397,23.843,21.705,18.089,20.764,25.316,17.704,15.548,28.029,29.383,36.438,32.034,22.679,24.319,18.004,17.537,20.366,22.782,19.169,13.807,29.743,25.591,29.096,26.482,22.405,27.044,17.970,18.730,19.684,19.785,18.479,10.698,31.956,29.506,34.506,27.165,26.736,23.691,18.157,17.328,18.205,20.995,17.382,9.367,31.124,26.551,30.651,25.859,25.100,25.778,20.418,18.688,20.424,24.776,19.814,12.738,31.566,30.111,30.019,31.934,25.826,26.835,20.205,17.789,20.520,22.518,15.572,11.509,25.447,24.090,27.786,26.195,20.516,22.759,19.028,16.971,20.036,22.485,18.730,14.538,27.561,25.985,34.670,32.066,27.186,29.586,21.359,21.553,19.573,24.256,22.380,16.167,27.297,28.287)
> 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] 98
> (np <- floor(n / par1))
[1] 8
> 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] 9 9 8 8 8 8 8 8 8 8 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,] 31.514 28.029 29.743 31.956 31.124 31.566 25.447 27.561 27.297
 [2,] 27.071 29.383 25.591 29.506 26.551 30.111 24.090 25.985 28.287
 [3,] 29.462 36.438 29.096 34.506 30.651 30.019 27.786 34.670     NA
 [4,] 26.105 32.034 26.482 27.165 25.859 31.934 26.195 32.066     NA
 [5,] 22.397 22.679 22.405 26.736 25.100 25.826 20.516 27.186     NA
 [6,] 23.843 24.319 27.044 23.691 25.778 26.835 22.759 29.586     NA
 [7,] 21.705 18.004 17.970 18.157 20.418 20.205 19.028 21.359     NA
 [8,] 18.089 17.537 18.730 17.328 18.688 17.789 16.971 21.553     NA
 [9,] 20.764 20.366 19.684 18.205 20.424 20.520 20.036 19.573     NA
[10,] 25.316 22.782 19.785 20.995 24.776 22.518 22.485 24.256     NA
[11,] 17.704 19.169 18.479 17.382 19.814 15.572 18.730 22.380     NA
[12,] 15.548 13.807 10.698  9.367 12.738 11.509 14.538 16.167     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8] [,9]
 [1,] -4.443  1.354 -4.152 -2.450 -4.573 -1.455 -1.357 -1.576 0.99
 [2,]  2.391  7.055  3.505  5.000  4.100 -0.092  3.696  8.685   NA
 [3,] -3.357 -4.404 -2.614 -7.341 -4.792  1.915 -1.591 -2.604   NA
 [4,] -3.708 -9.355 -4.077 -0.429 -0.759 -6.108 -5.679 -4.880   NA
 [5,]  1.446  1.640  4.639 -3.045  0.678  1.009  2.243  2.400   NA
 [6,] -2.138 -6.315 -9.074 -5.534 -5.360 -6.630 -3.731 -8.227   NA
 [7,] -3.616 -0.467  0.760 -0.829 -1.730 -2.416 -2.057  0.194   NA
 [8,]  2.675  2.829  0.954  0.877  1.736  2.731  3.065 -1.980   NA
 [9,]  4.552  2.416  0.101  2.790  4.352  1.998  2.449  4.683   NA
[10,] -7.612 -3.613 -1.306 -3.613 -4.962 -6.946 -3.755 -1.876   NA
[11,] -2.156 -5.362 -7.781 -8.015 -7.076 -4.063 -4.192 -6.213   NA
[12,] 12.481 15.936 21.258 21.757 18.828 13.938 13.023 11.130   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/16vqj1305836454.ps",horizontal=F,onefile=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/2wuby1305836455.ps",horizontal=F,onefile=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/3j2va1305836455.ps",horizontal=F,onefile=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/41ztm1305836455.ps",horizontal=F,onefile=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]
[1,] 25.447 24.090 27.786 25.8590 20.5160 22.7590 17.9700 16.9710 19.5730
[2,] 27.561 25.985 29.279 26.1500 22.4010 23.7670 18.0805 17.4325 19.6285
[3,] 29.743 27.071 30.335 26.8235 23.8895 25.0485 19.6165 17.9390 20.2010
[4,] 31.514 29.383 34.588 31.9840 26.2810 26.9395 20.8885 18.7090 20.4720
[5,] 31.956 30.111 36.438 32.0660 27.1860 29.5860 21.7050 18.7300 20.7640
      [,10]   [,11]   [,12]
[1,] 19.785 15.5720  9.3670
[2,] 21.740 17.5430 11.1035
[3,] 22.650 18.6045 13.2725
[4,] 24.516 19.4915 15.0430
[5,] 25.316 22.3800 16.1670

$n
 [1] 9 9 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 27.66109 25.28139 27.36932 23.56454 21.72208 23.27630 18.04791 17.22593
[2,] 31.82491 28.86061 33.30068 30.08246 26.05692 26.82070 21.18509 18.65207
         [,9]    [,10]    [,11]    [,12]
[1,] 19.72981 21.09929 17.51604 11.07184
[2,] 20.67219 24.20071 19.69296 15.47316

$out
[1] 21.553 18.205

$group
[1] 8 9

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

Warning message:
In bxp(list(stats = c(25.447, 27.561, 29.743, 31.514, 31.956, 24.09,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5v1d51305836455.ps",horizontal=F,onefile=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]
[1,] -4.573 -0.0920 -7.3410 -9.3550 0.6780 -9.0740 -3.6160 0.8770 0.1010
[2,] -4.152  2.9480 -4.5980 -5.8935 0.8435 -7.4285 -2.2365 0.9155 2.2070
[3,] -1.576  3.8980 -2.9855 -4.4785 1.5430 -5.9245 -1.2795 2.2055 2.6195
[4,] -1.357  6.0275 -2.0975 -2.2335 2.3215 -4.5455 -0.1365 2.7800 4.4520
[5,]  1.354  8.6850 -1.5910 -0.4290 2.4000 -2.1380  0.7600 3.0650 4.6830
       [,10]   [,11]  [,12]
[1,] -7.6120 -8.0150 11.130
[2,] -5.9540 -7.4285 12.752
[3,] -3.6840 -5.7875 14.937
[4,] -2.7445 -4.1275 20.043
[5,] -1.3060 -2.1560 21.757

$n
 [1] 9 8 8 8 8 8 8 8 8 8 8 8

$conf
           [,1]     [,2]      [,3]      [,4]     [,5]      [,6]       [,7]
[1,] -3.0480333 2.177747 -4.382315 -6.523029 0.717368 -7.534985 -2.4525901
[2,] -0.1039667 5.618253 -1.588685 -2.433971 2.368632 -4.314015 -0.1064099
         [,8]     [,9]     [,10]     [,11]    [,12]
[1,] 1.163964 1.365411 -5.476873 -7.631486 10.86414
[2,] 3.247036 3.873589 -1.891127 -3.943514 19.00986

$out
[1]  1.915  4.639 -3.045 -1.980

$group
[1] 3 5 5 8

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

Warning message:
In bxp(list(stats = c(-4.573, -4.152, -1.576, -1.357, 1.354, -0.0920000000000023,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6ixz71305836455.ps",horizontal=F,onefile=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]
[1,] 15.5480 13.8070 10.6980  9.3670 12.738 11.509 14.5380 16.1670 27.297
[2,] 19.4265 18.5865 18.6045 17.7695 20.116 18.997 18.8790 21.4560 27.297
[3,] 23.1200 22.7305 21.0950 22.3430 24.938 24.172 21.5005 25.1205 27.792
[4,] 26.5880 28.7060 26.7630 28.3355 26.205 30.065 24.7685 28.5735 28.287
[5,] 31.5140 36.4380 29.7430 34.5060 31.124 31.934 27.7860 34.6700 28.287

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 19.85359 18.11493 17.37385 17.52378 22.16077 19.12381 18.81426 21.87416
[2,] 26.38641 27.34607 24.81615 27.16222 27.71523 29.22019 24.18674 28.36684
         [,9]
[1,] 26.68594
[2,] 28.89806

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15.548, 19.4265, 23.12, 26.588, 31.514, 13.807,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/74jws1305836455.ps",horizontal=F,onefile=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,] 13.04650 13.27250 13.04838
[2,] 19.12975 19.11050 18.93144
[3,] 23.48488 23.26975 23.73875
[4,] 27.93861 26.94725 28.37487
[5,] 31.57850 30.33500 31.95875

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 19.46709 19.69536 19.43153
[2,] 27.50266 26.84414 28.04597

$out
numeric(0)

$group
numeric(0)

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

> dev.off()
null device 
          1 
> 
> try(system("convert tmp/16vqj1305836454.ps tmp/16vqj1305836454.png",intern=TRUE))
character(0)
> try(system("convert tmp/2wuby1305836455.ps tmp/2wuby1305836455.png",intern=TRUE))
character(0)
> try(system("convert tmp/3j2va1305836455.ps tmp/3j2va1305836455.png",intern=TRUE))
character(0)
> try(system("convert tmp/41ztm1305836455.ps tmp/41ztm1305836455.png",intern=TRUE))
character(0)
> try(system("convert tmp/5v1d51305836455.ps tmp/5v1d51305836455.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ixz71305836455.ps tmp/6ixz71305836455.png",intern=TRUE))
character(0)
> try(system("convert tmp/74jws1305836455.ps tmp/74jws1305836455.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.728   0.508   3.222