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(21863,20403,18792,17931,16475,16205,25134,31896,26537,22801,20200,19666,19809,18799,17884,17512,16327,16880,26537,31867,29427,25800,22041,21759,21333,20462,19594,18564,17640,18614,32562,35640,31865,28117,25508,25006,24452,22643,21474,20500,19505,21769,36062,38633,34629,30184,27271,26841,26482,25538,23789,22386,21087,22891,36192,38922,34669,30197,27001,25891,24879,23662,22741,21615,20305,21877,35369,37941,33480,29757,26323,25359,22207,21763,19944,19662,18624,19902,31726,32860,28894,22949,19758,18420,18245,16761,15341,14271,13418,15218,26485,27457,21402,17879,15607,15626,15303,14296,13686,12948,11609,14602,23629,24680)
> 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,] 21863 19809 21333 24452 26482 24879 22207 18245 15303
 [2,] 20403 18799 20462 22643 25538 23662 21763 16761 14296
 [3,] 18792 17884 19594 21474 23789 22741 19944 15341 13686
 [4,] 17931 17512 18564 20500 22386 21615 19662 14271 12948
 [5,] 16475 16327 17640 19505 21087 20305 18624 13418 11609
 [6,] 16205 16880 18614 21769 22891 21877 19902 15218 14602
 [7,] 25134 26537 32562 36062 36192 35369 31726 26485 23629
 [8,] 31896 31867 35640 38633 38922 37941 32860 27457 24680
 [9,] 26537 29427 31865 34629 34669 33480 28894 21402    NA
[10,] 22801 25800 28117 30184 30197 29757 22949 17879    NA
[11,] 20200 22041 25508 27271 27001 26323 19758 15607    NA
[12,] 19666 21759 25006 26841 25891 25359 18420 15626    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9]
 [1,] -1460 -1010  -871 -1809  -944 -1217  -444 -1484 -1007
 [2,] -1611  -915  -868 -1169 -1749  -921 -1819 -1420  -610
 [3,]  -861  -372 -1030  -974 -1403 -1126  -282 -1070  -738
 [4,] -1456 -1185  -924  -995 -1299 -1310 -1038  -853 -1339
 [5,]  -270   553   974  2264  1804  1572  1278  1800  2993
 [6,]  8929  9657 13948 14293 13301 13492 11824 11267  9027
 [7,]  6762  5330  3078  2571  2730  2572  1134   972  1051
 [8,] -5359 -2440 -3775 -4004 -4253 -4461 -3966 -6055    NA
 [9,] -3736 -3627 -3748 -4445 -4472 -3723 -5945 -3523    NA
[10,] -2601 -3759 -2609 -2913 -3196 -3434 -3191 -2272    NA
[11,]  -534  -282  -502  -430 -1110  -964 -1338    19    NA
[12,]   143  -426  -554  -359 -1012 -3152  -175  -323    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/1n9hf1229080651.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/2353c1229080651.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/3vhje1229080651.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/4n6av1229080651.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,] 15303 14296 13686 14271 11609 14602 23629 24680 21402.0 17879.0 15607.0
[2,] 19809 18799 17884 17512 16327 16205 26485 31867 27715.5 22875.0 19979.0
[3,] 21863 20462 19594 18564 17640 18614 31726 32860 30646.0 26958.5 23774.5
[4,] 24452 22643 21474 20500 19505 21769 35369 37941 34054.5 29970.5 26662.0
[5,] 26482 25538 23789 22386 21087 22891 36192 38922 34669.0 30197.0 27271.0
       [,12]
[1,] 15626.0
[2,] 19043.0
[3,] 23382.5
[4,] 25625.0
[5,] 26841.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,] 19417.69 18437.49 17703.27 16990.32 15966.25 15683.63 27047.09 29661.03
[2,] 24308.31 22486.51 21484.73 20137.68 19313.75 21544.37 36404.91 36058.97
         [,9]    [,10]    [,11]   [,12]
[1,] 27104.94 22994.85 20041.28 19705.7
[2,] 34187.06 30922.15 27507.72 27059.3

$out
[1] 12948

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(15303, 19809, 21863, 24452, 26482, 14296,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5nftj1229080651.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,] -1809 -1819 -1403 -1456 -270  8929  972 -6055.0 -4472.0 -3759 -1338
[2,] -1460 -1611 -1070 -1310  974  9657 1134 -4910.0 -4458.5 -3315 -1037
[3,] -1010 -1169  -974 -1185 1572 11824 2572 -4128.5 -3742.0 -3052  -518
[4,]  -944  -915  -738  -995 1804 13492 3078 -3870.5 -3675.0 -2605  -356
[5,]  -444  -610  -282  -853 2993 14293 5330 -2440.0 -3523.0 -2272    19
       [,12]
[1,] -1012.0
[2,]  -783.0
[3,]  -392.5
[4,]  -249.0
[5,]   143.0

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

$conf
         [,1]     [,2]       [,3]    [,4]     [,5]      [,6]    [,7]      [,8]
[1,] -1281.76 -1535.56 -1148.8533 -1350.9 1134.867  9804.233 1548.16 -4709.180
[2,]  -738.24  -802.44  -799.1467 -1019.1 2009.133 13843.767 3595.84 -3547.820
          [,9]     [,10]     [,11]      [,12]
[1,] -4179.674 -3448.616 -898.4164 -690.80007
[2,] -3304.326 -2655.384 -137.5836  -94.19993

$out
[1]  6762 -5945 -3152

$group
[1]  7  9 12

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

Warning message:
In bxp(list(stats = c(-1809, -1460, -1010, -944, -444, -1819, -1611,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wvb61229080651.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,] 16205.0 16327.0 17640.0 19505.0 21087.0 20305.0 18420.0 13418.0 11609
[2,] 18361.5 17698.0 19104.0 21621.5 23340.0 22309.0 19710.0 15279.5 13317
[3,] 20301.5 20784.0 23169.5 25646.5 26186.5 25119.0 20853.5 16193.5 14449
[4,] 23967.5 26168.5 29991.0 32406.5 32433.0 31618.5 25921.5 19823.5 19466
[5,] 31896.0 31867.0 35640.0 38633.0 38922.0 37941.0 32860.0 26485.0 24680

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 17744.57 16920.55 18203.87 20727.39 22039.12 20872.87 18020.39 14120.95
[2,] 22858.43 24647.45 28135.13 30565.61 30333.88 29365.13 23686.61 18266.05
         [,9]
[1,] 11014.08
[2,] 17883.92

$out
[1] 27457

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(16205, 18361.5, 20301.5, 23967.5, 31896, 16327,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/73nhl1229080651.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,] 17221.11 17640.00 17916.00
[2,] 18955.72 19104.00 19342.50
[3,] 21970.11 22622.75 22276.88
[4,] 28036.69 28802.25 28657.44
[5,] 33321.78 32860.00 34904.00

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 17828.22 18199.31 18028.27
[2,] 26112.00 27046.19 26525.48

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(17221.1111111111, 18955.7222222222, 21970.1111111111,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1n9hf1229080651.ps tmp/1n9hf1229080651.png")
> system("convert tmp/2353c1229080651.ps tmp/2353c1229080651.png")
> system("convert tmp/3vhje1229080651.ps tmp/3vhje1229080651.png")
> system("convert tmp/4n6av1229080651.ps tmp/4n6av1229080651.png")
> system("convert tmp/5nftj1229080651.ps tmp/5nftj1229080651.png")
> system("convert tmp/6wvb61229080651.ps tmp/6wvb61229080651.png")
> system("convert tmp/73nhl1229080651.ps tmp/73nhl1229080651.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.224   0.964   1.475