R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(12.008,9.169,8.788,8.417,8.247,8.197,8.236,8.253,7.733,8.366,8.626,8.863,10.102,8.463,9.114,8.563,8.872,8.301,8.301,8.278,7.736,7.973,8.268,9.476,11.100,8.962,9.173,8.738,8.459,8.078,8.411,8.291,7.810,8.616,8.312,9.692,9.911,8.915,9.452,9.112,8.472,8.230,8.384,8.625,8.221,8.649,8.625,10.443,10.357,8.586,8.892,8.329,8.101,7.922,8.120,7.838,7.735,8.406,8.209,9.451,10.041,9.411,10.405,8.467,8.464,8.102,7.627,7.513,7.510,8.291,8.064,9.383,9.706,8.579,9.474,8.318,8.213,8.059,9.111,7.708,7.680,8.014,8.007,8.718,9.486,9.113,9.025,8.476,7.952,7.759,7.835,7.600,7.651,8.319,8.812,8.630)
> 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] 96
> (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] 8 8 8 8 8 8 8 8 8 8 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]  [,7]  [,8] [,9]
 [1,] 12.008 10.102 11.100  9.911 10.357 10.041 9.706 9.486   NA
 [2,]  9.169  8.463  8.962  8.915  8.586  9.411 8.579 9.113   NA
 [3,]  8.788  9.114  9.173  9.452  8.892 10.405 9.474 9.025   NA
 [4,]  8.417  8.563  8.738  9.112  8.329  8.467 8.318 8.476   NA
 [5,]  8.247  8.872  8.459  8.472  8.101  8.464 8.213 7.952   NA
 [6,]  8.197  8.301  8.078  8.230  7.922  8.102 8.059 7.759   NA
 [7,]  8.236  8.301  8.411  8.384  8.120  7.627 9.111 7.835   NA
 [8,]  8.253  8.278  8.291  8.625  7.838  7.513 7.708 7.600   NA
 [9,]  7.733  7.736  7.810  8.221  7.735  7.510 7.680 7.651   NA
[10,]  8.366  7.973  8.616  8.649  8.406  8.291 8.014 8.319   NA
[11,]  8.626  8.268  8.312  8.625  8.209  8.064 8.007 8.812   NA
[12,]  8.863  9.476  9.692 10.443  9.451  9.383 8.718 8.630   NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8] [,9]
 [1,] -2.839 -1.639 -2.138 -0.996 -1.771 -0.630 -1.127 -0.373   NA
 [2,] -0.381  0.651  0.211  0.537  0.306  0.994  0.895 -0.088   NA
 [3,] -0.371 -0.551 -0.435 -0.340 -0.563 -1.938 -1.156 -0.549   NA
 [4,] -0.170  0.309 -0.279 -0.640 -0.228 -0.003 -0.105 -0.524   NA
 [5,] -0.050 -0.571 -0.381 -0.242 -0.179 -0.362 -0.154 -0.193   NA
 [6,]  0.039  0.000  0.333  0.154  0.198 -0.475  1.052  0.076   NA
 [7,]  0.017 -0.023 -0.120  0.241 -0.282 -0.114 -1.403 -0.235   NA
 [8,] -0.520 -0.542 -0.481 -0.404 -0.103 -0.003 -0.028  0.051   NA
 [9,]  0.633  0.237  0.806  0.428  0.671  0.781  0.334  0.668   NA
[10,]  0.260  0.295 -0.304 -0.024 -0.197 -0.227 -0.007  0.493   NA
[11,]  0.237  1.208  1.380  1.818  1.242  1.319  0.711 -0.182   NA
[12,]  1.239  1.624  0.219 -0.086  0.590  0.323  0.768     NA   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/1io2t1290625209.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/2io2t1290625209.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/3io2t1290625209.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/4bf2e1290625209.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]
[1,]  9.4860 8.4630 8.7880 8.3180 7.952 7.7590 7.6270 7.5130 7.5100 7.9730
[2,]  9.8085 8.5825 8.9585 8.3730 8.157 7.9905 7.9775 7.6540 7.6655 8.1525
[3,] 10.0715 8.9385 9.1435 8.4715 8.353 8.0900 8.2685 8.0455 7.7340 8.3425
[4,] 10.7285 9.1410 9.4630 8.6505 8.468 8.2135 8.3975 8.2845 7.7730 8.5110
[5,] 12.0080 9.4110 9.4740 8.7380 8.872 8.3010 8.4110 8.6250 7.8100 8.6490
      [,11]   [,12]
[1,] 8.0070  8.6300
[2,] 8.1365  8.7905
[3,] 8.2900  9.4170
[4,] 8.6255  9.5840
[5,] 8.8120 10.4430

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

$conf
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,]  9.557575 8.626514 8.861679 8.316485 8.179271 7.965429 8.033882 7.693294
[2,] 10.585425 9.250486 9.425321 8.626515 8.526729 8.214571 8.503118 8.397706
         [,9]    [,10]    [,11]   [,12]
[1,] 7.673949 8.142237 8.016838 8.97374
[2,] 7.794051 8.542763 8.563162 9.86026

$out
[1] 10.405  9.112  9.111  8.221

$group
[1] 3 4 7 9

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

Warning message:
In bxp(list(stats = c(9.486, 9.8085, 10.0715, 10.7285, 12.008, 8.463,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5bf2e1290625209.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]
[1,] -2.8390 -0.3810 -1.1560 -0.6400 -0.5710 0.0000 -0.2820 -0.5420 0.2370
[2,] -1.9545  0.0615 -0.8595 -0.4015 -0.3715 0.0195 -0.2585 -0.5005 0.3810
[3,] -1.3830  0.4215 -0.5500 -0.1990 -0.2175 0.1150 -0.1170 -0.2535 0.6505
[4,] -0.8130  0.7730 -0.4030 -0.0540 -0.1665 0.2655 -0.0030 -0.0155 0.7260
[5,] -0.3730  0.9940 -0.3400  0.3090 -0.0500 0.3330  0.2410  0.0510 0.8060
       [,10]   [,11]   [,12]
[1,] -0.3040 -0.1820 -0.0860
[2,] -0.2120  0.4740  0.2710
[3,] -0.0155  1.2250  0.5900
[4,]  0.2775  1.3495  1.0035
[5,]  0.4930  1.8180  1.6240

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

$conf
           [,1]       [,2]       [,3]         [,4]       [,5]        [,6]
[1,] -2.0206583 0.02404588 -0.8050075 -0.393118489 -0.3320159 -0.02241913
[2,] -0.7453417 0.81895412 -0.2949925 -0.004881511 -0.1029841  0.25241913
            [,7]        [,8]      [,9]      [,10]     [,11]     [,12]
[1,] -0.25972597 -0.52442796 0.4577780 -0.2889417 0.7359331 0.1525628
[2,]  0.02572597  0.01742796 0.8432220  0.2579417 1.7140669 1.0274372

$out
[1] -1.938 -0.475  1.052 -1.403

$group
[1] 3 6 6 7

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

Warning message:
In bxp(list(stats = c(-2.839, -1.9545, -1.383, -0.813, -0.373000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6bf2e1290625209.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] [,9]
[1,] 7.7330 7.736 7.8100  8.221 7.7350  7.5100 7.6800 7.6000   NA
[2,] 8.2415 8.273 8.3015  8.428 8.0115  7.8455 8.0105 7.7970   NA
[3,] 8.3915 8.382 8.5375  8.637 8.2690  8.3775 8.2655 8.3975   NA
[4,] 8.8255 8.993 9.0675  9.282 8.7390  9.3970 8.9145 8.9185   NA
[5,] 9.1690 9.476 9.6920 10.443 9.4510 10.4050 9.7060 9.4860   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]    [,7]     [,8] [,9]
[1,] 8.125134 8.053603 8.188122 8.247485 7.937182 7.66985 7.85318 7.885976   NA
[2,] 8.657866 8.710397 8.886878 9.026515 8.600818 9.08515 8.67782 8.909024   NA

$out
[1] 12.008 10.102 11.100 10.357

$group
[1] 1 2 3 5

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

Warning message:
In bxp(list(stats = c(7.733, 8.2415, 8.3915, 8.8255, 9.169, 7.736,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/73pjh1290625209.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,]  7.759500  7.73400  7.713625
[2,]  8.167063  8.17925  8.167375
[3,]  8.356438  8.34775  8.369563
[4,]  9.095063  9.04100  9.017000
[5,] 10.338875 10.07150 10.201250

$n
[1] 12 12 12

$conf
         [,1]   [,2]     [,3]
[1,] 7.933170 7.9547 7.982043
[2,] 8.779705 8.7408 8.757082

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(7.7595, 8.1670625, 8.3564375, 9.0950625, 10.338875,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1io2t1290625209.ps tmp/1io2t1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/2io2t1290625209.ps tmp/2io2t1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/3io2t1290625209.ps tmp/3io2t1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/4bf2e1290625209.ps tmp/4bf2e1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/5bf2e1290625209.ps tmp/5bf2e1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bf2e1290625209.ps tmp/6bf2e1290625209.png",intern=TRUE))
character(0)
> try(system("convert tmp/73pjh1290625209.ps tmp/73pjh1290625209.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.235   0.971  10.492