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(10519.2,10414.9,12476.8,12384.6,12266.7,12919.9,11497.3,12142,13919.4,12656.8,12034.1,13199.7,10881.3,11301.2,13643.9,12517,13981.1,14275.7,13435,13565.7,16216.3,12970,14079.9,14235,12213.4,12581,14130.4,14210.8,14378.5,13142.8,13714.7,13621.9,15379.8,13306.3,14391.2,14909.9,14025.4,12951.2,14344.3,16093.4,15413.6,14705.7,15972.8,16241.4,16626.4,17136.2,15622.9,18003.9,16136.1,14423.7,16789.4,16782.2,14133.8,12607,12004.5,12175.4,13268,12299.3,11800.6,13873.3,12315)
> 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] 61
> (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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]  [,6]
 [1,] 10519.2 10881.3 12213.4 14025.4 16136.1 12315
 [2,] 10414.9 11301.2 12581.0 12951.2 14423.7    NA
 [3,] 12476.8 13643.9 14130.4 14344.3 16789.4    NA
 [4,] 12384.6 12517.0 14210.8 16093.4 16782.2    NA
 [5,] 12266.7 13981.1 14378.5 15413.6 14133.8    NA
 [6,] 12919.9 14275.7 13142.8 14705.7 12607.0    NA
 [7,] 11497.3 13435.0 13714.7 15972.8 12004.5    NA
 [8,] 12142.0 13565.7 13621.9 16241.4 12175.4    NA
 [9,] 13919.4 16216.3 15379.8 16626.4 13268.0    NA
[10,] 12656.8 12970.0 13306.3 17136.2 12299.3    NA
[11,] 12034.1 14079.9 14391.2 15622.9 11800.6    NA
[12,] 13199.7 14235.0 14909.9 18003.9 13873.3    NA
> darr
         [,1]    [,2]    [,3]    [,4]    [,5] [,6]
 [1,]  -104.3   419.9   367.6 -1074.2 -1712.4   NA
 [2,]  2061.9  2342.7  1549.4  1393.1  2365.7   NA
 [3,]   -92.2 -1126.9    80.4  1749.1    -7.2   NA
 [4,]  -117.9  1464.1   167.7  -679.8 -2648.4   NA
 [5,]   653.2   294.6 -1235.7  -707.9 -1526.8   NA
 [6,] -1422.6  -840.7   571.9  1267.1  -602.5   NA
 [7,]   644.7   130.7   -92.8   268.6   170.9   NA
 [8,]  1777.4  2650.6  1757.9   385.0  1092.6   NA
 [9,] -1262.6 -3246.3 -2073.5   509.8  -968.7   NA
[10,]  -622.7  1109.9  1084.9 -1513.3  -498.7   NA
[11,]  1165.6   155.1   518.7  2381.0  2072.7   NA
[12,] -2318.4 -2021.6  -884.5 -1867.8 -1558.3   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/118io1258045872.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/2qkwd1258045872.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/35qll1258045872.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/4ba8d1258045872.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,] 10519.2 10414.9 13643.9 12384.6 13981.1 12607.0 11497.3 12142.0 13268.0
[2,] 10881.3 11301.2 13643.9 12517.0 13981.1 12919.9 12004.5 12175.4 13919.4
[3,] 12264.2 12581.0 14130.4 14210.8 14133.8 13142.8 13435.0 13565.7 15379.8
[4,] 14025.4 12951.2 14344.3 16093.4 14378.5 14275.7 13714.7 13621.9 16216.3
[5,] 16136.1 14423.7 14344.3 16782.2 14378.5 14705.7 15972.8 13621.9 16626.4
       [,10]   [,11]   [,12]
[1,] 12299.3 11800.6 13199.7
[2,] 12656.8 12034.1 13873.3
[3,] 12970.0 14079.9 14235.0
[4,] 13306.3 14391.2 14909.9
[5,] 13306.3 15622.9 14909.9

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

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 10236.15 11415.11 13635.5 11683.72 13853.00 12184.80 12226.58 12543.61
[2,] 14292.25 13746.89 14625.3 16737.88 14414.60 14100.80 14643.42 14587.79
         [,9]    [,10]    [,11]    [,12]
[1,] 13756.82 12511.06 12414.38 13502.54
[2,] 17002.78 13428.94 15745.42 14967.46

$out
[1] 12476.8 16789.4 12266.7 15413.6 16241.4 17136.2 18003.9

$group
[1]  3  3  5  5  8 10 12

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

Warning message:
In bxp(list(stats = c(10519.2, 10881.3, 12264.2, 14025.4, 16136.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5ahvp1258045872.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,] -1712.4 1393.1 -92.2 -679.8 -1526.8 -1422.6 130.7  385.0 -3246.3 -1513.3
[2,] -1074.2 1549.4 -92.2 -679.8 -1235.7  -840.7 130.7 1092.6 -2073.5  -622.7
[3,]  -104.3 2061.9  -7.2 -117.9  -707.9  -602.5 170.9 1757.9 -1262.6  -498.7
[4,]   367.6 2342.7  80.4  167.7   294.6   571.9 268.6 1777.4  -968.7  1084.9
[5,]   419.9 2365.7  80.4  167.7   653.2  1267.1 268.6 2650.6   509.8  1109.9
      [,11]   [,12]
[1,]  155.1 -2318.4
[2,]  518.7 -2021.6
[3,] 1165.6 -1867.8
[4,] 2072.7 -1558.3
[5,] 2381.0  -884.5

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

$conf
           [,1]     [,2]      [,3]      [,4]       [,5]       [,6]      [,7]
[1,] -1123.0722 1501.356 -129.1587 -716.7414 -1789.2061 -1600.6396  73.46021
[2,]   914.4722 2622.444  114.7587  480.9414   373.4061   395.6396 268.33979
         [,8]       [,9]      [,10]      [,11]     [,12]
[1,] 1274.022 -2043.2489 -1705.2859   67.54751 -2195.167
[2,] 2241.778  -481.9511   707.8859 2263.65249 -1540.433

$out
[1] -1126.9  1749.1  1464.1 -2648.4   644.7   -92.8

$group
[1] 3 3 4 4 7 7

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

Warning message:
In bxp(list(stats = c(-1712.4, -1074.20000000000, -104.300000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/65naf1258045872.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,] 10414.90 10881.30 12213.40 12951.20 11800.60 12315
[2,] 11765.70 12743.50 13224.55 14525.00 12237.35 12315
[3,] 12325.65 13604.80 13922.55 15797.85 13570.65 12315
[4,] 12788.35 14157.45 14384.85 16433.90 15279.90 12315
[5,] 13919.40 16216.30 15379.80 18003.90 16789.40 12315

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]  [,6]
[1,] 11859.21 12959.89 13393.33 14927.19 12182.92 12315
[2,] 12792.09 14249.71 14451.77 16668.51 14958.38 12315

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(10414.9, 11765.7, 12325.65, 12788.35, 13919.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7k04k1258045872.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,] 12334.40 12264.2 12126.20
[2,] 13427.54 13056.4 12879.12
[3,] 13629.73 13822.8 13405.23
[4,] 14337.28 14172.3 14242.50
[5,] 15081.98 15379.8 15067.85

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 13214.79 13313.83 12783.38
[2,] 14044.67 14331.77 14027.07

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(12334.4, 13427.54, 13629.73, 14337.28, 15081.98,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/118io1258045872.ps tmp/118io1258045872.png")
> system("convert tmp/2qkwd1258045872.ps tmp/2qkwd1258045872.png")
> system("convert tmp/35qll1258045872.ps tmp/35qll1258045872.png")
> system("convert tmp/4ba8d1258045872.ps tmp/4ba8d1258045872.png")
> system("convert tmp/5ahvp1258045872.ps tmp/5ahvp1258045872.png")
> system("convert tmp/65naf1258045872.ps tmp/65naf1258045872.png")
> system("convert tmp/7k04k1258045872.ps tmp/7k04k1258045872.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.894   1.502   2.208