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(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)
> 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] 84
> (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] 7 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   NA
 [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     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/freestat/rcomp/tmp/1c9q11289660478.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/2c9q11289660478.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/3c9q11289660478.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/4nj7l1289660478.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,] 68184.0 34077.0 44720.0 36432.0 30936.0 44953.0 36222 27733.0 39401.0
[2,] 68540.5 36776.5 45064.5 37480.5 37860.5 46545.5 37035 32486.5 39920.0
[3,] 68906.0 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] 7 7 7 7 7 7 7 7 7 7 7 7

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 68224.02 37006.99 44742.63 37272.12 37564.39 46470.43 37077.04 30245.89
[2,] 69587.98 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(68184, 68540.5, 68906, 69682.5, 70938, 34077,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5nj7l1289660478.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 56459
[2,]  -7088 -23528 56459
[3,]  -4812 -22505 58871
[4,]  -3578 -21322 59413
[5,]  -1923 -19268 60286

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

$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 56965.57
[2,] 8093.614   901.5861 -2715.885 -21187.61 60776.43

$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/freestat/rcomp/tmp/6ysp61289660478.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   NA
[2,] 37171.5 36639.5 33542.5 35511.5 38465.0 39413.5 28951   NA
[3,] 40062.5 39318.0 37430.5 38684.0 41208.5 41692.5 33447   NA
[4,] 45017.5 44378.5 44182.5 44659.5 46604.5 47679.5 35273   NA
[5,] 50380.0 49267.0 45409.0 50669.0 52073.0 50530.0 35699   NA

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

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

$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/freestat/rcomp/tmp/7ysp61289660478.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] 67566.571  9816.714 68906.000  8771.000 69111.500 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/1c9q11289660478.ps tmp/1c9q11289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/2c9q11289660478.ps tmp/2c9q11289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/3c9q11289660478.ps tmp/3c9q11289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/4nj7l1289660478.ps tmp/4nj7l1289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/5nj7l1289660478.ps tmp/5nj7l1289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ysp61289660478.ps tmp/6ysp61289660478.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ysp61289660478.ps tmp/7ysp61289660478.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.885   1.425   2.135