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(101.73
+ ,101.63
+ ,101.43
+ ,101.34
+ ,101.01
+ ,100.89
+ ,100.93
+ ,100.77
+ ,100.3
+ ,99.86
+ ,99.71
+ ,99.93
+ ,99.88
+ ,99.92
+ ,99.87
+ ,99.63
+ ,100.05
+ ,99.88
+ ,100.11
+ ,100.05
+ ,100.07
+ ,100.2
+ ,100.21
+ ,99.76
+ ,99.41
+ ,99.24
+ ,99.65
+ ,99.7
+ ,99.79
+ ,99.84
+ ,101
+ ,101.62
+ ,101.98
+ ,101.46
+ ,102.28
+ ,102.14
+ ,102.02
+ ,102.21
+ ,101.61
+ ,102.38
+ ,102.19
+ ,102.04
+ ,101.76
+ ,101.9
+ ,102.01
+ ,102.37
+ ,103.04
+ ,103.42
+ ,103.76
+ ,104.41
+ ,104.75
+ ,104.28
+ ,103.89
+ ,104.09
+ ,103.8
+ ,105.03
+ ,105.86
+ ,106.04
+ ,106.03
+ ,106.13
+ ,107.21
+ ,107.66
+ ,108.08
+ ,108.76
+ ,108.26
+ ,108.71
+ ,108.65
+ ,108.61
+ ,108.86
+ ,109.54
+ ,108.22
+ ,108.77
+ ,109.9
+ ,110.13
+ ,109.6
+ ,110.42
+ ,110.6
+ ,109.73
+ ,110.72
+ ,111.08
+ ,111.14
+ ,111.01
+ ,110.56
+ ,111.57
+ ,101.73
+ ,101.63
+ ,101.43
+ ,101.34
+ ,101.01
+ ,100.89
+ ,100.93
+ ,100.77
+ ,100.3
+ ,99.86
+ ,99.71
+ ,99.93
+ ,99.88
+ ,99.92
+ ,99.87
+ ,99.63
+ ,100.05
+ ,99.88
+ ,100.11
+ ,100.05
+ ,100.07
+ ,100.2
+ ,100.21
+ ,99.76
+ ,99.41
+ ,99.24
+ ,99.65
+ ,99.7
+ ,99.79
+ ,99.84
+ ,101
+ ,101.62
+ ,101.98
+ ,101.46
+ ,102.28
+ ,102.14
+ ,102.02
+ ,102.21
+ ,101.61
+ ,102.38
+ ,102.19
+ ,102.04
+ ,101.76
+ ,101.9
+ ,102.01
+ ,102.37
+ ,103.04
+ ,103.42
+ ,103.76
+ ,104.41
+ ,104.75
+ ,104.28
+ ,103.89
+ ,104.09
+ ,103.8
+ ,105.03
+ ,105.86
+ ,106.04
+ ,106.03
+ ,106.13
+ ,107.21
+ ,107.66
+ ,108.08
+ ,108.76
+ ,108.26
+ ,108.71
+ ,108.65
+ ,108.61
+ ,108.86
+ ,109.54
+ ,108.22
+ ,108.77
+ ,109.9
+ ,110.13
+ ,109.6
+ ,110.42
+ ,110.6
+ ,109.73
+ ,110.72
+ ,111.08
+ ,111.14
+ ,111.01
+ ,110.56
+ ,111.57)
> 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] 168
> (np <- floor(n / par1))
[1] 14
> 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] 14 14 14 14 14 14 14 14 14 14 14 14
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 101.73  99.88  99.41 102.02 103.76 107.21 109.90 101.73  99.88  99.41
 [2,] 101.63  99.92  99.24 102.21 104.41 107.66 110.13 101.63  99.92  99.24
 [3,] 101.43  99.87  99.65 101.61 104.75 108.08 109.60 101.43  99.87  99.65
 [4,] 101.34  99.63  99.70 102.38 104.28 108.76 110.42 101.34  99.63  99.70
 [5,] 101.01 100.05  99.79 102.19 103.89 108.26 110.60 101.01 100.05  99.79
 [6,] 100.89  99.88  99.84 102.04 104.09 108.71 109.73 100.89  99.88  99.84
 [7,] 100.93 100.11 101.00 101.76 103.80 108.65 110.72 100.93 100.11 101.00
 [8,] 100.77 100.05 101.62 101.90 105.03 108.61 111.08 100.77 100.05 101.62
 [9,] 100.30 100.07 101.98 102.01 105.86 108.86 111.14 100.30 100.07 101.98
[10,]  99.86 100.20 101.46 102.37 106.04 109.54 111.01  99.86 100.20 101.46
[11,]  99.71 100.21 102.28 103.04 106.03 108.22 110.56  99.71 100.21 102.28
[12,]  99.93  99.76 102.14 103.42 106.13 108.77 111.57  99.93  99.76 102.14
       [,11]  [,12]  [,13]  [,14] [,15]
 [1,] 102.02 103.76 107.21 109.90    NA
 [2,] 102.21 104.41 107.66 110.13    NA
 [3,] 101.61 104.75 108.08 109.60    NA
 [4,] 102.38 104.28 108.76 110.42    NA
 [5,] 102.19 103.89 108.26 110.60    NA
 [6,] 102.04 104.09 108.71 109.73    NA
 [7,] 101.76 103.80 108.65 110.72    NA
 [8,] 101.90 105.03 108.61 111.08    NA
 [9,] 102.01 105.86 108.86 111.14    NA
[10,] 102.37 106.04 109.54 111.01    NA
[11,] 103.04 106.03 108.22 110.56    NA
[12,] 103.42 106.13 108.77 111.57    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] -0.10  0.04 -0.17  0.19  0.65  0.45  0.23 -0.10  0.04 -0.17  0.19  0.65
 [2,] -0.20 -0.05  0.41 -0.60  0.34  0.42 -0.53 -0.20 -0.05  0.41 -0.60  0.34
 [3,] -0.09 -0.24  0.05  0.77 -0.47  0.68  0.82 -0.09 -0.24  0.05  0.77 -0.47
 [4,] -0.33  0.42  0.09 -0.19 -0.39 -0.50  0.18 -0.33  0.42  0.09 -0.19 -0.39
 [5,] -0.12 -0.17  0.05 -0.15  0.20  0.45 -0.87 -0.12 -0.17  0.05 -0.15  0.20
 [6,]  0.04  0.23  1.16 -0.28 -0.29 -0.06  0.99  0.04  0.23  1.16 -0.28 -0.29
 [7,] -0.16 -0.06  0.62  0.14  1.23 -0.04  0.36 -0.16 -0.06  0.62  0.14  1.23
 [8,] -0.47  0.02  0.36  0.11  0.83  0.25  0.06 -0.47  0.02  0.36  0.11  0.83
 [9,] -0.44  0.13 -0.52  0.36  0.18  0.68 -0.13 -0.44  0.13 -0.52  0.36  0.18
[10,] -0.15  0.01  0.82  0.67 -0.01 -1.32 -0.45 -0.15  0.01  0.82  0.67 -0.01
[11,]  0.22 -0.45 -0.14  0.38  0.10  0.55  1.01  0.22 -0.45 -0.14  0.38  0.10
[12,] -0.05 -0.35 -0.12  0.34  1.08  1.13 -9.84 -0.05 -0.35 -0.12  0.34  1.08
      [,13] [,14] [,15]
 [1,]  0.45  0.23    NA
 [2,]  0.42 -0.53    NA
 [3,]  0.68  0.82    NA
 [4,] -0.50  0.18    NA
 [5,]  0.45 -0.87    NA
 [6,] -0.06  0.99    NA
 [7,] -0.04  0.36    NA
 [8,]  0.25  0.06    NA
 [9,]  0.68 -0.13    NA
[10,] -1.32 -0.45    NA
[11,]  0.55  1.01    NA
[12,]  1.13    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/1eyru1229595867.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/2d9mm1229595867.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/3mkfz1229595867.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/4nqi71229595867.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,]  99.41  99.24  99.65  99.63  99.79  99.84 100.11 100.05 100.07  99.86
[2,]  99.88  99.92  99.87  99.70 100.05  99.88 100.93 100.77 100.30 100.20
[3,] 102.02 102.21 101.61 102.38 102.19 102.04 101.76 101.90 102.01 102.37
[4,] 107.21 107.66 108.08 108.76 108.26 108.71 108.65 108.61 108.86 109.54
[5,] 109.90 110.13 109.60 110.42 110.60 109.73 110.72 111.08 111.14 111.01
      [,11]  [,12]
[1,]  99.71  99.76
[2,] 100.21  99.93
[3,] 103.04 103.42
[4,] 108.22 108.77
[5,] 110.56 111.57

$n
 [1] 14 14 14 14 14 14 14 14 14 14 14 14

$conf
          [,1]     [,2]      [,3]     [,4]      [,5]      [,6]      [,7]
[1,]  98.92474  98.9416  98.14314  98.5542  98.72314  98.31133  98.50005
[2,] 105.11526 105.4784 105.07686 106.2058 105.65686 105.76867 105.01995
          [,8]      [,9]     [,10]    [,11]    [,12]
[1,]  98.58938  98.39535  98.42597  99.6576  99.6871
[2,] 105.21062 105.62465 106.31403 106.4224 107.1529

$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(99.41, 99.88, 102.02, 107.21, 109.9, 99.24,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/59zz11229595867.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,] -0.17 -0.60 -0.47 -0.50 -0.17 -0.29 -0.16 -0.47 -0.52 -1.32 -0.45 -0.35
[2,] -0.10 -0.53 -0.24 -0.39 -0.17 -0.28 -0.06  0.02 -0.44 -0.45 -0.14 -0.12
[3,]  0.19 -0.05  0.05 -0.19 -0.12  0.04  0.14  0.11  0.13 -0.01  0.22 -0.05
[4,]  0.45  0.41  0.77  0.18  0.20  0.99  0.62  0.36  0.36  0.67  0.55  1.08
[5,]  0.65  0.42  0.82  0.42  0.45  1.16  1.23  0.83  0.68  0.82  1.01  1.13

$n
 [1] 14 14 14 14 14 14 14 14 14 14 14 13

$conf
            [,1]       [,2]       [,3]        [,4]        [,5]       [,6]
[1,] -0.04225002 -0.4469364 -0.3764955 -0.43069547 -0.27624092 -0.4962864
[2,]  0.42225002  0.3469364  0.4764955  0.05069547  0.03624092  0.5762864
           [,7]        [,8]       [,9]      [,10]      [,11]      [,12]
[1,] -0.1471455 -0.03357274 -0.2078182 -0.4829455 -0.0713682 -0.5758558
[2,]  0.4271455  0.25357274  0.4678182  0.4629455  0.5113682  0.4758558

$out
[1] -0.87 -0.87 -9.84

$group
[1]  5  5 12

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

Warning message:
In bxp(list(stats = c(-0.170000000000002, -0.100000000000009, 0.189999999999998,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/697q61229595867.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]    [,7]    [,8]    [,9]
[1,]  99.710  99.630  99.240 101.610 103.760 107.660 109.600  99.710  99.630
[2,] 100.115  99.875  99.675 101.955 103.990 108.150 110.015 100.115  99.875
[3,] 100.910  99.985 100.420 102.115 104.580 108.630 110.580 100.910  99.985
[4,] 101.385 100.090 101.800 102.375 105.945 108.765 111.045 101.385 100.090
[5,] 101.730 100.210 102.280 102.380 106.130 109.540 111.570 101.730 100.210
       [,10]   [,11]   [,12]   [,13]   [,14] [,15]
[1,]  99.240 101.610 103.760 107.660 109.600    NA
[2,]  99.675 101.955 103.990 108.150 110.015    NA
[3,] 100.420 102.115 104.580 108.630 110.580    NA
[4,] 101.800 102.375 105.945 108.765 111.045    NA
[5,] 102.280 102.380 106.130 109.540 111.570    NA

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12  0

$conf
         [,1]      [,2]      [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 100.3307  99.88694  99.45077 101.9234 103.6883 108.3495 110.1102 100.3307
[2,] 101.4893 100.08306 101.38923 102.3066 105.4717 108.9105 111.0498 101.4893
          [,9]     [,10]    [,11]    [,12]    [,13]    [,14] [,15]
[1,]  99.88694  99.45077 101.9234 103.6883 108.3495 110.1102    NA
[2,] 100.08306 101.38923 102.3066 105.4717 108.9105 111.0498    NA

$out
[1] 103.04 103.42 107.21 103.04 103.42 107.21

$group
[1]  4  4  6 11 11 13

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

Warning message:
In bxp(list(stats = c(99.71, 100.115, 100.91, 101.385, 101.73, 99.63,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/78erd1229595867.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,] 103.4157 101.610 103.3450
[2,] 103.5986 101.955 103.7412
[3,] 103.8200 102.115 104.0338
[4,] 104.3050 102.375 104.3225
[5,] 104.5314 102.380 104.5900

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 103.4978 101.9234 103.7686
[2,] 104.1422 102.3066 104.2989

$out
[1] 103.04 103.42

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(103.415714285714, 103.598571428571, 103.82,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1eyru1229595867.ps tmp/1eyru1229595867.png")
> system("convert tmp/2d9mm1229595867.ps tmp/2d9mm1229595867.png")
> system("convert tmp/3mkfz1229595867.ps tmp/3mkfz1229595867.png")
> system("convert tmp/4nqi71229595867.ps tmp/4nqi71229595867.png")
> system("convert tmp/59zz11229595867.ps tmp/59zz11229595867.png")
> system("convert tmp/697q61229595867.ps tmp/697q61229595867.png")
> system("convert tmp/78erd1229595867.ps tmp/78erd1229595867.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.276   0.985   2.126