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.

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(60804,57907,54355,52536,49081,48877,64599,75314,71209,65210,59829,57656,57428,55315,52790,51050,48519,48354,65333,73990,72755,67424,59214,57427,56681,55437,53600,51641,49478,50124,71313,76208,74387,69520,64735,63413,62553,60109,57764,55667,53103,55301,76795,80928,79213,72759,67802,66940,66396,67539,67776,68014,68251,68488,68725,68962,69200,69437,68212,65444,63181,61198,59010,56388,53723,55340,75352,79817,78289,71892,66448,64167,61250,59580,56417,54662,53349,55385,73546,77683,74995,67282,60742,57283,57314,54704,51578,49962,46252,47234,64708,68753,62970,57474,52494,51831,51663,49637,46679,45557,41630,44417,60070,63157)
> 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] 104
> (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 9 9 9 9 9 9 8 8 8 8
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]
 [1,] 60804 57428 56681 62553 66396 63181 61250 57314 51663
 [2,] 57907 55315 55437 60109 67539 61198 59580 54704 49637
 [3,] 54355 52790 53600 57764 67776 59010 56417 51578 46679
 [4,] 52536 51050 51641 55667 68014 56388 54662 49962 45557
 [5,] 49081 48519 49478 53103 68251 53723 53349 46252 41630
 [6,] 48877 48354 50124 55301 68488 55340 55385 47234 44417
 [7,] 64599 65333 71313 76795 68725 75352 73546 64708 60070
 [8,] 75314 73990 76208 80928 68962 79817 77683 68753 63157
 [9,] 71209 72755 74387 79213 69200 78289 74995 62970    NA
[10,] 65210 67424 69520 72759 69437 71892 67282 57474    NA
[11,] 59829 59214 64735 67802 68212 66448 60742 52494    NA
[12,] 57656 57427 63413 66940 65444 64167 57283 51831    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]
 [1,] -2897 -2113 -1244 -2444  1143 -1983 -1670 -2610 -2026
 [2,] -3552 -2525 -1837 -2345   237 -2188 -3163 -3126 -2958
 [3,] -1819 -1740 -1959 -2097   238 -2622 -1755 -1616 -1122
 [4,] -3455 -2531 -2163 -2564   237 -2665 -1313 -3710 -3927
 [5,]  -204  -165   646  2198   237  1617  2036   982  2787
 [6,] 15722 16979 21189 21494   237 20012 18161 17474 15653
 [7,] 10715  8657  4895  4133   237  4465  4137  4045  3087
 [8,] -4105 -1235 -1821 -1715   238 -1528 -2688 -5783    NA
 [9,] -5999 -5331 -4867 -6454   237 -6397 -7713 -5496    NA
[10,] -5381 -8210 -4785 -4957 -1225 -5444 -6540 -4980    NA
[11,] -2173 -1787 -1322  -862 -2768 -2281 -3459  -663    NA
[12,]  -228  -746  -860  -544 -2263 -2917    31  -168    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/1bwi21229080287.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/2rlw51229080287.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/3u1ec1229080287.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/4e8do1229080287.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,] 51663 49637 46679 45557 41630 44417 60070 63157 62970.0 65210.0 52494.0
[2,] 57314 55315 52790 51050 48519 48354 64708 68962 70204.5 66246.0 59521.5
[3,] 60804 57907 54355 52536 49478 50124 68725 75314 73571.0 68430.5 62738.5
[4,] 62553 60109 57764 55667 53349 55340 73546 77683 76642.0 70706.0 67125.0
[5,] 66396 61198 59010 56388 53723 55385 76795 80928 79213.0 72759.0 68212.0
       [,12]
[1,] 51831.0
[2,] 57355.0
[3,] 60534.5
[4,] 64805.5
[5,] 66940.0

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 58044.79 55382.16 51735.36 50104.38 46934.2 46444.71 64070.32 70720.94
[2,] 63563.21 60431.84 56974.64 54967.62 52021.8 53803.29 73379.68 79907.06
         [,9]    [,10]    [,11]    [,12]
[1,] 69974.92 65939.08 58491.08 56372.54
[2,] 77167.08 70921.92 66985.92 64696.46

$out
[1] 67539 67776 68014 68251 68488 57474

$group
[1]  2  3  4  5  6 10

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

Warning message:
In bxp(list(stats = c(51663, 57314, 60804, 62553, 66396, 49637,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5a8131229080287.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,] -2897 -3552 -2097 -3927 -204 15653 3087 -5783.0 -7713.0 -6540.0 -3459.0
[2,] -2444 -3126 -1959 -3455  237 15722 4045 -3396.5 -6425.5 -5992.0 -2524.5
[3,] -2026 -2525 -1755 -2564  982 17474 4137 -1768.0 -5747.5 -5180.5 -1980.0
[4,] -1670 -2188 -1616 -2163 2036 20012 4895 -1381.5 -5099.0 -4871.0 -1092.0
[5,] -1244 -1837 -1122 -1313 2787 21494 4895   238.0 -4867.0 -4785.0  -663.0
       [,12]
[1,] -2917.0
[2,] -1561.5
[3,]  -645.0
[4,]  -198.0
[5,]    31.0

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

$conf
         [,1]      [,2]      [,3]      [,4]       [,5]    [,6]     [,7]
[1,] -2433.64 -3019.013 -1935.647 -3244.453   34.52667 15214.6 3689.333
[2,] -1618.36 -2030.987 -1574.353 -1883.547 1929.47333 19733.4 4584.667
           [,8]      [,9]     [,10]     [,11]      [,12]
[1,] -2893.6079 -6488.502 -5806.707 -2780.215 -1406.6707
[2,]  -642.3921 -5006.498 -4554.293 -1179.785   116.6707

$out
 [1]  1143   237   238 -2622   237   237 10715  8657   237   237 -8210 -1225

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

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

Warning message:
In bxp(list(stats = c(-2897, -2444, -2026, -1670, -1244, -3552,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6tr8m1229080287.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,] 48877.0 48354.0 49478.0 53103.0 66396.0 53723 53349 46252 41630.0
[2,] 53445.5 51920.0 52620.5 56715.5 67657.5 57699 55901 50770 44987.0
[3,] 58868.0 57427.5 60047.0 64746.5 68231.5 63674 60161 53599 48158.0
[4,] 64904.5 66378.5 70416.5 74777.0 68843.5 73622 70414 60222 55866.5
[5,] 75314.0 73990.0 76208.0 80928.0 69437.0 79817 77683 68753 63157.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 53641.47 50832.88 51930.12 56508.53 67690.56 56411.41 53541.52 49287.88
[2,] 64094.53 64022.12 68163.88 72984.47 68772.44 70936.59 66780.48 57910.12
         [,9]
[1,] 42080.56
[2,] 54235.44

$out
[1] 65444

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(48877, 53445.5, 58868, 64904.5, 75314, 48354,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7drok1229080287.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,] 51487.33 49478.00 50934.00
[2,] 54747.00 53445.50 54317.75
[3,] 60108.40 60669.25 60436.06
[4,] 68281.32 68577.75 68782.75
[5,] 73868.00 75314.00 73322.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 53935.30 53767.33 53838.48
[2,] 66281.49 67571.17 67033.65

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(51487.3333333333, 54747, 60108.3958333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1bwi21229080287.ps tmp/1bwi21229080287.png")
> system("convert tmp/2rlw51229080287.ps tmp/2rlw51229080287.png")
> system("convert tmp/3u1ec1229080287.ps tmp/3u1ec1229080287.png")
> system("convert tmp/4e8do1229080287.ps tmp/4e8do1229080287.png")
> system("convert tmp/5a8131229080287.ps tmp/5a8131229080287.png")
> system("convert tmp/6tr8m1229080287.ps tmp/6tr8m1229080287.png")
> system("convert tmp/7drok1229080287.ps tmp/7drok1229080287.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.204   0.956   1.754