#-- 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")
#-- Summarize Output
summ_fit <- summary(fit)
summ_fit
Run the code above in your browser using DataLab