## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>",
                      fig.width = 6, fig.height = 4)
set.seed(20260821)
library(metaGLMM)

## ----data-fit-----------------------------------------------------------------
dat <- data.frame(
  study = paste0("Study ", 1:8),
  estimate = c(-0.80, -0.30, 0.10, 0.70, 1.10, -0.50, 0.40, 1.30),
  moderator = c(0, 1, 0, 1, 0, 1, 0, 1),
  vi = c(0.04, 0.05, 0.03, 0.06, 0.04, 0.05, 0.03, 0.05),
  ni = c(100, 90, 120, 80, 110, 95, 130, 85)
)

fit <- metaGLMM(
  estimate ~ moderator,
  data = dat,
  vi = dat$vi,
  ni = dat$ni,
  tau2 = NA,
  tau2_var = TRUE,
  family = gaussian(link = "identity"),
  fast = TRUE
)
summary(fit)

## ----wald---------------------------------------------------------------------
confint(fit, method = "wald")
confint(fit, parm = "moderator", level = 0.90, method = "wald")

## ----profile, eval=FALSE------------------------------------------------------
# profile_ci <- confint(fit, parm = "moderator",
#                       method = "profile", level = 0.95)
# profile_ci
# 
# sbc_ci <- confint(fit, parm = "moderator",
#                   method = "SBC", level = 0.95)
# sbc_ci

## ----methods------------------------------------------------------------------
coef(fit)
vcov(fit)
logLik(fit)
nobs(fit)
fit$tau2
fit$tau
stopifnot(is.finite(fit$tau), fit$tau > 0)