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

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

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

$out
[1] 20359793 14384486 15290621 20547654 13779681 15686349

$group
[1]  2  5  8  9 10 11

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

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

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

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

$out
[1]  2181511.4  2140771.3 -4072972.6  -261266.9

$group
[1] 1 2 3 6

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

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

$n
[1] 12 12 12 12 12  0

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

$out
[1] 13768040 14851949 20359793 19711553 12522498 19858876

$group
[1] 1 2 4 4 5 5

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

Warning message:
In bxp(list(stats = c(15903762.33, 16198499.4, 16563004.765, 17769474.015,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7bzj51291129230.ps",horizontal=F,onefile=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/1qg3z1291129230.ps tmp/1qg3z1291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/2qg3z1291129230.ps tmp/2qg3z1291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jqk21291129230.ps tmp/3jqk21291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jqk21291129230.ps tmp/4jqk21291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/5jqk21291129230.ps tmp/5jqk21291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bzj51291129230.ps tmp/6bzj51291129230.png",intern=TRUE))
character(0)
> try(system("convert tmp/7bzj51291129230.ps tmp/7bzj51291129230.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.928   1.471   2.145