nt-williams

lmtp_type <- \(x) structure(x, class = class(x[[1]]))

rubins_rules <- \(x, ...) UseMethod("rubins_rules", lmtp_type(x))

rubins_rules.lmtp <- function(lmtps, label, alpha = 0.05) {
  thetas <- purrr::map_dbl(lmtps, "theta")
  vw <- mean(purrr::map_dbl(lmtps, \(x) x$standard_error^2))
  vb <- var(thetas)
  theta <- mean(thetas)
  se <- pooled_se(vw, vb, length(thetas))