首页 » Python辅导 » Python辅导 | MATH3871/MATH5970 Bayesian Inference and Computation

Python辅导 | MATH3871/MATH5970 Bayesian Inference and Computation

MATH3871/MATH5970
Bayesian Inference and Computation
Introduction
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:
1
y=c(4,5,4,1,0,4,3,4,0,6,
3,3,4,0,2,6,3,3,5,4,5,3,1,4,4,1,5,5,3,4,2,5,2,2,3,4,2,1,3,2,
1,1,1,1,1,3,0,0,1,0,1,1,0,0,3,1,0,3,2,2,0,1,1,1,0,1,0,1,0,0,
0,2,1,0,0,0,1,1,0,2,2,3,1,1,2,1,1,1,1,2,4,2,0,0,0,1,4,0,0,0,
1,0,0,0,0,0,1,0,0,1,0,0)
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;
2


程序辅导定制C/C++/JAVA/安卓/PYTHON/留学生/PHP/APP开发/MATLAB


本网站支持 Alipay WeChatPay PayPal等支付方式

E-mail: vipdue@outlook.com  微信号:vipnxx


如果您使用手机请先保存二维码,微信识别。如果用电脑,直接掏出手机果断扫描。

blank

发表评论

您的电子邮箱地址不会被公开。