estimate_lambda {SpatialEpi}

Estimate lambda values

Description

Internal function to estimate values of lambda needed for MCMC_simulation and prior probability of k clusters/anti-clusters for k=0,...,J

Usage

estimate_lambda(n.sim, J, prior.z, overlap, pi0)

Arguments

n.sim

number of importance sampling iterations

J

maximum number of clusters/anti-clusters to consider

prior.z

prior probability of each single zone

overlap

output of create_geo_objects(): list with two elements: presence which lists for each area all the single zones it is present in and cluster_list for each single zone its component areas

pi0

prior probability of no clusters

Value

estimates of lambda and prior.j

References

Wakefield J. and Kim A.Y. (2013) A Bayesian model for cluster detection. Biostatistics, 14, 752–765.