Estimate posteriors for version 7 model
1 2 3 4 5 | cause_mh_v7(dat, mix_grid, rho, seed, gamma_prior_func = function(b) {
dnorm(b, 0, 0.6, log = TRUE) }, gamma_prime_prior_func = function(b) {
dnorm(b, 0, 0.6, log = TRUE) }, q_prior_func = function(q) {
dbeta(q, 0.1, 1, log = TRUE) }, nchain = 5, niter_perchain = 10000,
burn_in = niter_perchain/2, thin = 5)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.