R/likelihood_phi.R

In BriceonWiley/IntegratedLikelihood.R:

likelihood_phi <- function(phi, data, init) {
  ub <- init + 5
  pcubature(
    likelihood_orig_phi,
    lowerLimit = c(1e-5, 1e-5, 1e-5),
    upperLimit = c(ub, 1 - 1e-5, 1 - 1e-5),
    phi = phi, data = data, vectorInterface = TRUE
  )$integral
}