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(2357402,2199181,2603060,2629720,2638792,2717481,3810804,3871664,2998364,2923432,2712359,2996099,2395029,2483862,3120231,3360606,3177203,3062783,4242509,4026394,3192481,3118695,2782482,3209833,2630190,2592882,3785309,3231539,3421369,3312134,4647303,4289177,3463853,3304422,3006121,3464238,2921118,2624018,3500718,3939351,3467672,3343628,4852445,4597807,3653145,3572079,3334861,3695369,3075704,2852998,3942704,4004560,3822145,3760085,5267816,5271333,4144142,4109749,3896808,4211074,3402318,3279817,4706628,4079499,4344530,4048625,5394915,5611967,4145481,4025610,3552218,3910443)
> 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] 72
> (np <- floor(n / par1))
[1] 6
> 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] 6 6 6 6 6 6 6 6 6 6 6 6
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6] [,7]
 [1,] 2357402 2395029 2630190 2921118 3075704 3402318   NA
 [2,] 2199181 2483862 2592882 2624018 2852998 3279817   NA
 [3,] 2603060 3120231 3785309 3500718 3942704 4706628   NA
 [4,] 2629720 3360606 3231539 3939351 4004560 4079499   NA
 [5,] 2638792 3177203 3421369 3467672 3822145 4344530   NA
 [6,] 2717481 3062783 3312134 3343628 3760085 4048625   NA
 [7,] 3810804 4242509 4647303 4852445 5267816 5394915   NA
 [8,] 3871664 4026394 4289177 4597807 5271333 5611967   NA
 [9,] 2998364 3192481 3463853 3653145 4144142 4145481   NA
[10,] 2923432 3118695 3304422 3572079 4109749 4025610   NA
[11,] 2712359 2782482 3006121 3334861 3896808 3552218   NA
[12,] 2996099 3209833 3464238 3695369 4211074 3910443   NA
> darr
         [,1]    [,2]    [,3]    [,4]     [,5]     [,6] [,7]
 [1,] -158221   88833  -37308 -297100  -222706  -122501   NA
 [2,]  403879  636369 1192427  876700  1089706  1426811   NA
 [3,]   26660  240375 -553770  438633    61856  -627129   NA
 [4,]    9072 -183403  189830 -471679  -182415   265031   NA
 [5,]   78689 -114420 -109235 -124044   -62060  -295905   NA
 [6,] 1093323 1179726 1335169 1508817  1507731  1346290   NA
 [7,]   60860 -216115 -358126 -254638     3517   217052   NA
 [8,] -873300 -833913 -825324 -944662 -1127191 -1466486   NA
 [9,]  -74932  -73786 -159431  -81066   -34393  -119871   NA
[10,] -211073 -336213 -298301 -237218  -212941  -473392   NA
[11,]  283740  427351  458117  360508   314266   358225   NA
[12,] -601070 -579643 -543120 -619665  -808756       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/13h3v1259518854.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/26ljh1259518854.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/3kizs1259518854.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/4gpia1259518854.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,] 2357402 2199181 2603060 2629720 2638792 2717481 3810804 3871664 2998364
[2,] 2395029 2483862 3120231 3231539 3177203 3062783 4242509 4026394 3192481
[3,] 2775654 2608450 3643014 3649978 3444520 3327881 4749874 4443492 3558499
[4,] 3075704 2852998 3942704 4004560 3822145 3760085 5267816 5271333 4144142
[5,] 3402318 3279817 4706628 4079499 4344530 4048625 5394915 5611967 4145481
       [,10]   [,11]   [,12]
[1,] 2923432 2712359 2996099
[2,] 3118695 2782482 3209833
[3,] 3438250 3170491 3579804
[4,] 4025610 3552218 3910443
[5,] 4109749 3896808 4211074

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

$conf
        [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
[1,] 2336597 2370345 3112492 3151355 3028512 2878099 4088518 3640466 2944647
[2,] 3214711 2846555 4173535 4148602 3860529 3777663 5411230 5246518 4172351
       [,10]   [,11]   [,12]
[1,] 2853261 2673986 3127887
[2,] 4023240 3666996 4031720

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2357402, 2395029, 2775654, 3075704, 3402318,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/59wrj1259518854.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]
[1,] -297100  403879 -627129 -471679.0 -124044.0 1093323 -358126 -1466486
[2,] -222706  636369 -553770 -183403.0 -124044.0 1179726 -254638 -1127191
[3,] -140361  983203   44258  -86671.5 -111827.5 1340730 -106299  -908981
[4,]  -37308 1192427  240375  189830.0  -62060.0 1507731   60860  -833913
[5,]   88833 1426811  438633  265031.0  -62060.0 1508817  217052  -825324
        [,9]     [,10]    [,11]   [,12]
[1,] -159431 -473392.0 283740.0 -619665
[2,] -119871 -336213.0 314266.0 -619665
[3,]  -77999 -267759.5 359366.5 -601070
[4,]  -73786 -212941.0 427351.0 -579643
[5,]  -34393 -211073.0 458117.0 -543120

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

$conf
          [,1]      [,2]      [,3]      [,4]       [,5]    [,6]      [,7]
[1,] -259948.7  624527.6 -467991.2 -327418.8 -151809.18 1129156 -309805.4
[2,]  -20773.3 1341878.4  556507.2  154075.8  -71845.82 1552303   97207.4
           [,8]       [,9]     [,10]  [,11]     [,12]
[1,] -1098154.8 -107725.31 -347273.9 286423 -629349.4
[2,]  -719807.2  -48272.69 -188245.1 432310 -572790.6

$out
[1]   78689 -295905 -808756

$group
[1]  5  5 12

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

Warning message:
In bxp(list(stats = c(-297100, -222706, -140361, -37308, 88833,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/638yp1259518854.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]
[1,] 2199181 2395029 2592882 2624018 3760085 3279817   NA
[2,] 2616390 2922632 3118830 3339244 3791115 3731330   NA
[3,] 2714920 3148717 3366752 3536398 3973632 4064062   NA
[4,] 2997232 3285220 3624774 3817360 4177608 4525579   NA
[5,] 2998364 3360606 4289177 3939351 4211074 5611967   NA

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

$conf
        [,1]    [,2]    [,3]    [,4]    [,5]    [,6] [,7]
[1,] 2541216 2983339 3135987 3318327 3797350 3701800   NA
[2,] 2888624 3314095 3597516 3754470 4149914 4426324   NA

$out
 [1] 3810804 3871664 4242509 4026394 4647303 4852445 4597807 3075704 2852998
[10] 5267816 5271333

$group
 [1] 1 1 2 2 3 4 4 5 5 5 5

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

Warning message:
In bxp(list(stats = c(2199181, 2616390, 2714920, 2997231.5, 2998364,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/78e0u1259518854.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,] 3214142 2775654 3168135
[2,] 3294132 3249186 3279340
[3,] 3524938 3501510 3549016
[4,] 3604676 3646496 3633445
[5,] 3609775 3649978 3640858

$n
[1] 12 12 12

$conf
        [,1]    [,2]    [,3]
[1,] 3383297 3320294 3387506
[2,] 3666580 3682726 3710525

$out
 [1] 2796960 2672126 4702632 4611390 2608450 4749874 4443492 2745438 2653435
[10] 4753840 4597521

$group
 [1] 1 1 1 1 2 2 2 3 3 3 3

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

Warning message:
In bxp(list(stats = c(3214141.5, 3294132.08333333, 3524938.5, 3604676.33333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/13h3v1259518854.ps tmp/13h3v1259518854.png")
> system("convert tmp/26ljh1259518854.ps tmp/26ljh1259518854.png")
> system("convert tmp/3kizs1259518854.ps tmp/3kizs1259518854.png")
> system("convert tmp/4gpia1259518854.ps tmp/4gpia1259518854.png")
> system("convert tmp/59wrj1259518854.ps tmp/59wrj1259518854.png")
> system("convert tmp/638yp1259518854.ps tmp/638yp1259518854.png")
> system("convert tmp/78e0u1259518854.ps tmp/78e0u1259518854.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.221   0.954   1.512