MCMC simulation to estimate prior and posterior quantities by sampling configurations.
1 2  MCMC_simulation(n_sim, pattern, theta_init, overlap, cluster_coords,
p_moves_orig, J, lkhd_z, lambda)

n_sim 
number of MCMC iterations 
pattern 
alternating pattern between unif and prop prior on single zones 
theta_init 
initial configuration 
overlap 
output of 
cluster_coords 
output of 
p_moves_orig 
probability of sampling each of the 5 possible moves to explore sample space 
J 
maximum number of clusters/anticlusters to consider 
lkhd_z 
values associated with each single zone to use in MetropolisHastings ratio 
lambda 
lambda from definition of prior on single zones 
sample 
sampled configurations 
move_trace 
trace of moves (1 = growth, 2 = trim, 3 = recenter, 4 = death, 5 = birth) 
accpt_trace 
trace of acceptance (0 = not accepted) 
ratio_trace 
trace of MetropolisHastings ratio 
Albert Y. Kim
Wakefield J. and Kim A.Y. (2013) A Bayesian model for cluster detection. Biostatistics, 14, 752–765.
create_geo_objects
, process_MCMC_sample
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.