In the first semester of my MSc. studies, we developed a phyton version of the Gibbs Sampler and Metropolis-Hastings Algorithm from the scratch. We described our results and analysis in a report.
For implementing Metropolis-Hastings algorithm simulated data were used with k=3 covariates and a constant term. A multivariate normal distribution was taken as a prior for , N(0,) where is a diagonal matrix with on the
diagonal. After considering a number of values for , = 1 was chosen. The posterior distribution was simulated using as proposal a multivariate normal distribution centered at the current update of and with a covariance matrix given by the inverse of Fisher information evaluated at the current update.
Gibbs Sampler for binary data
For this analyse, we simulated N independent binary random variables where each yi comes from a Bernoulli distribution with probability of success .
About
In the first semester of my MSc. studies, we developed a phyton version of the Gibbs Sampler and Metropolis-Hastings Algorithm from the scratch. We described our results and analysis in a report.