R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(2529,2196,3202,2718,2728,2354,2697,2651,2067,2641,2539,2294,2712,2314,3092,2677,2813,2668,2939,2617,2231,2481,2421,2408,2560,2100,3315,2801,2403,3024,2507,2980,2211,2471,2594,2452,2232,2373,3127,2802,2641,2787,2619,2806,2193,2323,2529,2412,2262,2154,3230,2295,2715,2733,2317,2730,1913,2390,2484,1960)
> 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,] 2529 2712 2560 2232 2262   NA
 [2,] 2196 2314 2100 2373 2154   NA
 [3,] 3202 3092 3315 3127 3230   NA
 [4,] 2718 2677 2801 2802 2295   NA
 [5,] 2728 2813 2403 2641 2715   NA
 [6,] 2354 2668 3024 2787 2733   NA
 [7,] 2697 2939 2507 2619 2317   NA
 [8,] 2651 2617 2980 2806 2730   NA
 [9,] 2067 2231 2211 2193 1913   NA
[10,] 2641 2481 2471 2323 2390   NA
[11,] 2539 2421 2594 2529 2484   NA
[12,] 2294 2408 2452 2412 1960   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,] -333 -398 -460  141 -108   NA
 [2,] 1006  778 1215  754 1076   NA
 [3,] -484 -415 -514 -325 -935   NA
 [4,]   10  136 -398 -161  420   NA
 [5,] -374 -145  621  146   18   NA
 [6,]  343  271 -517 -168 -416   NA
 [7,]  -46 -322  473  187  413   NA
 [8,] -584 -386 -769 -613 -817   NA
 [9,]  574  250  260  130  477   NA
[10,] -102  -60  123  206   94   NA
[11,] -245  -13 -142 -117 -524   NA
[12,]  418  152 -220 -150   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/1cge91262093117.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/2al601262093117.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/3m56t1262093117.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/4fw0w1262093117.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] [,11] [,12]
[1,] 2232 2100 3092 2677 2641 2668 2317 2617 1913  2323  2421  2294
[2,] 2262 2154 3127 2677 2641 2668 2507 2651 2067  2390  2484  2294
[3,] 2529 2196 3202 2718 2715 2733 2619 2730 2193  2471  2529  2408
[4,] 2560 2314 3230 2801 2728 2787 2697 2806 2211  2481  2539  2412
[5,] 2712 2373 3315 2802 2813 2787 2939 2980 2231  2481  2594  2452

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2318.434 2082.944 3129.220 2630.382 2653.526 2648.915 2484.746 2620.477
[2,] 2739.566 2309.056 3274.780 2805.618 2776.474 2817.085 2753.254 2839.523
        [,9]    [,10]    [,11]    [,12]
[1,] 2091.25 2406.700 2490.137 2324.621
[2,] 2294.75 2535.300 2567.863 2491.379

$out
[1] 2295 2403 2354 3024 2641 1960

$group
[1]  4  5  6  6 10 12

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

Warning message:
In bxp(list(stats = c(2232, 2262, 2529, 2560, 2712, 2100, 2154,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5r21i1262093117.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] [,12]
[1,] -460  754 -514 -398 -374 -517 -322 -817  130  -102  -245  -220
[2,] -398  778 -514 -161 -145 -416  -46 -769  250   -60  -245  -185
[3,] -333 1006 -484   10   18 -168  187 -613  260    94  -142     1
[4,] -108 1076 -415  136  146  271  413 -584  477   123  -117   285
[5,]  141 1215 -325  420  146  343  473 -386  574   206   -13   418

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

$conf
          [,1]     [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] -537.9133  795.434 -553.9532 -199.8595 -187.6199 -653.4325 -137.3282
[2,] -128.0867 1216.566 -414.0468  219.8595  223.6199  317.4325  511.3282
          [,8]      [,9]     [,10]      [,11]  [,12]
[1,] -743.7205  99.60237 -35.30734 -232.44448 -370.3
[2,] -482.2795 420.39763 223.30734  -51.55552  372.3

$out
[1] -935  621 -524

$group
[1]  3  5 11

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

Warning message:
In bxp(list(stats = c(-460, -398, -333, -108, 141, 754, 778, 1006,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6i12u1262093117.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,] 2067.0 2231.0 2100.0 2193.0 1913.0   NA
[2,] 2324.0 2414.5 2427.5 2348.0 2208.0   NA
[3,] 2590.0 2642.5 2533.5 2574.0 2353.5   NA
[4,] 2707.5 2762.5 2890.5 2794.5 2722.5   NA
[5,] 3202.0 3092.0 3315.0 3127.0 3230.0   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 2415.083 2483.775 2322.323 2370.348 2118.833   NA
[2,] 2764.917 2801.225 2744.677 2777.652 2588.167   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2067, 2324, 2590, 2707.5, 3202, 2231, 2414.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7pk6v1262093117.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,] 2123.0 2193.0 2139.00
[2,] 2382.1 2439.5 2382.00
[3,] 2564.6 2574.0 2556.75
[4,] 2686.6 2724.0 2728.00
[5,] 2756.8 2733.0 3178.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 2425.716 2444.238 2398.937
[2,] 2703.484 2703.762 2714.563

$out
[1] 3193.2 3202.0

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(2123, 2382.1, 2564.6, 2686.6, 2756.8, 2193,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> try(system("convert tmp/1cge91262093117.ps tmp/1cge91262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/2al601262093117.ps tmp/2al601262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/3m56t1262093117.ps tmp/3m56t1262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/4fw0w1262093117.ps tmp/4fw0w1262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/5r21i1262093117.ps tmp/5r21i1262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/6i12u1262093117.ps tmp/6i12u1262093117.png",intern=TRUE))
character(0)
> try(system("convert tmp/7pk6v1262093117.ps tmp/7pk6v1262093117.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.208   0.958   2.434