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(990633,1047696,835567,867169,1161937,846354,852288,1033321,865605,1447252,187910,631074,630867,774318,637859,764244,591077,6990528,693887,807691,802010,1303548,162242,683606,957467,945691,904843,908178,1107174,1049559,960388,1086083,1162190,1568232,201733,871654,963424,930097,940654,883223,1197679,955581,886167,1180383,954937,1589668,291111,800423,873241,906223,768775,833106,1113844,940493,837936,1165387,909788,1584838)
> 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] 58
> (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 5 5 5 5 5 5 5 5 4 4
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]
 [1,]  990633  630867  957467  963424  873241
 [2,] 1047696  774318  945691  930097  906223
 [3,]  835567  637859  904843  940654  768775
 [4,]  867169  764244  908178  883223  833106
 [5,] 1161937  591077 1107174 1197679 1113844
 [6,]  846354 6990528 1049559  955581  940493
 [7,]  852288  693887  960388  886167  837936
 [8,] 1033321  807691 1086083 1180383 1165387
 [9,]  865605  802010 1162190  954937  909788
[10,] 1447252 1303548 1568232 1589668 1584838
[11,]  187910  162242  201733  291111      NA
[12,]  631074  683606  871654  800423      NA
> darr
          [,1]     [,2]     [,3]     [,4]    [,5]
 [1,]    57063   143451   -11776   -33327   32982
 [2,]  -212129  -136459   -40848    10557 -137448
 [3,]    31602   126385     3335   -57431   64331
 [4,]   294768  -173167   198996   314456  280738
 [5,]  -315583  6399451   -57615  -242098 -173351
 [6,]     5934 -6296641   -89171   -69414 -102557
 [7,]   181033   113804   125695   294216  327451
 [8,]  -167716    -5681    76107  -225446 -255599
 [9,]   581647   501538   406042   634731  675050
[10,] -1259342 -1141306 -1366499 -1298557      NA
[11,]   443164   521364   669921   509312      NA
[12,]     -207   273861    91770    72818      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/17g8e1290532193.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/27g8e1290532193.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/37g8e1290532193.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/4i77z1290532193.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]
[1,] 873241 906223 637859 764244 1107174  846354 837936 1033321 802010 1303548
[2,] 873241 906223 768775 833106 1107174  940493 837936 1033321 865605 1447252
[3,] 957467 930097 835567 867169 1113844  955581 852288 1086083 909788 1568232
[4,] 963424 945691 904843 883223 1161937 1049559 886167 1165387 954937 1584838
[5,] 990633 945691 940654 908178 1197679 1049559 886167 1180383 954937 1589668
        [,11]    [,12]
[1,] 162242.0 631074.0
[2,] 175076.0 657340.0
[3,] 194821.5 742014.5
[4,] 246422.0 836038.5
[5,] 291111.0 871654.0

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

$conf
        [,1]   [,2]     [,3]     [,4]    [,5]      [,6]     [,7]      [,8]
[1,]  893744 902209 739421.7 831756.5 1075149  878515.2 818208.1  992765.5
[2,] 1021190 957985 931712.3 902581.5 1152539 1032646.8 886367.9 1179400.5
         [,9]   [,10]    [,11]    [,12]
[1,] 846666.2 1471014 138458.2 600842.7
[2,] 972909.8 1665450 251184.8 883186.3

$out
[1]  630867 1047696  774318  591077 6990528  693887  960388  807691 1162190

$group
[1] 1 2 2 5 6 7 7 8 9

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

Warning message:
In bxp(list(stats = c(873241, 873241, 957467, 963424, 990633, 906223,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5i77z1290532193.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]
[1,] -33327 -212129 -57431 198996 -315583 -102557 113804 -255599 406042
[2,] -11776 -137448   3335 198996 -242098 -102557 125695 -225446 501538
[3,]  32982 -136459  31602 280738 -173351  -89171 181033 -167716 581647
[4,]  57063  -40848  64331 294768  -57615  -69414 294216   -5681 634731
[5,] 143451   10557 126385 314456  -57615  -69414 327451   76107 675050
        [,10]    [,11]    [,12]
[1,] -1366499 443164.0   -207.0
[2,] -1332528 476238.0  36305.5
[3,] -1278950 515338.0  82294.0
[4,] -1200324 595642.5 182815.5
[5,] -1141306 669921.0 273861.0

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

$conf
          [,1]       [,2]      [,3]     [,4]       [,5]       [,6]      [,7]
[1,] -15659.46 -204716.32 -11497.62 213065.7 -303706.22 -112589.76  61956.49
[2,]  81623.46  -68201.68  74701.62 348410.3  -42995.78  -65752.24 300109.51
           [,8]     [,9]    [,10]    [,11]    [,12]
[1,] -323001.40 487533.2 -1383391 421008.4 -33448.9
[2,]  -12430.60 675760.8 -1174508 609667.6 198036.9

$out
[1]  -173167  6399451     5934 -6296641

$group
[1] 4 5 6 6

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

Warning message:
In bxp(list(stats = c(-33327, -11776, 32982, 57063, 143451, -212129,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6bgpk1290532193.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,]  631074.0 591077.0  871654.0  800423.0  768775.0
[2,]  840960.5 634363.0  906510.5  884695.0  837936.0
[3,]  866387.0 729065.5  958927.5  947795.5  908005.5
[4,] 1040508.5 804850.5 1096628.5 1071903.5 1113844.0
[5,] 1161937.0 807691.0 1162190.0 1197679.0 1165387.0

$n
[1] 12 12 12 12 10

$conf
         [,1]   [,2]      [,3]      [,4]      [,5]
[1,] 775371.8 651305  872213.4  862408.4  770150.9
[2,] 957402.2 806826 1045641.6 1033182.6 1045860.1

$out
 [1] 1447252  187910 6990528 1303548  162242 1568232  201733 1589668  291111
[10] 1584838

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

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

Warning message:
In bxp(list(stats = c(631074, 840960.5, 866387, 1040508.5, 1161937,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7bgpk1290532193.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,]  746689.2  742014.5  744351.9
[2,]  831836.4  843927.5  847486.8
[3,]  901965.7  919942.5  914301.8
[4,] 1044457.6 1021775.0 1047190.0
[5,] 1054573.0 1113844.0 1134555.5

$n
[1] 12 12 12

$conf
         [,1]    [,2]      [,3]
[1,] 804987.7  838825  823215.8
[2,] 998943.7 1001060 1005387.7

$out
[1] 2156503.0 1498707.6  210749.0 1568232.0  194821.5 1516045.0  202785.2

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

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

Warning message:
In bxp(list(stats = c(746689.25, 831836.4, 901965.7, 1044457.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/17g8e1290532193.ps tmp/17g8e1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/27g8e1290532193.ps tmp/27g8e1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/37g8e1290532193.ps tmp/37g8e1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/4i77z1290532193.ps tmp/4i77z1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/5i77z1290532193.ps tmp/5i77z1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/6bgpk1290532193.ps tmp/6bgpk1290532193.png",intern=TRUE))
character(0)
> try(system("convert tmp/7bgpk1290532193.ps tmp/7bgpk1290532193.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.898   1.435   3.717