Python辅导 | MATH3871/MATH5970 Bayesian Inference and Computation

Bayesian Inference and Computation
This is the second assignment for course MATH3871/5960.
The assignment consists in two parts: a quiz and a programming part. The deadline for both
is the 6th of November 2019 at 6PM (Sydney time). The deadline is strict.
For this assignment, it is possible to work in groups of maximum 4 people. However, it is not
compulsory: groups may be formed by 1, 2, 3 or 4 people.
Each student should complete the quiz individually, while the programming part must be
submitted once for each group. The members of the group should be clearly stated at the
beginning of the code. Every results should be properly commented. Only .py or .R (or .rmd)
files will be accepted. The programming part should run without errors to be evaluated.
The quiz worths 6 points, while the programming part worths 14 points.
The weight for this assignment on the final mark is 20%.
These instructions are strict.
Assignement 2
The annual number of serious accidents in the mines of Great Britain were recorded every year
from 1851 to 1962. We want to build a statistical model to evaluate if during those 112 years
there has been a change in the rate of occurrence of events, which could perhaps be due to
legislative changes to protect security.
Let (Y1, . . . , Ym−1, Ym, Ym+1, . . . , Yn) be the number of annual incidents, with
Yi ∼ P oi(λ) i = 1, 2, . . . , m
Yj ∼ P oi(φ) j = m + 1, . . . , n
The unknown parameters of the model are (φ, λ, m), while n is known and equal to 112.
Since the observations are supposed to follow a Poisson distribution, it is reasonable to consider
Gamma priors for (φ, λ):
λ ∼ Gamma(α, β)
φ ∼ Gamma(a, b)
On the other side, the parameter m can have a uniform distribution:
m ∼ Unif(0, n − 1)
The data are the following:
Write an MCMC algorithm to perform a Bayesian analysis for this model. It could be useful
to procede with the following steps
• write the likelihood function;
• write the full posterior distribution;
• code a Gibbs sampler with N = 10000 simulations and with a burnin of 1000 values;
• plot the simulates marginal chain for each parameter and comment;
• analyse the convergence of the chains with at least two tools;