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(41.84,42.94,49.14,44.61,40.22,44.23,45.85,53.38,53.26,51.80,55.30,57.81,63.96,63.77,59.15,56.12,57.42,63.52,61.71,63.01,68.18,72.03,69.75,74.41,74.33,64.24,60.03,59.44,62.50,55.04,58.34,61.92,67.65,67.68,70.30,75.26,71.44,76.36,81.71,92.60,90.60,92.23,94.09,102.79,109.65,124.05,132.69,135.81,116.07,101.42)
> 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] 50
> (np <- floor(n / par1))
[1] 4
> 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 4 4 4 4 4 4 4 4 4 4
> arr
       [,1]  [,2]  [,3]   [,4]   [,5]
 [1,] 41.84 63.96 74.33  71.44 116.07
 [2,] 42.94 63.77 64.24  76.36 101.42
 [3,] 49.14 59.15 60.03  81.71     NA
 [4,] 44.61 56.12 59.44  92.60     NA
 [5,] 40.22 57.42 62.50  90.60     NA
 [6,] 44.23 63.52 55.04  92.23     NA
 [7,] 45.85 61.71 58.34  94.09     NA
 [8,] 53.38 63.01 61.92 102.79     NA
 [9,] 53.26 68.18 67.65 109.65     NA
[10,] 51.80 72.03 67.68 124.05     NA
[11,] 55.30 69.75 70.30 132.69     NA
[12,] 57.81 74.41 75.26 135.81     NA
> darr
       [,1]  [,2]   [,3]   [,4]   [,5]
 [1,]  1.10 -0.19 -10.09   4.92 -14.65
 [2,]  6.20 -4.62  -4.21   5.35     NA
 [3,] -4.53 -3.03  -0.59  10.89     NA
 [4,] -4.39  1.30   3.06  -2.00     NA
 [5,]  4.01  6.10  -7.46   1.63     NA
 [6,]  1.62 -1.81   3.30   1.86     NA
 [7,]  7.53  1.30   3.58   8.70     NA
 [8,] -0.12  5.17   5.73   6.86     NA
 [9,] -1.46  3.85   0.03  14.40     NA
[10,]  3.50 -2.28   2.62   8.64     NA
[11,]  2.51  4.66   4.96   3.12     NA
[12,]  6.15 -0.08  -3.82 -19.74     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/161yz1229337508.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/27oqf1229337508.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/38taz1229337508.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/4w61x1229337508.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,] 63.96 63.77 49.140 44.610 40.22 44.230 45.850  53.380  53.260  51.800
[2,] 63.96 63.77 54.145 50.365 48.82 49.635 52.095  57.650  60.455  59.740
[3,] 71.44 64.24 59.590 57.780 59.96 59.280 60.025  62.465  67.915  69.855
[4,] 74.33 76.36 70.870 76.020 76.55 77.875 77.900  82.900  88.915  98.040
[5,] 74.33 76.36 81.710 92.600 90.60 92.230 94.090 102.790 109.650 124.050
       [,11]   [,12]
[1,]  55.300  57.810
[2,]  62.525  66.110
[3,]  70.025  74.835
[4,] 101.495 105.535
[5,] 132.690 135.810

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]    [,6]     [,7]    [,8]
[1,] 64.11258 55.34394 46.37725 37.51255 38.0533 36.9704 39.63905 42.5175
[2,] 78.76742 73.13606 72.80275 78.04745 81.8667 81.5896 80.41095 82.4125
        [,9]   [,10]    [,11]     [,12]
[1,] 45.4316  39.598  39.2387  43.68925
[2,] 90.3984 100.112 100.8113 105.98075

$out
[1]  41.84 116.07  42.94 101.42

$group
[1] 1 1 2 2

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

Warning message:
In bxp(list(stats = c(63.96, 63.96, 71.44, 74.33, 74.33, 63.77,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5ajog1229337508.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,] -14.65 -4.620 -4.53 -4.390 -7.460 -1.810 1.300 -0.120 -1.460 -2.28 2.510
[2,] -10.09 -4.415 -3.78 -3.195 -2.915 -0.095 2.440  2.525 -0.715  0.17 2.815
[3,]  -0.19  0.570 -1.81 -0.350  2.820  1.740 5.555  5.450  1.940  3.06 3.890
[4,]   1.10  5.775  5.15  2.180  5.055  2.580 8.115  6.295  9.125  6.07 4.810
[5,]   4.92  6.200 10.89  3.060  6.100  3.300 8.700  6.860 14.400  8.64 4.960
       [,12]
[1,] -19.740
[2,] -11.780
[3,]  -1.950
[4,]   3.035
[5,]   6.150

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

$conf
          [,1]    [,2]    [,3]     [,4]    [,5]     [,6]     [,7]   [,8]
[1,] -8.096826 -7.4801 -8.8647 -4.59625 -3.4763 -0.37325  1.07175 2.4717
[2,]  7.716826  8.6201  5.2447  3.89625  9.1163  3.85325 10.03825 8.4283
        [,9]  [,10]   [,11]     [,12]
[1,] -5.8336 -1.601 2.31395 -13.65385
[2,]  9.7136  7.721 5.46605   9.75385

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(-14.65, -10.09, -0.189999999999998, 1.09999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6oiq31229337508.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]
[1,] 40.220 56.120 55.040  71.440 101.420
[2,] 43.585 60.430 59.735  86.155 101.420
[3,] 47.495 63.645 63.370  93.345 108.745
[4,] 53.320 68.965 68.990 116.850 116.070
[5,] 57.810 74.410 75.260 135.810 116.070

$n
[1] 12 12 12 12  2

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]
[1,] 43.0548 59.75213 59.14873  79.3448  92.3776
[2,] 51.9352 67.53787 67.59127 107.3452 125.1124

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(40.22, 43.585, 47.495, 53.32, 57.81, 56.12,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7esie1229337508.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,] 62.50750 57.7800 60.48625
[2,] 63.47375 59.7750 61.42000
[3,] 70.01050 63.3525 67.75750
[4,] 76.78750 69.9400 72.83625
[5,] 85.82250 74.8350 80.32875

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 63.93801 58.71618 62.55047
[2,] 76.08299 67.98882 72.96453

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(62.5075, 63.47375, 70.0105, 76.7875, 85.8225,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/161yz1229337508.ps tmp/161yz1229337508.png")
> system("convert tmp/27oqf1229337508.ps tmp/27oqf1229337508.png")
> system("convert tmp/38taz1229337508.ps tmp/38taz1229337508.png")
> system("convert tmp/4w61x1229337508.ps tmp/4w61x1229337508.png")
> system("convert tmp/5ajog1229337508.ps tmp/5ajog1229337508.png")
> system("convert tmp/6oiq31229337508.ps tmp/6oiq31229337508.png")
> system("convert tmp/7esie1229337508.ps tmp/7esie1229337508.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.199   0.941   1.648