R version 2.7.2 (2008-08-25)
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(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,17422,16704.5,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,20769.2249)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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 5 5
> arr
         [,1]    [,2]    [,3]    [,4]     [,5] [,6]
 [1,] 15370.6 16379.4 17838.7 17170.0 18594.60   NA
 [2,] 14956.9 15420.3 17209.4 19439.5 19823.10   NA
 [3,] 15469.7 17870.5 18586.5 19795.4 20844.40   NA
 [4,] 15101.8 15912.8 16258.1 17574.9 19640.20   NA
 [5,] 11703.7 13866.5 15141.6 16165.4 17735.40   NA
 [6,] 16283.6 17823.2 19202.1 19464.6 19813.60   NA
 [7,] 16726.5 17872.0 17746.5 19932.1 22238.50   NA
 [8,] 14968.9 17422.0 19090.1 19961.2 20682.20   NA
 [9,] 14861.0 16704.5 18040.3 17343.4 17818.60   NA
[10,] 14583.3 15991.2 17515.5 18924.2 21872.10   NA
[11,] 15305.8 16583.6 17751.8 18574.1 22117.00   NA
[12,] 17903.9 19123.5 21072.4 21350.6 20769.22   NA
> darr
         [,1]    [,2]    [,3]    [,4]      [,5] [,6]
 [1,]  -413.7  -959.1  -629.3  2269.5  1228.500   NA
 [2,]   512.8  2450.2  1377.1   355.9  1021.300   NA
 [3,]  -367.9 -1957.7 -2328.4 -2220.5 -1204.200   NA
 [4,] -3398.1 -2046.3 -1116.5 -1409.5 -1904.800   NA
 [5,]  4579.9  3956.7  4060.5  3299.2  2078.200   NA
 [6,]   442.9    48.8 -1455.6   467.5  2424.900   NA
 [7,] -1757.6  -450.0  1343.6    29.1 -1556.300   NA
 [8,]  -107.9  -717.5 -1049.8 -2617.8 -2863.600   NA
 [9,]  -277.7  -713.3  -524.8  1580.8  4053.500   NA
[10,]   722.5   592.4   236.3  -350.1   244.900   NA
[11,]  2598.1  2539.9  3320.6  2776.5 -1347.775   NA
[12,] -1524.5 -1284.8 -3902.4 -2756.0        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/freestat/rcomp/tmp/19xv41225554476.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/21jia1225554476.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/3ti2a1225554476.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/4n1jf1225554476.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]
[1,] 15370.6 14956.9 15469.7 15101.8 11703.7 16283.6 16726.5 14968.9 16704.5
[2,] 16379.4 15420.3 17870.5 15912.8 13866.5 17823.2 17746.5 17422.0 16704.5
[3,] 17170.0 17209.4 18586.5 16258.1 15141.6 19202.1 17872.0 19090.1 17343.4
[4,] 17838.7 19439.5 19795.4 17574.9 16165.4 19464.6 19932.1 19961.2 17818.6
[5,] 18594.6 19823.1 20844.4 19640.2 17735.4 19813.6 22238.5 20682.2 18040.3
       [,10]   [,11]    [,12]
[1,] 14583.3 15305.8 17903.90
[2,] 15991.2 16583.6 19123.50
[3,] 17515.5 17751.8 20769.22
[4,] 18924.2 18574.1 21072.40
[5,] 21872.1 18574.1 21350.60

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 16138.86 14369.44 17226.37 15083.66 13517.20 18042.29 16327.66 17295.91
[2,] 18201.14 20049.36 19946.63 17432.54 16766.00 20361.91 19416.34 20884.29
         [,9]    [,10]    [,11]    [,12]
[1,] 16556.18 15443.05 16345.32 19392.14
[2,] 18130.62 19587.95 19158.28 22146.31

$out
[1] 14861 22117

$group
[1]  9 11

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

Warning message:
In bxp(list(stats = c(15370.6, 16379.4, 17170, 17838.7, 18594.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5r9we1225554476.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]
[1,] -959.1  355.9 -2328.4 -2046.3 3299.2  48.8 -1757.6 -2863.6 -713.3 236.3
[2,] -629.3  512.8 -2220.5 -2046.3 3299.2  48.8 -1556.3 -2617.8 -524.8 236.3
[3,] -413.7 1021.3 -1957.7 -1904.8 3956.7 442.9  -450.0 -1049.8 -277.7 244.9
[4,] 1228.5 1377.1 -1204.2 -1409.5 4060.5 467.5    29.1  -717.5 1580.8 592.4
[5,] 2269.5 2450.2  -367.9 -1116.5 4579.9 467.5  1343.6  -107.9 4053.5 722.5
      [,11]    [,12]
[1,] 2539.9 -3902.40
[2,] 2539.9 -3329.20
[3,] 2598.1 -2140.25
[4,] 2776.5 -1404.65
[5,] 2776.5 -1284.80

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

$conf
           [,1]      [,2]      [,3]      [,4]     [,5]     [,6]       [,7]
[1,] -1726.4168  410.5878 -2675.815 -2354.761 3418.767 147.0476 -1570.2396
[2,]   899.0168 1632.0122 -1239.585 -1454.839 4494.633 738.7524   670.2396
           [,8]      [,9]      [,10]    [,11]      [,12]
[1,] -2392.5472 -1765.512  -6.719363 2430.919 -3660.6445
[2,]   292.9472  1210.112 496.519363 2765.281  -619.8555

$out
[1] -3398.100  2078.200 -1455.600  2424.900  -350.100  3320.600 -1347.775

$group
[1]  4  5  6  6 10 11 11

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

Warning message:
In bxp(list(stats = c(-959.1, -629.299999999999, -413.700000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6i9cr1225554476.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]
[1,] 14583.30 13866.50 16258.10 16165.40 17735.40   NA
[2,] 14908.95 15952.00 17362.45 17459.15 19117.40   NA
[3,] 15203.80 16644.05 17795.25 19181.85 20252.65   NA
[4,] 15876.65 17846.85 18838.30 19863.75 21358.25   NA
[5,] 16726.50 19123.50 19202.10 21350.60 22238.50   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 14762.43 15779.80 17122.10 18085.10 19230.58   NA
[2,] 15645.17 17508.30 18468.40 20278.60 21274.72   NA

$out
[1] 11703.7 17903.9 15141.6 21072.4

$group
[1] 1 1 3 3

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

Warning message:
In bxp(list(stats = c(14583.3, 14908.95, 15203.8, 15876.65, 16726.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/706qr1225554476.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,] 14922.52 15141.60 15015.95
[2,] 17012.11 17189.70 17185.30
[3,] 17921.86 17633.65 17518.28
[4,] 18515.36 18838.30 18762.28
[5,] 20043.92 20769.22 20097.95

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 17236.22 16881.71 16799.01
[2,] 18607.50 18385.59 18237.54

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(14922.52, 17012.11, 17921.86, 18515.36, 20043.92498,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/19xv41225554476.ps tmp/19xv41225554476.png")
> system("convert tmp/21jia1225554476.ps tmp/21jia1225554476.png")
> system("convert tmp/3ti2a1225554476.ps tmp/3ti2a1225554476.png")
> system("convert tmp/4n1jf1225554476.ps tmp/4n1jf1225554476.png")
> system("convert tmp/5r9we1225554476.ps tmp/5r9we1225554476.png")
> system("convert tmp/6i9cr1225554476.ps tmp/6i9cr1225554476.png")
> system("convert tmp/706qr1225554476.ps tmp/706qr1225554476.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.142   1.450   2.393