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.

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(43657,42811,45419,50846,54500,51035,38675,36214,38763,39486,40540,40719,40471,39947,42683,47090,51520,48823,36122,33812,36928,37737,40123,41713,42025,42169,46352,50939,56139,52713,38532,37860,40880,41988,44576,46728,46913,49357,54709,60819,63695,60109,45544,43596,44431,45575,47980,49211,51374,52954,57529,62960,64530,61008,44964,43480,45429,47616,49364,51010,53188,55317,60106,65845,67028,63617,47605,45844,47925,50156,52258,53476,54327,55214,59347,64718,66208,62744,45587,43684,45676,47088,48907,50964,51798)
> 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] 85
> (np <- floor(n / par1))
[1] 7
> 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] 8 7 7 7 7 7 7 7 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 43657 40471 42025 46913 51374 53188 54327 51798
 [2,] 42811 39947 42169 49357 52954 55317 55214    NA
 [3,] 45419 42683 46352 54709 57529 60106 59347    NA
 [4,] 50846 47090 50939 60819 62960 65845 64718    NA
 [5,] 54500 51520 56139 63695 64530 67028 66208    NA
 [6,] 51035 48823 52713 60109 61008 63617 62744    NA
 [7,] 38675 36122 38532 45544 44964 47605 45587    NA
 [8,] 36214 33812 37860 43596 43480 45844 43684    NA
 [9,] 38763 36928 40880 44431 45429 47925 45676    NA
[10,] 39486 37737 41988 45575 47616 50156 47088    NA
[11,] 40540 40123 44576 47980 49364 52258 48907    NA
[12,] 40719 41713 46728 49211 51010 53476 50964    NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7] [,8]
 [1,]   -846   -524    144   2444   1580   2129    887   NA
 [2,]   2608   2736   4183   5352   4575   4789   4133   NA
 [3,]   5427   4407   4587   6110   5431   5739   5371   NA
 [4,]   3654   4430   5200   2876   1570   1183   1490   NA
 [5,]  -3465  -2697  -3426  -3586  -3522  -3411  -3464   NA
 [6,] -12360 -12701 -14181 -14565 -16044 -16012 -17157   NA
 [7,]  -2461  -2310   -672  -1948  -1484  -1761  -1903   NA
 [8,]   2549   3116   3020    835   1949   2081   1992   NA
 [9,]    723    809   1108   1144   2187   2231   1412   NA
[10,]   1054   2386   2588   2405   1748   2102   1819   NA
[11,]    179   1590   2152   1231   1646   1218   2057   NA
[12,]   -248    312    185   2163   2178    851    834   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/1c1nu1272027565.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/2c1nu1272027565.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/3c1nu1272027565.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/45smw1272027565.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,] 40471.0 39947 42683.0 47090.0 51520.0 48823 36122.0 33812 36928.0 37737
[2,] 42841.0 42490 45885.5 50892.5 55319.5 51874 38603.5 37037 39821.5 40737
[3,] 49143.5 49357 54709.0 60819.0 63695.0 60109 44964.0 43480 44431.0 45575
[4,] 52493.0 54084 58438.0 63839.0 65369.0 61876 45565.5 43640 45552.5 47352
[5,] 54327.0 55317 60106.0 65845.0 67028.0 63617 47605.0 45844 47925.0 50156
       [,11]   [,12]
[1,] 40123.0 40719.0
[2,] 42558.0 44220.5
[3,] 47980.0 49211.0
[4,] 49135.5 50987.0
[5,] 52258.0 53476.0

$n
 [1] 8 7 7 7 7 7 7 7 7 7 7 7

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 43751.75 42433.25 47212.85 53087.56 57693.6 54135.97 40806.41 39536.79
[2,] 54535.25 56280.75 62205.15 68550.44 69696.4 66082.03 49121.59 47423.21
         [,9]    [,10]    [,11]    [,12]
[1,] 41008.54 41624.63 44052.02 45170.16
[2,] 47853.46 49525.37 51907.98 53251.84

$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(40471, 42841, 49143.5, 52493, 54327, 39947,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/55smw1272027565.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]  [,10]  [,11]
[1,] -846.0 2608.0 4407 1183 -3586.0 -17157 -2461.0  835.0  723.0 1054.0 1218.0
[2,] -190.0 3434.5 4979 1530 -3493.5 -16028 -2129.0 1970.5  958.5 1783.5 1224.5
[3,]  887.0 4183.0 5427 2876 -3464.0 -14565 -1903.0 2081.0 1144.0 2102.0 1590.0
[4,] 1854.5 4682.0 5585 4042 -3418.5 -13441 -1622.5 2784.5 1799.5 2395.5 1851.5
[5,] 2444.0 5352.0 6110 5200 -3411.0 -12360 -1484.0 3116.0 2231.0 2588.0 2152.0
      [,12]
[1,] -248.0
[2,]  248.5
[3,]  834.0
[4,] 1507.0
[5,] 2178.0

$n
 [1] 7 7 7 7 7 7 7 7 7 7 7 7

$conf
          [,1]     [,2]     [,3]     [,4]      [,5]      [,6]      [,7]
[1,] -333.9424 3438.013 5065.107 1375.874 -3508.789 -16109.91 -2205.474
[2,] 2107.9424 4927.987 5788.893 4376.126 -3419.211 -13020.09 -1600.526
         [,8]      [,9]    [,10]    [,11]     [,12]
[1,] 1594.892  641.7684 1736.523 1215.566   82.4441
[2,] 2567.108 1646.2316 2467.477 1964.434 1585.5559

$out
[1] -2697  -672   179

$group
[1]  5  7 11

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

Warning message:
In bxp(list(stats = c(-846, -190, 887, 1854.5, 2444, 2608, 3434.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6x14z1272027565.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]  [,8]
[1,] 36214.0 33812.0 37860.0 43596.0 43480.0 45844.0 43684.0 51798
[2,] 39124.5 37332.5 41434.0 45559.5 46522.5 49040.5 46382.0 51798
[3,] 41765.0 40297.0 43372.5 48595.5 51192.0 53332.0 52645.5 51798
[4,] 48132.5 44886.5 48833.5 57409.0 59268.5 61861.5 61045.5 51798
[5,] 54500.0 51520.0 56139.0 63695.0 64530.0 67028.0 66208.0 51798

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]  [,8]
[1,] 37656.39 36851.57 39997.54 43190.86 45378.46 47484.26 45957.38 51798
[2,] 45873.61 43742.43 46747.46 54000.14 57005.54 59179.74 59333.62 51798

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(36214, 39124.5, 41765, 48132.5, 54500, 33812,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7x14z1272027565.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,] 40641.43 43480.00 40338.50
[2,] 43548.43 45269.50 43365.75
[3,] 47828.92 49177.25 47650.50
[4,] 54728.14 57409.00 54518.38
[5,] 60517.14 63695.00 60344.25

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 42729.78 43640.34 42563.71
[2,] 52928.06 54714.16 52737.29

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(40641.4285714286, 43548.4285714286, 47828.9196428571,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1c1nu1272027565.ps tmp/1c1nu1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/2c1nu1272027565.ps tmp/2c1nu1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/3c1nu1272027565.ps tmp/3c1nu1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/45smw1272027565.ps tmp/45smw1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/55smw1272027565.ps tmp/55smw1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/6x14z1272027565.ps tmp/6x14z1272027565.png",intern=TRUE))
character(0)
> try(system("convert tmp/7x14z1272027565.ps tmp/7x14z1272027565.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.886   1.479   2.210