R/likelihood_ratio_gam.R

In BriceonWiley/comparingCIs.R:

Defines functions likelihood_ratio_gam

Documented in likelihood_ratio_gam

likelihood_ratio_gam <- function(gam, num, data1, data2, crit) {
  denom <- likelihood_gam(gam, data1, data2)
  2 * (log(num) - log(denom)) - crit
}