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(8365275,8336926,8300788,8268641,8227829,8171966,8121423,8081957,8042293,8011566,7992324,7976789,7968041,7959017,7948278,7936118,7905633,7840709,7754891,7677850,7619567,7585317,7574586,7569794,7564985,7561434,7561910,7574140,7568710,7549433,7549425,7562305,7568430,7565525,7578903,7599038,7586115,7544201,7500482,7467086,7441055,7415403,7376998,7322066,7270889,7223801,7175811,7129864,7086299,7047539)
> 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] 50
> (np <- floor(n / par1))
[1] 4
> 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 4 4 4 4 4 4 4 4 4 4
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]
 [1,] 8365275 7968041 7564985 7586115 7086299
 [2,] 8336926 7959017 7561434 7544201 7047539
 [3,] 8300788 7948278 7561910 7500482      NA
 [4,] 8268641 7936118 7574140 7467086      NA
 [5,] 8227829 7905633 7568710 7441055      NA
 [6,] 8171966 7840709 7549433 7415403      NA
 [7,] 8121423 7754891 7549425 7376998      NA
 [8,] 8081957 7677850 7562305 7322066      NA
 [9,] 8042293 7619567 7568430 7270889      NA
[10,] 8011566 7585317 7565525 7223801      NA
[11,] 7992324 7574586 7578903 7175811      NA
[12,] 7976789 7569794 7599038 7129864      NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]
 [1,] -28349  -9024  -3551 -41914 -38760
 [2,] -36138 -10739    476 -43719     NA
 [3,] -32147 -12160  12230 -33396     NA
 [4,] -40812 -30485  -5430 -26031     NA
 [5,] -55863 -64924 -19277 -25652     NA
 [6,] -50543 -85818     -8 -38405     NA
 [7,] -39466 -77041  12880 -54932     NA
 [8,] -39664 -58283   6125 -51177     NA
 [9,] -30727 -34250  -2905 -47088     NA
[10,] -19242 -10731  13378 -47990     NA
[11,] -15535  -4792  20135 -45947     NA
[12,]  -8748  -4809 -12923 -43565     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/1j3101290712158.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/freestat/rcomp/tmp/2bu0l1290712158.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/freestat/rcomp/tmp/3bu0l1290712158.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/freestat/rcomp/tmp/444i61290712158.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,] 7086299 7047539 7500482 7467086 7441055 7415403 7376998 7322066 7270889
[2,] 7564985 7544201 7531196 7520613 7504882 7482418 7463212 7442186 7419660
[3,] 7586115 7561434 7755094 7755129 7737172 7695071 7652158 7620078 7593998
[4,] 7968041 7959017 8124533 8102380 8066731 8006338 7938157 7879904 7830930
[5,] 8365275 8336926 8300788 8268641 8227829 8171966 8121423 8081957 8042293
       [,10]   [,11]   [,12]
[1,] 7223801 7175811 7129864
[2,] 7394663 7375198 7349829
[3,] 7575421 7576744 7584416
[4,] 7798442 7785614 7787914
[5,] 8011566 7992324 7976789

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

$conf
        [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
[1,] 7301317 7268326 7286358 7295533 7293311 7281175 7276951 7274280 7269095
[2,] 7870913 7854542 8223830 8214725 8181032 8108967 8027365 7965875 7918902
       [,10]   [,11]   [,12]
[1,] 7256436 7252517 7238329
[2,] 7894406 7900972 7930503

$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(7086299, 7564985, 7586115, 7968041, 8365275,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/544i61290712158.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,] -41914 -43719.0 -33396.0 -40812.0 -64924.0 -85818.0 -77041.0 -58283.0
[2,] -38760 -39928.5 -32771.5 -35648.5 -60393.5 -68180.5 -65986.5 -54730.0
[3,] -28349 -23438.5 -22153.5 -28258.0 -40757.5 -44474.0 -47199.0 -45420.5
[4,]  -9024  -5131.5     35.0 -15730.5 -22464.5 -19206.5 -13293.0 -16769.5
[5,]  -3551    476.0  12230.0  -5430.0 -19277.0     -8.0  12880.0   6125.0
         [,9]    [,10]    [,11]    [,12]
[1,] -47088.0 -47990.0 -45947.0 -43565.0
[2,] -40669.0 -33616.0 -30741.0 -28244.0
[3,] -32488.5 -14986.5 -10163.5 -10835.5
[4,] -16816.0   1323.5   7671.5  -6778.5
[5,]  -2905.0  13378.0  20135.0  -4809.0

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

$conf
           [,1]      [,2]       [,3]      [,4]      [,5]      [,6]       [,7]
[1,] -49360.383 -50928.13 -48070.635 -43993.22 -70721.41 -83163.46 -88826.865
[2,]  -7337.617   4051.13   3763.635 -12522.78 -10793.59  -5784.54  -5571.135
          [,8]      [,9]     [,10]     [,11]      [,12]
[1,] -75409.29 -51332.37 -42588.71 -40509.38 -27793.245
[2,] -15431.70 -13644.63  12615.71  20182.38   6122.245

$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(-41914, -38760, -28349, -9024, -3551, -43719,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/644i61290712158.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]
[1,] 7976789 7569794 7549425 7129864 7047539
[2,] 8026930 7602442 7561672 7247345 7047539
[3,] 8146694 7797800 7565255 7396200 7066919
[4,] 8284714 7942198 7571425 7483784 7086299
[5,] 8365275 7968041 7578903 7586115 7086299

$n
[1] 12 12 12 12  2

$conf
        [,1]    [,2]    [,3]    [,4]    [,5]
[1,] 8029117 7642835 7560807 7288359 7023615
[2,] 8264272 7952765 7569703 7504042 7110223

$out
[1] 7599038

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(7976789, 8026929.5, 8146694.5, 8284714.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7fdh91290712158.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,] 7568871 7561434 7576644
[2,] 7610924 7580580 7597817
[3,] 7695254 7607038 7698073
[4,] 7765092 7716121 7764001
[5,] 7827864 7755129 7791479

$n
[1] 12 12 12

$conf
        [,1]    [,2]    [,3]
[1,] 7624936 7545217 7622275
[2,] 7765571 7668859 7773871

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(7568871.25, 7610923.5, 7695253.825, 7765092.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1j3101290712158.ps tmp/1j3101290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/2bu0l1290712158.ps tmp/2bu0l1290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/3bu0l1290712158.ps tmp/3bu0l1290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/444i61290712158.ps tmp/444i61290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/544i61290712158.ps tmp/544i61290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/644i61290712158.ps tmp/644i61290712158.png",intern=TRUE))
character(0)
> try(system("convert tmp/7fdh91290712158.ps tmp/7fdh91290712158.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.064   1.617   2.448