R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(14097.80,14776.80,16833.30,15385.50,15172.60,16858.90,14143.50,14731.80,16471.60,15214.00,17637.40,17972.40,16896.20,16698.00,19691.60,15930.70,17444.60,17699.40,15189.80,15672.70,17180.80,17664.90,17862.90,16162.30,17463.60,16772.10,19106.90,16721.30,18161.30,18509.90,17802.70,16409.90,17967.70,20286.60,19537.30,18021.90,20194.30,19049.60,20244.70,21473.30,19673.60,21053.20,20159.50,18203.60,21289.50,20432.30,17180.40,15816.80,15076.60,14531.60,15761.30,14345.50,13916.80,15496.80,14285.60,13597.30,16263.10,16773.30,15986.90,16842.60,15911.90,15782.90,18622.80,17422.50,16989.80,18990.50,16849.30,16511.30,18704.50,19111.10,19420.70,18985.10)
> 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,] 14097.8 16896.2 17463.6 20194.3 15076.6 15911.9   NA
 [2,] 14776.8 16698.0 16772.1 19049.6 14531.6 15782.9   NA
 [3,] 16833.3 19691.6 19106.9 20244.7 15761.3 18622.8   NA
 [4,] 15385.5 15930.7 16721.3 21473.3 14345.5 17422.5   NA
 [5,] 15172.6 17444.6 18161.3 19673.6 13916.8 16989.8   NA
 [6,] 16858.9 17699.4 18509.9 21053.2 15496.8 18990.5   NA
 [7,] 14143.5 15189.8 17802.7 20159.5 14285.6 16849.3   NA
 [8,] 14731.8 15672.7 16409.9 18203.6 13597.3 16511.3   NA
 [9,] 16471.6 17180.8 17967.7 21289.5 16263.1 18704.5   NA
[10,] 15214.0 17664.9 20286.6 20432.3 16773.3 19111.1   NA
[11,] 17637.4 17862.9 19537.3 17180.4 15986.9 19420.7   NA
[12,] 17972.4 16162.3 18021.9 15816.8 16842.6 18985.1   NA
> darr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6] [,7]
 [1,]   679.0  -198.2  -691.5 -1144.7  -545.0  -129.0   NA
 [2,]  2056.5  2993.6  2334.8  1195.1  1229.7  2839.9   NA
 [3,] -1447.8 -3760.9 -2385.6  1228.6 -1415.8 -1200.3   NA
 [4,]  -212.9  1513.9  1440.0 -1799.7  -428.7  -432.7   NA
 [5,]  1686.3   254.8   348.6  1379.6  1580.0  2000.7   NA
 [6,] -2715.4 -2509.6  -707.2  -893.7 -1211.2 -2141.2   NA
 [7,]   588.3   482.9 -1392.8 -1955.9  -688.3  -338.0   NA
 [8,]  1739.8  1508.1  1557.8  3085.9  2665.8  2193.2   NA
 [9,] -1257.6   484.1  2318.9  -857.2   510.2   406.6   NA
[10,]  2423.4   198.0  -749.3 -3251.9  -786.4   309.6   NA
[11,]   335.0 -1700.6 -1515.4 -1363.6   855.7  -435.6   NA
[12,] -1076.2  1301.3  2172.4  -740.2  -930.7      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/rcomp/tmp/1dkij1319651899.ps",horizontal=F,onefile=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/rcomp/tmp/2euk91319651899.ps",horizontal=F,onefile=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/rcomp/tmp/38w7v1319651899.ps",horizontal=F,onefile=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/rcomp/tmp/4aa731319651899.ps",horizontal=F,onefile=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,] 14097.80 14531.60 15761.30 14345.5 13916.8 15496.80 14143.50 13597.3
[2,] 15076.60 14776.80 16833.30 15385.5 15172.6 16858.90 14285.60 14731.8
[3,] 16404.05 16240.45 18864.85 16326.0 17217.2 18104.65 16019.55 16041.3
[4,] 17463.60 16772.10 19691.60 17422.5 18161.3 18990.50 17802.70 16511.3
[5,] 20194.30 19049.60 20244.70 17422.5 19673.6 21053.20 20159.50 18203.6
         [,9]   [,10]    [,11]   [,12]
[1,] 16263.10 15214.0 15986.90 15816.8
[2,] 16471.60 16773.3 17180.40 16162.3
[3,] 17574.25 18388.0 17750.15 17407.5
[4,] 18704.50 20286.6 19420.70 18021.9
[5,] 21289.50 20432.3 19537.30 18985.1

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]    [,6]     [,7]     [,8]
[1,] 14864.36 14953.42 17021.15 15012.07 15289.39 16729.7 13750.91 14893.47
[2,] 17943.74 17527.48 20708.55 17639.93 19145.01 19479.6 18288.19 17189.13
         [,9]    [,10]    [,11]    [,12]
[1,] 16133.96 16121.81 16305.08 16208.00
[2,] 19014.54 20654.19 19195.22 18607.00

$out
[1] 21473.3

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(14097.8, 15076.6, 16404.05, 17463.6, 20194.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5x4j81319651899.ps",horizontal=F,onefile=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,] -1144.7 1195.10 -3760.9 -1799.7  254.8 -2715.4 -1955.90 1508.1 -1257.60
[2,]  -691.5 1229.70 -2385.6  -432.7  348.6 -2509.6 -1392.80 1557.8  -857.20
[3,]  -371.6 2195.65 -1431.8  -320.8 1479.8 -1676.2  -513.15 1966.5   445.35
[4,]  -129.0 2839.90 -1200.3  1440.0 1686.3  -893.7   482.90 2665.8   510.20
[5,]   679.0 2993.60 -1200.3  1513.9 2000.7  -707.2   588.30 3085.9  2318.90
       [,10]   [,11]   [,12]
[1,] -786.40 -1700.6 -1076.2
[2,] -786.40 -1515.4  -930.7
[3,] -275.65  -899.6  -740.2
[4,]  309.60   335.0  1301.3
[5,]  309.60   855.7  2172.4

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

$conf
            [,1]     [,2]       [,3]      [,4]      [,5]       [,6]       [,7]
[1,] -734.430668 1157.019 -2196.3568 -1528.752  616.9403 -2718.5077 -1723.0371
[2,]   -8.769332 3234.281  -667.2432   887.152 2342.6597  -633.8923   696.7371
         [,8]      [,9]     [,10]      [,11]      [,12]
[1,] 1251.804 -436.6672 -982.6054 -2093.1678 -2317.3256
[2,] 2681.196 1327.3672  431.3054   293.9678   836.9256

$out
[1]  1228.6  2423.4 -3251.9

$group
[1]  3 10 10

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

Warning message:
In bxp(list(stats = c(-1144.7, -691.5, -371.6, -129, 679, 1195.10000000000,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6pnet1319651899.ps",horizontal=F,onefile=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,] 14097.80 15189.80 16409.90 15816.80 13597.30 15782.90   NA
[2,] 14754.30 16046.50 17117.85 18626.60 14315.55 16680.30   NA
[3,] 15299.75 17038.50 17994.80 20176.90 15286.70 18022.65   NA
[4,] 16846.10 17682.15 18808.40 20742.75 16125.00 18987.80   NA
[5,] 17972.40 19691.60 20286.60 21473.30 16842.60 19420.70   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 14345.67 16292.47 17223.73 19211.71 14461.40 16970.18   NA
[2,] 16253.83 17784.53 18765.87 21142.09 16112.00 19075.12   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(14097.8, 14754.3, 15299.75, 16846.1, 17972.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/7c8n81319651899.ps",horizontal=F,onefile=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,] 15854.43 16019.55 15726.49
[2,] 16505.90 16283.23 16170.79
[3,] 17096.65 17312.35 16987.73
[4,] 18040.49 17927.40 18066.32
[5,] 18376.77 18864.85 18494.46

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 16396.71 16562.43 16123.17
[2,] 17796.59 18062.27 17852.29

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15854.4333333333, 16505.9, 17096.65, 18040.4916666667,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1dkij1319651899.ps tmp/1dkij1319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/2euk91319651899.ps tmp/2euk91319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/38w7v1319651899.ps tmp/38w7v1319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/4aa731319651899.ps tmp/4aa731319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/5x4j81319651899.ps tmp/5x4j81319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/6pnet1319651899.ps tmp/6pnet1319651899.png",intern=TRUE))
character(0)
> try(system("convert tmp/7c8n81319651899.ps tmp/7c8n81319651899.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.088   0.508   2.579