ea_phi_mixG_DL_LB: Obtain the global lower bound for phi function

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

Obtain the global lower bound for phi function

Description

Finds the global bound of the phi function between a given interval

Usage

ea_phi_mixG_DL_LB(
  n_comp,
  weights,
  means,
  sds,
  beta,
  precondition,
  bounds_multiplier
)

Arguments

n_comp	integer number of components of mixture Gaussian
weights	vector: weights of mixture Gaussian
means	vector: means of mixture Gassuan
sds	vector: st.devs of mixture Gaussian
beta	real value
normalised	boolean value to determine if normalisation constant is calculated

Value

The global lower bound of phi

Examples

weights <- c(0.4, 0.6)
means <- c(-3, 6)
sds <- c(1, 2)
beta <- 1
precondition <- 1
curve(ea_phi_mixG_DL(x,
                     n_comp = 2,
                     weights = weights,
                     means = means,
                     sds = sds,
                     beta = beta,
                     precondition = precondition),
      from = -20, to = 25, n = 10000,
      ylab = 'phi')
PHI <- ea_phi_mixG_DL_LB(n_comp = 2,
                         weights = weights,
                         means = means,
                         sds = sds,
                         beta = beta,
                         precondition = precondition,
                         bounds_multiplier = 1)
abline(h = PHI, col = 'red')

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