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(66857.2,64722.8,68489.6,71342.9,63542.5,69425.0,58927.9,61009.0,66837.0,66147.6,65982.3,65527.5,65914.6,59189.9,66211.4,66400.8,60167.7,64547.9,57706.2,58642.6,60082.1,63414.8,66044.0,57628.5,62838.8,55758.6,61004.5,66173.4,57489.0,59552.2,57061.8,55895.3,56314.7,61232.8,60014.1,57685.4,60403.1,52349.7,55693.3,65676.1,54898.8,55518.2,53779.1,52340.9,55704.4,60330.3,52837.4,55388.1,60383.4,52070.3,54077.0,62887.8,49212.8,57722.0,53936.8,46991.0,54984.2,56485.1,51277.8,53596.4,54252.5,49413.0,53213.2,58695.3,48723.5,54510.0,49454.1,46136.6,54622.5,50583.0,53224.3,53056.4)
> 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,] 66857.2 65914.6 62838.8 60403.1 60383.4 54252.5   NA
 [2,] 64722.8 59189.9 55758.6 52349.7 52070.3 49413.0   NA
 [3,] 68489.6 66211.4 61004.5 55693.3 54077.0 53213.2   NA
 [4,] 71342.9 66400.8 66173.4 65676.1 62887.8 58695.3   NA
 [5,] 63542.5 60167.7 57489.0 54898.8 49212.8 48723.5   NA
 [6,] 69425.0 64547.9 59552.2 55518.2 57722.0 54510.0   NA
 [7,] 58927.9 57706.2 57061.8 53779.1 53936.8 49454.1   NA
 [8,] 61009.0 58642.6 55895.3 52340.9 46991.0 46136.6   NA
 [9,] 66837.0 60082.1 56314.7 55704.4 54984.2 54622.5   NA
[10,] 66147.6 63414.8 61232.8 60330.3 56485.1 50583.0   NA
[11,] 65982.3 66044.0 60014.1 52837.4 51277.8 53224.3   NA
[12,] 65527.5 57628.5 57685.4 55388.1 53596.4 53056.4   NA
> darr
          [,1]    [,2]    [,3]     [,4]     [,5]    [,6] [,7]
 [1,]  -2134.4 -6724.7 -7080.2  -8053.4  -8313.1 -4839.5   NA
 [2,]   3766.8  7021.5  5245.9   3343.6   2006.7  3800.2   NA
 [3,]   2853.3   189.4  5168.9   9982.8   8810.8  5482.1   NA
 [4,]  -7800.4 -6233.1 -8684.4 -10777.3 -13675.0 -9971.8   NA
 [5,]   5882.5  4380.2  2063.2    619.4   8509.2  5786.5   NA
 [6,] -10497.1 -6841.7 -2490.4  -1739.1  -3785.2 -5055.9   NA
 [7,]   2081.1   936.4 -1166.5  -1438.2  -6945.8 -3317.5   NA
 [8,]   5828.0  1439.5   419.4   3363.5   7993.2  8485.9   NA
 [9,]   -689.4  3332.7  4918.1   4625.9   1500.9 -4039.5   NA
[10,]   -165.3  2629.2 -1218.7  -7492.9  -5207.3  2641.3   NA
[11,]   -454.8 -8415.5 -2328.7   2550.7   2318.6  -167.9   NA
[12,]    387.1  5210.3  2717.7   4995.3    656.1      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/1gouj1272303989.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/2gouj1272303989.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/39xt41272303989.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/49xt41272303989.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]
[1,] 54252.50 49413.00 53213.2 58695.30 48723.5 54510.0 49454.1 46136.6
[2,] 60383.40 52070.30 54077.0 62887.80 49212.8 55518.2 53779.1 46991.0
[3,] 61620.95 54054.15 58348.9 65924.75 56193.9 58637.1 55499.3 54118.1
[4,] 65914.60 59189.90 66211.4 66400.80 60167.7 64547.9 57706.2 58642.6
[5,] 66857.20 64722.80 68489.6 71342.90 63542.5 69425.0 58927.9 61009.0
         [,9]    [,10]   [,11]   [,12]
[1,] 54622.50 50583.00 51277.8 53056.4
[2,] 54984.20 56485.10 52837.4 53596.4
[3,] 56009.55 60781.55 56619.2 56508.3
[4,] 60082.10 63414.80 65982.3 57685.4
[5,] 66837.00 66147.60 66044.0 57685.4

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 58053.15 49461.78 50521.82 63658.75 49127.64 52812.65 52966.19 46602.44
[2,] 65188.75 58646.52 66175.98 68190.75 63260.16 64461.55 58032.41 61633.76
         [,9]    [,10]    [,11]    [,12]
[1,] 52721.24 56311.67 48140.31 53870.76
[2,] 59297.86 65251.43 65098.09 59145.84

$out
[1] 65527.5

$group
[1] 12

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

Warning message:
In bxp(list(stats = c(54252.5, 60383.4, 61620.95, 65914.6, 66857.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/59xt41272303989.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,] -8313.10 2006.7  189.4 -13675.0  619.40 -10497.10 -6945.80  419.40 -4039.5
[2,] -8053.40 3343.6 2853.3 -10777.3 2063.20  -6841.70 -3317.50 1439.50  -689.4
[3,] -6902.45 3783.5 5325.5  -9328.1 5083.35  -4420.55 -1302.35 4595.75  2416.8
[4,] -4839.50 5245.9 8810.8  -7800.4 5882.50  -2490.40   936.40 7993.20  4625.9
[5,] -2134.40 7021.5 9982.8  -6233.1 8509.20  -1739.10  2081.10 8485.90  4918.1
       [,10]    [,11]  [,12]
[1,] -7492.9 -8415.50  387.1
[2,] -5207.3 -2328.70  656.1
[3,]  -692.0  -311.35 2717.7
[4,]  2629.2  2318.60 4995.3
[5,]  2641.3  2550.70 5210.3

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

$conf
          [,1]     [,2]    [,3]       [,4]     [,5]      [,6]      [,7]
[1,] -8975.519 2556.455 1482.72 -11248.297 2619.778 -7227.279 -4046.253
[2,] -4829.381 5010.545 9168.28  -7407.903 7546.922 -1613.821  1441.553
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,]  368.4018 -1011.740 -5746.796 -3309.009 -348.3678
[2,] 8823.0982  5845.340  4362.796  2686.309 5783.7678

$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(-8313.1, -8053.4, -6902.45, -4839.5, -2134.39999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6f3lj1272303989.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,] 58927.90 57628.50 55758.60 52340.90 46991.00 46136.60   NA
[2,] 64132.65 58916.25 56688.25 53308.25 51674.05 49433.55   NA
[3,] 66064.95 61791.25 58618.80 55453.15 54006.90 53134.80   NA
[4,] 67673.40 65979.30 61118.65 58017.35 57103.55 54381.25   NA
[5,] 71342.90 66400.80 66173.40 60403.10 62887.80 58695.30   NA

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6] [,7]
[1,] 64449.99 58569.75 56598.06 53305.3 51530.47 50878.12   NA
[2,] 67679.91 65012.75 60639.54 57601.0 56483.33 55391.48   NA

$out
[1] 65676.1

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(58927.9, 64132.65, 66064.95, 67673.4, 71342.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7f3lj1272303989.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,] 53502.57 54054.15 53142.12
[2,] 55628.22 55754.43 55458.96
[3,] 58160.40 56563.75 57932.22
[4,] 59997.03 59709.32 59926.61
[5,] 65196.05 61620.95 64964.41

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 56167.76 54759.89 55894.49
[2,] 60153.04 58367.61 59969.94

$out
[1] 65924.75

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(53502.5666666667, 55628.2166666667, 58160.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1gouj1272303989.ps tmp/1gouj1272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/2gouj1272303989.ps tmp/2gouj1272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/39xt41272303989.ps tmp/39xt41272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/49xt41272303989.ps tmp/49xt41272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/59xt41272303989.ps tmp/59xt41272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/6f3lj1272303989.ps tmp/6f3lj1272303989.png",intern=TRUE))
character(0)
> try(system("convert tmp/7f3lj1272303989.ps tmp/7f3lj1272303989.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.205   0.932   1.582