Full Bayesian Parameter Estimation - sampling
Bayesian parameter estimation The previous post introduced the basic concept of Bayesian parameter estimation. If you haven't read it, please check that post before this one. When Bayesian estimation is not simple There are two cases when Bayesian estimation is tractable. The simplest is when a prior distribution is a conjugate prior. Computing the normalising constant of a single or low-dimensional parameter of \( \theta \) is also tractable with a grid approximation. In this post, I will stick to the Pikachu encounter rate example like in the previous post too. Imagine we have analysed 10 routes, each with its own Pikachu encounter rate: \[ \boldsymbol{\theta} = [\theta_1, \dots, \theta_{10}], \quad \theta_k \in [0, 1] \] On each route \( k \), we observed \( n_k \) Pikachu out of \( N_k \) Pokemon, so the likelihood is a product of Binomials, \( P(D|\boldsymbol{\theta}) \propto \prod_{k=1}^{10} \theta_...