Compute the bounds for the likelihood ratio under the marginal sensitivity model

Compute the bounds for the likelihood ratio under the marginal sensitivity model

lr_bnds(
  estimand = c("treated", "control"),
  type = c("ate", "att", "atc"),
  Gamma,
  ps,
  p1,
  p0
)