Run a Bayesian functional data model under an instrinsic GMRF prior whose precision parameters employ a DP prior for a COUNT data response type where: y ~ poisson(E*exp(Psi)) Psi ~ N(gamma,tau_e^-1) which is a Poisson-lognormal model
Useful links