ea_phi_uniGaussian_DL_bounds: Obtain bounds for phi function

In rchan26/hierarchicalFusion: Divide-and-Conquer Monte Carlo Fusion

Obtain bounds for phi function

Description

Finds the lower and upper bounds of the phi function between an interval

Usage

ea_phi_uniGaussian_DL_bounds(mean, sd, beta, precondition, lower, upper)

Arguments

mean	mean value
sd	standard deviation value
beta	real value
precondition	precondition value
lower	lower end of interval
upper	upper end of interval

Value

A list of components

LB

lower bound of phi

UB

upper bound of phi

Examples

mu <- 0.423
sd <- 3.231
beta <- 0.8693
precondition <- 1.243
lower <- -2.823
upper <- 4.322
curve(ea_phi_uniGaussian_DL(x, mu, sd, beta, precondition), lower, upper)
abline(h=ea_phi_uniGaussian_DL_bounds(mean = mu,
                                      sd = sd,
                                      beta = beta,
                                      lower = lower,
                                      upper = upper,
                                      precondition = precondition),
       col = 'red', lty = 2)

# another example where the mean is not in the interval
mu <- 0.423
sd <- 3.231
beta <- 0.8693
precondition <- 1.243
lower <- 2.823
upper <- 5.322
curve(ea_phi_uniGaussian_DL(x, mu, sd, beta, precondition), lower, upper)
abline(h=ea_phi_uniGaussian_DL_bounds(mean = mu,
                                      sd = sd,
                                      beta = beta,
                                      precondition = precondition,
                                      lower = lower,
                                      upper = upper),
       col = 'red', lty = 2)

rchan26/hierarchicalFusion documentation built on Sept. 11, 2022, 10:30 p.m.