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(68897,38683,44720,39525,45315,50380,40600,36279,42438,38064,31879,11379,70249,39253,47060,41697,38708,49267,39018,32228,40870,39383,34571,12066,70938,34077,45409,40809,37013,44953,37848,32745,43412,34931,33008,8620,68906,39556,50669,36432,40891,48428,36222,33425,39401,37967,34801,12657,69116,41519,51321,38529,41547,52073,38401,40898,40439,41888,37898,8771,68184,50530,47221,41756,45633,48138,39486,39341,41117,41629,29722,7054,56676,34870,35117,30169,30936,35699,33228,27733,33666,35429,27438,8170,62557)
> 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,] 68897 70249 70938 68906 69116 68184 56676 62557
 [2,] 38683 39253 34077 39556 41519 50530 34870    NA
 [3,] 44720 47060 45409 50669 51321 47221 35117    NA
 [4,] 39525 41697 40809 36432 38529 41756 30169    NA
 [5,] 45315 38708 37013 40891 41547 45633 30936    NA
 [6,] 50380 49267 44953 48428 52073 48138 35699    NA
 [7,] 40600 39018 37848 36222 38401 39486 33228    NA
 [8,] 36279 32228 32745 33425 40898 39341 27733    NA
 [9,] 42438 40870 43412 39401 40439 41117 33666    NA
[10,] 38064 39383 34931 37967 41888 41629 35429    NA
[11,] 31879 34571 33008 34801 37898 29722 27438    NA
[12,] 11379 12066  8620 12657  8771  7054  8170    NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7] [,8]
 [1,] -30214 -30996 -36861 -29350 -27597 -17654 -21806   NA
 [2,]   6037   7807  11332  11113   9802  -3309    247   NA
 [3,]  -5195  -5363  -4600 -14237 -12792  -5465  -4948   NA
 [4,]   5790  -2989  -3796   4459   3018   3877    767   NA
 [5,]   5065  10559   7940   7537  10526   2505   4763   NA
 [6,]  -9780 -10249  -7105 -12206 -13672  -8652  -2471   NA
 [7,]  -4321  -6790  -5103  -2797   2497   -145  -5495   NA
 [8,]   6159   8642  10667   5976   -459   1776   5933   NA
 [9,]  -4374  -1487  -8481  -1434   1449    512   1763   NA
[10,]  -6185  -4812  -1923  -3166  -3990 -11907  -7991   NA
[11,] -20500 -22505 -24388 -22144 -29127 -22668 -19268   NA
[12,]  58870  58872  60286  56459  59413  49622  54387   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/1npph1272476329.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/2npph1272476329.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/3npph1272476329.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/4ggok1272476329.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]
[1,] 62557.0 34077.0 44720.0 36432.0 30936.0 44953.0 36222 27733.0 39401.0
[2,] 65370.5 36776.5 45064.5 37480.5 37860.5 46545.5 37035 32486.5 39920.0
[3,] 68901.5 39253.0 47060.0 39525.0 40891.0 48428.0 38401 33425.0 40870.0
[4,] 69682.5 40537.5 48945.0 41253.0 43431.0 49823.5 39252 37810.0 41777.5
[5,] 70938.0 41519.0 51321.0 41756.0 45633.0 52073.0 40600 40898.0 43412.0
     [,10]   [,11]   [,12]
[1,] 34931 27438.0  7054.0
[2,] 36698 30800.5  8395.0
[3,] 38064 33008.0  8771.0
[4,] 40506 34686.0 11722.5
[5,] 41888 37898.0 12657.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 66492.75 37006.99 44742.63 37272.12 37564.39 46470.43 37077.04 30245.89
[2,] 71310.25 41499.01 49377.37 41777.88 44217.61 50385.57 39724.96 36604.11
         [,9]    [,10]    [,11]     [,12]
[1,] 39760.73 35789.92 30687.64  6783.871
[2,] 41979.27 40338.08 35328.36 10758.129

$out
[1] 56676 50530 35117 30169 35699 33228 33666

$group
[1] 1 2 3 4 6 7 9

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

Warning message:
In bxp(list(stats = c(62557, 65370.5, 68901.5, 69682.5, 70938, 34077,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5ggok1272476329.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,] -36861.0 -3309.0 -14237.0 -3796  2505 -13672.0 -6790  -459.0 -8481.0
[2,] -30605.0  3142.0  -9128.5 -1111  4914 -11227.5 -5299  3854.5 -2930.5
[3,] -29350.0  7807.0  -5363.0  3018  7537  -9780.0 -4321  5976.0 -1434.0
[4,] -24701.5 10457.5  -5071.5  4168  9233  -7878.5 -1471  7400.5   980.5
[5,] -17654.0 11332.0  -4600.0  5790 10559  -7105.0  2497 10667.0  1763.0
      [,10]  [,11]   [,12]
[1,] -11907 -24388 54387.0
[2,]  -7088 -23528 55423.0
[3,]  -4812 -22505 58870.0
[4,]  -3578 -21322 59142.5
[5,]  -1923 -19268 60286.0

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]       [,6]     [,7]
[1,] -32875.47  3438.301 -7785.775 -134.5336  4957.763 -11779.969 -6607.02
[2,] -25824.53 12175.699 -2940.225 6170.5336 10116.237  -7780.031 -2034.98
         [,8]       [,9]     [,10]     [,11]    [,12]
[1,] 3858.386 -3769.5861 -6908.115 -23822.39 56648.77
[2,] 8093.614   901.5861 -2715.885 -21187.61 61091.23

$out
[1]  -2471 -29127  49622

$group
[1]  6 11 12

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

Warning message:
In bxp(list(stats = c(-36861, -30605, -29350, -24701.5, -17654,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6q7nn1272476329.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,] 31879.0 32228.0 32745.0 33425.0 37898.0 29722.0 27438 62557
[2,] 37171.5 36639.5 33542.5 35511.5 38465.0 39413.5 28951 62557
[3,] 40062.5 39318.0 37430.5 38684.0 41208.5 41692.5 33447 62557
[4,] 45017.5 44378.5 44182.5 44659.5 46604.5 47679.5 35273 62557
[5,] 50380.0 49267.0 45409.0 50669.0 52073.0 50530.0 35699 62557

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]  [,8]
[1,] 36483.89 35788.19 32577.52 34511.54 37496.02 37922.32 30563.49 62557
[2,] 43641.11 42847.81 42283.48 42856.46 44920.98 45462.68 36330.51 62557

$out
 [1] 68897 11379 70249 12066 70938  8620 68906 12657 69116  8771 68184  7054
[13] 56676  8170

$group
 [1] 1 1 2 2 3 3 4 4 5 5 6 6 7 7

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

Warning message:
In bxp(list(stats = c(31879, 37171.5, 40062.5, 45017.5, 50380, 32228,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7q7nn1272476329.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,] 32759.57 33008.0 32743.25
[2,] 36246.57 35744.5 36645.88
[3,] 39127.07 39389.0 39011.88
[4,] 43061.43 43975.5 43926.75
[5,] 46991.14 48428.0 48184.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 36018.77 35634.79 35691.02
[2,] 42235.37 43143.21 42332.73

$out
[1] 66940.375  9816.714 68901.500  8771.000 68088.250 10058.750

$group
[1] 1 1 2 2 3 3

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

Warning message:
In bxp(list(stats = c(32759.5714285714, 36246.5714285714, 39127.0714285714,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1npph1272476329.ps tmp/1npph1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/2npph1272476329.ps tmp/2npph1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/3npph1272476329.ps tmp/3npph1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/4ggok1272476329.ps tmp/4ggok1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/5ggok1272476329.ps tmp/5ggok1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/6q7nn1272476329.ps tmp/6q7nn1272476329.png",intern=TRUE))
character(0)
> try(system("convert tmp/7q7nn1272476329.ps tmp/7q7nn1272476329.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.173   0.932   2.357