#-- 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")
#-- Augment Data
augmented_df <- augment(fit)
head(augmented_df)
Run the code above in your browser using DataLab