#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
#
#-----

source('utils.R')

# Ex 16.27:
#
mu_0 = 16
n = 5 # selected every hour
delta = 0.1
h = 0.2
data = c( 15.992, 16.051, 16.066, 15.912, 16.030,
          16.060, 15.982, 15.899, 16.038, 16.074,
          16.029, 15.935, 16.032, 15.960, 16.055 )

CUMSUM( data, mu_0, delta/2 )


# Ex 16.28:
#
mu_0 = 0.75
n = 4
delta = 0.002
h = 0.003

data = c( .7507, .7504, .7492, .7501, .7503, .7510, .7490, .7497, .7488, .7504, .7516,
          .7472, .7489, .7483, .7471, .7498, .7460, .7482, .7470, .7493, .7462, .7481 )

CUMSUM( data, mu_0, delta/2 )