R version 2.7.0 (2008-04-22)
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(2236,2084.9,2409.5,2199.3,2203.5,2254.1,1975.8,1742.2,2520.6,2438.1,2126.3,2267.5,2201.1,2128.5,2596,2458.2,2210.5,2621.2,2231.4,2103.6,2685.8,2539.3,2462.4,2693.3,2307.7,2385.9,2737.6,2653.9,2545.4,2848.8,2359.5,2488.3,2861.1,2717.9,2844,2749,2652.9,2660.2,3187.1,2774.1,3158.2,3244.6,2665.5,2820.8,2983.4,3077.4,3024.8,2731.8,3046.2,2834.8,3292.8,2946.1,3196.9,3284.2,3003,2979,3137.4,3630.2,3270.7,2942.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] 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,] 2236.0 2201.1 2307.7 2652.9 3046.2   NA
 [2,] 2084.9 2128.5 2385.9 2660.2 2834.8   NA
 [3,] 2409.5 2596.0 2737.6 3187.1 3292.8   NA
 [4,] 2199.3 2458.2 2653.9 2774.1 2946.1   NA
 [5,] 2203.5 2210.5 2545.4 3158.2 3196.9   NA
 [6,] 2254.1 2621.2 2848.8 3244.6 3284.2   NA
 [7,] 1975.8 2231.4 2359.5 2665.5 3003.0   NA
 [8,] 1742.2 2103.6 2488.3 2820.8 2979.0   NA
 [9,] 2520.6 2685.8 2861.1 2983.4 3137.4   NA
[10,] 2438.1 2539.3 2717.9 3077.4 3630.2   NA
[11,] 2126.3 2462.4 2844.0 3024.8 3270.7   NA
[12,] 2267.5 2693.3 2749.0 2731.8 2942.3   NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5] [,6]
 [1,] -151.1  -72.6   78.2    7.3 -211.4   NA
 [2,]  324.6  467.5  351.7  526.9  458.0   NA
 [3,] -210.2 -137.8  -83.7 -413.0 -346.7   NA
 [4,]    4.2 -247.7 -108.5  384.1  250.8   NA
 [5,]   50.6  410.7  303.4   86.4   87.3   NA
 [6,] -278.3 -389.8 -489.3 -579.1 -281.2   NA
 [7,] -233.6 -127.8  128.8  155.3  -24.0   NA
 [8,]  778.4  582.2  372.8  162.6  158.4   NA
 [9,]  -82.5 -146.5 -143.2   94.0  492.8   NA
[10,] -311.8  -76.9  126.1  -52.6 -359.5   NA
[11,]  141.2  230.9  -95.0 -293.0 -328.4   NA
[12,]  -66.4 -385.6  -96.1  314.4     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/14i2c1229091863.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/2y9uy1229091863.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/36w7b1229091863.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/4a9mq1229091863.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]  [,10]
[1,] 2201.1 2084.9 2409.5 2199.3 2203.5 2254.1 1975.8 1742.2 2520.6 2438.1
[2,] 2236.0 2128.5 2596.0 2458.2 2210.5 2621.2 2231.4 2103.6 2685.8 2539.3
[3,] 2307.7 2385.9 2737.6 2653.9 2545.4 2848.8 2359.5 2488.3 2861.1 2717.9
[4,] 2652.9 2660.2 3187.1 2774.1 3158.2 3244.6 2665.5 2820.8 2983.4 3077.4
[5,] 3046.2 2834.8 3292.8 2946.1 3196.9 3284.2 3003.0 2979.0 3137.4 3630.2
      [,11]  [,12]
[1,] 2126.3 2693.3
[2,] 2462.4 2693.3
[3,] 2844.0 2731.8
[4,] 3024.8 2749.0
[5,] 3270.7 2749.0

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

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2013.120 2010.202 2319.93 2430.686 1875.758 2408.307 2052.766 1981.528
[2,] 2602.280 2761.598 3155.27 2877.114 3215.042 3289.293 2666.234 2995.072
         [,9]   [,10]    [,11]    [,12]
[1,] 2650.817 2337.68 2446.610 2692.443
[2,] 3071.383 3098.12 3241.390 2771.157

$out
[1] 2267.5 2942.3

$group
[1] 12 12

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

Warning message:
In bxp(list(stats = c(2201.1, 2236, 2307.7, 2652.9, 3046.2, 2084.9,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5pog61229091863.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]  [,11]
[1,] -211.4 324.6 -413.0 -247.7  50.6 -579.1 -233.6 158.4 -146.5 -359.5 -328.4
[2,] -151.1 351.7 -346.7 -108.5  86.4 -489.3 -127.8 162.6 -143.2 -311.8 -293.0
[3,]  -72.6 458.0 -210.2    4.2  87.3 -389.8  -24.0 372.8  -82.5  -76.9  -95.0
[4,]    7.3 467.5 -137.8  250.8 303.4 -281.2  128.8 582.2   94.0  -52.6  141.2
[5,]   78.2 526.9  -83.7  384.1 410.7 -278.3  155.3 778.4   94.0  126.1  230.9
       [,12]
[1,] -385.60
[2,] -240.85
[3,]  -81.25
[4,]  124.00
[5,]  314.40

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

$conf
           [,1]    [,2]       [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -184.52504 376.176 -357.80821 -249.6805 -66.03165 -536.8429 -205.3129
[2,]   39.32504 539.824  -62.59179  258.0805 240.63165 -242.7571  157.3129
         [,8]       [,9]     [,10]     [,11]     [,12]
[1,]  76.3117 -250.10492 -260.0501 -401.8046 -369.4815
[2,] 669.2883   85.10492  106.2501  211.8046  206.9815

$out
[1] 492.8

$group
[1] 9

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

Warning message:
In bxp(list(stats = c(-211.400000000000, -151.1, -72.5999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6936d1229091863.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,] 1975.80 2103.6 2307.7 2652.90 2834.80   NA
[2,] 2105.60 2205.8 2437.1 2698.65 2962.55   NA
[3,] 2219.75 2460.3 2685.9 2902.10 3091.80   NA
[4,] 2338.50 2608.6 2796.5 3117.80 3277.45   NA
[5,] 2520.60 2693.3 2861.1 3244.60 3630.20   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]    [,2]     [,3]     [,4]     [,5] [,6]
[1,] 2113.523 2276.58 2521.975 2710.923 2948.172   NA
[2,] 2325.977 2644.02 2849.825 3093.277 3235.428   NA

$out
[1] 1742.2

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(1975.8, 2105.6, 2219.75, 2338.5, 2520.6, 2103.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7wut81229091863.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,] 2418.86 2307.7 2394.350
[2,] 2467.91 2437.1 2455.325
[3,] 2669.84 2685.9 2702.750
[4,] 2841.13 2790.8 2821.475
[5,] 2880.58 2861.1 2932.900

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2499.612 2524.575 2535.747
[2,] 2840.068 2847.225 2869.753

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2418.86, 2467.91, 2669.84, 2841.13, 2880.58,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/14i2c1229091863.ps tmp/14i2c1229091863.png")
> system("convert tmp/2y9uy1229091863.ps tmp/2y9uy1229091863.png")
> system("convert tmp/36w7b1229091863.ps tmp/36w7b1229091863.png")
> system("convert tmp/4a9mq1229091863.ps tmp/4a9mq1229091863.png")
> system("convert tmp/5pog61229091863.ps tmp/5pog61229091863.png")
> system("convert tmp/6936d1229091863.ps tmp/6936d1229091863.png")
> system("convert tmp/7wut81229091863.ps tmp/7wut81229091863.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.576   1.602   2.837