ace.bnds.lipid <- ACE_bounds_posterior(158, 14, 0, 0, 52, 12, 23, 78,
prior = c( rep(1, 2), rep(0, 2), rep(1, 4)))
ACE_bounds_triangle.plot(ace.bnds.lipid, "Bounds on ACE for Lipid Data")
## Not run:
# ace.bnds.lipid <- ACE_bounds_posterior(158, 14, 0, 0, 52, 12, 23, 78,
# prior = c( rep(1, 2), rep(0, 2), rep(1, 4)), num.sims = 2e4)
# ACE_bounds_triangle.plot(ace.bnds.lipid, "Bounds on ACE for Lipid Data") ## End(Not run)
Run the code above in your browser using DataLab