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(17681170.13,19858875.65,16997477.78,16969881.42,18908869.11,15692144.17,15159951.13,15806842.33,16007123.26,16059123.25,16189383.57,12522497.94,14733828.17,15686348.94,13779681.49,14423755.53,15290621.44,14308336.46,13855616.12,14384486.00,15638580.70,19711553.27,20359793.22,16141449.72,20056915.06,20605542.58,19325754.03,20547653.75,19211178.55,19009453.59,18746813.27,16471529.53,18957217.20,20515191.95,18374420.60,16192909.22,18147463.68,19301440.71,18344657.85,17183629.01,19629987.60,17167191.42,17428458.32,16016524.60,18466459.42,18406552.23,18174068.44,14851949.20,16260733.42,18329610.38,18003781.65,15903762.33,19554176.37,16554237.93,16198892.67,16571771.60,17535166.38,16198106.13,17487530.67,13768040.14)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
          [,1]     [,2]     [,3]     [,4]     [,5] [,6]
 [1,] 17681170 14733828 20056915 18147464 16260733   NA
 [2,] 19858876 15686349 20605543 19301441 18329610   NA
 [3,] 16997478 13779681 19325754 18344658 18003782   NA
 [4,] 16969881 14423756 20547654 17183629 15903762   NA
 [5,] 18908869 15290621 19211179 19629988 19554176   NA
 [6,] 15692144 14308336 19009454 17167191 16554238   NA
 [7,] 15159951 13855616 18746813 17428458 16198893   NA
 [8,] 15806842 14384486 16471530 16016525 16571772   NA
 [9,] 16007123 15638581 18957217 18466459 17535166   NA
[10,] 16059123 19711553 20515192 18406552 16198106   NA
[11,] 16189384 20359793 18374421 18174068 17487531   NA
[12,] 12522498 16141450 16192909 14851949 13768040   NA
> darr
             [,1]       [,2]       [,3]        [,4]       [,5] [,6]
 [1,]  2177705.52   952520.8   548627.5  1153977.03  2068877.0   NA
 [2,] -2861397.87 -1906667.4 -1279788.5  -956782.86  -325828.7   NA
 [3,]   -27596.36   644074.0  1221899.7 -1161028.84 -2100019.3   NA
 [4,]  1938987.69   866865.9 -1336475.2  2446358.59  3650414.0   NA
 [5,] -3216724.94  -982285.0  -201725.0 -2462796.18 -2999938.4   NA
 [6,]  -532193.04  -452720.3  -262640.3   261266.90  -355345.3   NA
 [7,]   646891.20   528869.9 -2275283.7 -1411933.72   372878.9   NA
 [8,]   200280.93  1254094.7  2485687.7  2449934.82   963394.8   NA
 [9,]    51999.99  4072972.6  1557974.8   -59907.19 -1337060.2   NA
[10,]   130260.32   648239.9 -2140771.3  -232483.79  1289424.5   NA
[11,] -3666885.63 -4218343.5 -2181511.4 -3322119.24 -3719490.5   NA
[12,]  2211330.23  3915465.3  1954554.5  1408784.22         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/19kj91289844401.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/29kj91289844401.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/39kj91289844401.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/42bic1289844401.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]
[1,] 14733828 18329610 16997478 14423756 18908869 14308336 13855616 15806842
[2,] 16260733 18329610 16997478 15903762 18908869 15692144 15159951 15806842
[3,] 17681170 19301441 18003782 16969881 19211179 16554238 16198893 16016525
[4,] 18147464 19858876 18344658 17183629 19554176 17167191 17428458 16471530
[5,] 20056915 20605543 19325754 17183629 19629988 19009454 18746813 16571772
         [,9]    [,10]    [,11]    [,12]
[1,] 15638581 16059123 16189384 12522498
[2,] 16007123 16198106 17487531 13768040
[3,] 17535166 18406552 18174068 14851949
[4,] 18466459 19711553 18374421 16141450
[5,] 18957217 20515192 18374421 16192909

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 16348011 18220866 17051868 16065531 18755206 15511973 14595971 15546858
[2,] 19014329 20382016 18955696 17874232 19667151 17596503 17801814 16486191
         [,9]    [,10]    [,11]    [,12]
[1,] 15797406 15923959 17547394 13174904
[2,] 19272927 20889145 18800743 16528994

$out
[1] 15686349 13779681 20547654 15290621 14384486 20359793

$group
[1]  2  3  4  5  8 11

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

Warning message:
In bxp(list(stats = c(14733828.17, 16260733.42, 17681170.13, 18147463.68,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/52bic1289844401.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]
[1,]  548627.5 -2861397.9 -2100019.32 -1336475.2 -3216724.9 -532193.0
[2,]  952520.8 -1906667.4 -1161028.84   866865.9 -2999938.4 -452720.3
[3,] 1153977.0 -1279788.5   -27596.36  1938987.7 -2462796.2 -355345.3
[4,] 2068877.0  -956782.9   644074.04  2446358.6  -982285.0 -262640.3
[5,] 2177705.5  -325828.7  1221899.72  3650414.0  -201725.0 -262640.3
           [,7]      [,8]        [,9]     [,10]    [,11]   [,12]
[1,] -2275283.7  200280.9 -1337060.25 -232483.8 -4218343 1408784
[2,] -1411933.7  963394.8   -59907.19 -232483.8 -3719491 1681669
[3,]   372878.9 1254094.7    51999.99  130260.3 -3666886 2082942
[4,]   528869.9 2449934.8  1557974.75  648239.9 -3322119 3063398
[5,]   646891.2 2485687.7  1557974.75 1289424.5 -3322119 3915465

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 4

$conf
          [,1]       [,2]     [,3]      [,4]     [,5]      [,6]    [,7]
[1,]  365162.6 -1950974.6 -1303078  822922.1 -3888465 -489655.3 -998488
[2,] 1942791.5  -608602.5  1247885 3055053.2 -1037127 -221035.2 1744246
          [,8]     [,9]     [,10]    [,11]     [,12]
[1,]  203709.3 -1091191 -492056.9 -3947667  991376.9
[2,] 2304480.1  1195191  752577.5 -3386104 3174507.8

$out
[1]   261266.9  4072972.6 -2140771.3 -2181511.4

$group
[1]  6  9 10 11

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

Warning message:
In bxp(list(stats = c(548627.52, 952520.77, 1153977.0300, 2068876.96,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6vl0f1289844401.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]
[1,] 15159951 13779681 16192909 16016525 15903762   NA
[2,] 15749493 14346411 18560617 17175410 16198499   NA
[3,] 16124253 15012225 19110316 18160766 16563005   NA
[4,] 17339324 15913899 20286054 18436506 17769474   NA
[5,] 18908869 16141450 20605543 19629988 19554176   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 15399121 14297283 18323333 17585572 15846473   NA
[2,] 16849386 15727167 19897299 18735960 17279537   NA

$out
[1] 19858876 12522498 19711553 20359793 14851949 13768040

$group
[1] 1 1 2 2 4 5

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

Warning message:
In bxp(list(stats = c(15159951.13, 15749493.25, 16124253.41, 17339323.955,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7vl0f1289844401.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,] 14695369 14851949 14954745
[2,] 16412110 16376565 16361936
[3,] 17305590 17608168 17220445
[4,] 18147572 18290310 17942903
[5,] 18756364 19301441 19231523

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 16514034 16735296 16499356
[2,] 18097146 18481040 17941534

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(14695369.244, 16412109.508, 17305589.976,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/19kj91289844401.ps tmp/19kj91289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/29kj91289844401.ps tmp/29kj91289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/39kj91289844401.ps tmp/39kj91289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/42bic1289844401.ps tmp/42bic1289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/52bic1289844401.ps tmp/52bic1289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/6vl0f1289844401.ps tmp/6vl0f1289844401.png",intern=TRUE))
character(0)
> try(system("convert tmp/7vl0f1289844401.ps tmp/7vl0f1289844401.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.160   0.903   1.610