dat <- data.frame(
cov1 = rnorm(50, 10, 1),
cov2 = rnorm(50, 7, 1),
cov3 = rnorm(50, 5, 1),
treated = rep(c(0, 1), c(35, 15))
)
out <- VCG_sampler(treated ~ cov1 + cov2 + cov3, data=dat, n=5, plot=FALSE)
plot_var(out, what='cov1', group='VCG')
plot_var(out, what='cov2', group='VCG')
Run the code above in your browser using DataLab