data <- data.table::data.table(
x = rnorm(1e3), y = as.factor(rbinom(1e3, size = 1, prob=0.5))
)
m <- estimate_covar_model_full_cond(data)
samples <- sample_covar_parametric_model(n=10, model = m)
print(head(samples))
Run the code above in your browser using DataLab