#-- Generate Example Data
set.seed(2023)
dat <- simdat(n = c(300, 300, 300), effect = 1, sigma_Y = 1)
head(dat)
formula <- Y - f ~ X1
#-- Fit IPD
fit <- ipd(formula, method = "postpi_analytic", model = "ols",
data = dat, label = "set_label")
#-- Tidy Output
tidy(fit)
Run the code above in your browser using DataLab