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(15916.4,16535.9,15796,14418.6,15044.5,14944.2,16754.8,14254,15454.9,15644.8,14568.3,12520.2,14803,15873.2,14755.3,12875.1,14291.1,14205.3,15859.4,15258.9,15498.6,15106.5,15023.6,12083,15761.3,16943,15070.3,13659.6,14768.9,14725.1,15998.1,15370.6,14956.9,15469.7,15101.8,11703.7,16283.6,16726.5,14968.9,14861,14583.3,15305.8,17903.9,16379.4,15420.3,17870.5,15912.8,13866.5,17823.2,17872,17420.4,16704.4,15991.2,16583.6,19123.5,17838.7,17209.4,18586.5,16258.1,15141.6,19202.1,17746.5,19090.1,18040.3,17515.5,17751.8,21072.4,17170,19439.5,19795.4,17574.9,16165.4,19464.6,19932.1,19961.2,17343.4,18924.2,18574.1,21350.6,18594.6,19823.1,20844.4,19640.2,17735.4,19813.6,22238.5,20682.2,17818.6,21872.1,22117,21865.9,23451.3,20953.7,22497.3)
> 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] 94
> (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 8 8 8 8 8 8 8 8 8 7 7
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
 [1,] 15916.4 14803.0 15761.3 16283.6 17823.2 19202.1 19464.6 19813.6
 [2,] 16535.9 15873.2 16943.0 16726.5 17872.0 17746.5 19932.1 22238.5
 [3,] 15796.0 14755.3 15070.3 14968.9 17420.4 19090.1 19961.2 20682.2
 [4,] 14418.6 12875.1 13659.6 14861.0 16704.4 18040.3 17343.4 17818.6
 [5,] 15044.5 14291.1 14768.9 14583.3 15991.2 17515.5 18924.2 21872.1
 [6,] 14944.2 14205.3 14725.1 15305.8 16583.6 17751.8 18574.1 22117.0
 [7,] 16754.8 15859.4 15998.1 17903.9 19123.5 21072.4 21350.6 21865.9
 [8,] 14254.0 15258.9 15370.6 16379.4 17838.7 17170.0 18594.6 23451.3
 [9,] 15454.9 15498.6 14956.9 15420.3 17209.4 19439.5 19823.1 20953.7
[10,] 15644.8 15106.5 15469.7 17870.5 18586.5 19795.4 20844.4 22497.3
[11,] 14568.3 15023.6 15101.8 15912.8 16258.1 17574.9 19640.2      NA
[12,] 12520.2 12083.0 11703.7 13866.5 15141.6 16165.4 17735.4      NA
> darr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
 [1,]   619.5  1070.2  1181.7   442.9    48.8 -1455.6   467.5  2424.9
 [2,]  -739.9 -1117.9 -1872.7 -1757.6  -451.6  1343.6    29.1 -1556.3
 [3,] -1377.4 -1880.2 -1410.7  -107.9  -716.0 -1049.8 -2617.8 -2863.6
 [4,]   625.9  1416.0  1109.3  -277.7  -713.2  -524.8  1580.8  4053.5
 [5,]  -100.3   -85.8   -43.8   722.5   592.4   236.3  -350.1   244.9
 [6,]  1810.6  1654.1  1273.0  2598.1  2539.9  3320.6  2776.5  -251.1
 [7,] -2500.8  -600.5  -627.5 -1524.5 -1284.8 -3902.4 -2756.0  1585.4
 [8,]  1200.9   239.7  -413.7  -959.1  -629.3  2269.5  1228.5 -2497.6
 [9,]   189.9  -392.1   512.8  2450.2  1377.1   355.9  1021.3  1543.6
[10,] -1076.5   -82.9  -367.9 -1957.7 -2328.4 -2220.5 -1204.2      NA
[11,] -2048.1 -2940.6 -3398.1 -2046.3 -1116.5 -1409.5 -1904.8      NA
[12,]  2282.8  3678.3  4579.9  3956.7  4060.5  3299.2  2078.2      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/190yd1229889532.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/2o61v1229889532.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/3bgmf1229889532.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/4ypki1229889532.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,] 14803.00 15873.20 14755.30 12875.1 14291.10 14205.30 15859.40 14254.00
[2,] 15838.85 16631.20 15019.60 14039.1 14676.10 14834.65 16376.45 15314.75
[3,] 17053.40 17344.75 16608.20 15782.7 15517.85 15944.70 18513.70 16774.70
[4,] 19333.35 18902.05 19525.65 17581.0 18219.85 18162.95 21211.50 18216.65
[5,] 19813.60 22238.50 20682.20 18040.3 21872.10 22117.00 21865.90 18594.60
        [,9]    [,10]   [,11]   [,12]
[1,] 14956.9 15106.50 14568.3 11703.7
[2,] 15437.6 15557.25 15062.7 12301.6
[3,] 16354.0 18228.50 15912.8 13866.5
[4,] 19631.3 20319.90 16916.5 15653.5
[5,] 20953.7 22497.30 19640.2 17735.4

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 15101.32 16076.22 14091.06 13804.14 13538.26 14085.46 15812.77 15153.66
[2,] 19005.48 18613.28 19125.34 17761.26 17497.44 17803.94 21214.63 18395.74
         [,9]    [,10]    [,11]   [,12]
[1,] 14011.34 15568.02 14805.74 11864.8
[2,] 18696.66 20888.98 17019.86 15868.2

$out
[1] 23451.3

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(14803, 15838.85, 17053.4, 19333.35, 19813.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5a4ox1229889532.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,]   48.80 -1872.70 -2863.60 -713.20 -350.10 -251.10 -3902.40 -2497.6 -392.10
[2,]  245.85 -1656.95 -2249.00 -401.25  -93.05 1463.55 -2628.40  -794.2  272.90
[3,]  543.50  -928.90 -1394.05  867.60   96.25 2175.25 -1404.65   -87.0  767.05
[4,] 1125.95  -211.25  -882.90 1498.40  418.65 2687.30  -614.00  1214.7 1460.35
[5,] 2424.90  1343.60  -107.90 4053.50  722.50 3320.60  1585.40  2269.5 2450.20
       [,10]    [,11]  [,12]
[1,] -2328.4 -3398.10 2078.2
[2,] -2089.1 -2494.35 2791.0
[3,] -1204.2 -2046.30 3678.3
[4,]  -722.2 -1657.15 4008.6
[5,]   -82.9 -1116.50 4579.9

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

$conf
          [,1]       [,2]       [,3]      [,4]      [,5]     [,6]       [,7]
[1,]   51.8635 -1736.4888 -2157.1731 -193.5718 -189.5930 1491.646 -2529.9228
[2,] 1035.1365  -121.3112  -630.9269 1928.7718  382.0930 2858.854  -279.3772
          [,8]      [,9]      [,10]     [,11]    [,12]
[1,] -1209.200  103.7234 -2020.4906 -2546.262 2951.169
[2,]  1035.200 1430.3766  -387.9094 -1546.338 4405.431

$out
[1] -1455.6

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(48.7999999999993, 245.849999999999, 543.500000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/63c5y1229889532.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,] 12520.20 12875.10 13659.60 13866.50 15141.60 16165.40 17343.40 19813.6
[2,] 14493.45 14248.20 14747.00 14914.95 16420.85 17545.20 18584.35 20682.2
[3,] 15249.70 14913.30 15086.05 15666.55 17314.90 17896.05 19552.40 21869.0
[4,] 15856.20 15378.75 15615.50 16552.95 17855.35 19320.80 19946.65 22238.5
[5,] 16754.80 15873.20 15998.10 17903.90 19123.50 21072.40 21350.60 23451.3

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 14628.14 14397.65 14689.92 14919.45 16660.61 17086.19 18931.05 21091.41
[2,] 15871.26 15428.95 15482.18 16413.65 17969.19 18705.91 20173.75 22646.59

$out
[1] 12083.0 16943.0 11703.7 17818.6

$group
[1] 2 3 3 8

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

Warning message:
In bxp(list(stats = c(12520.2, 14493.45, 15249.7, 15856.2, 16754.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7a7fr1229889532.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,] 15715.12 13866.50 15845.53
[2,] 16460.47 15847.75 16142.34
[3,] 17253.87 16481.10 16930.79
[4,] 17683.47 17199.08 17552.81
[5,] 18741.08 18513.70 18853.79

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 16696.05 15864.75 16287.47
[2,] 17811.68 17097.45 17574.12

$out
[1] 14173.69 13977.55

$group
[1] 1 3

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

Warning message:
In bxp(list(stats = c(15715.125, 16460.475, 17253.86875, 17683.46875,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/190yd1229889532.ps tmp/190yd1229889532.png")
> system("convert tmp/2o61v1229889532.ps tmp/2o61v1229889532.png")
> system("convert tmp/3bgmf1229889532.ps tmp/3bgmf1229889532.png")
> system("convert tmp/4ypki1229889532.ps tmp/4ypki1229889532.png")
> system("convert tmp/5a4ox1229889532.ps tmp/5a4ox1229889532.png")
> system("convert tmp/63c5y1229889532.ps tmp/63c5y1229889532.png")
> system("convert tmp/7a7fr1229889532.ps tmp/7a7fr1229889532.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.887   1.485   2.186