Homework is at the time and date indicated. Please make sure your name is on every page of the assignment, and it is clearly marked which question you are answering. Your response will be graded for correctness and clarity.
1. This problem continues with the HCB pollution data from the Assignments 6 and 7.
|
Surface |
Bottom |
|
3.74 |
5.44 |
|
4.61 |
6.88 |
|
4.00 |
5.37 |
|
4.67 |
5.44 |
|
4.87 |
5.03 |
|
5.12 |
6.48 |
|
4.52 |
3.89 |
|
5.29 |
5.85 |
|
5.74 |
6.85 |
|
5.48 |
7.16 |
As with Assignments 6 and 7, assume the observations are normally distributed with unknown depth- specific means Θ” and Θ# and precisions Rs and Rb. Assume that experts have provided the following prior information based on previous studies.
● The unknown means Θ” and Θ# are independent and normally distributed with mean μ and standard deviation τ. The unknown precisions 𝛲 and 𝛲 are independent of Θs and Θ and
“## have gamma distributions with shape 𝛼 and scale 𝛽.
● Experts specified a 95% prior credible interval of [3, 9] for Θ” and Θ#. A good fit to this credible interval is obtained by setting the prior mean to μ =6 and the prior standard deviation to τ=1.5.
● A 95% prior credible interval of [0.75, 2.0] is given for the unknown standard deviations 𝛴” and 𝛴 . This translates to a credible interval of [0.25, 1.8] for 𝛲 = Σ*+ and 𝛲 = Σ*+. A
# “”##
good fit to this credible interval is obtained by setting the prior shape to 𝛼 = 4.5 and the prior
scale to 𝛽 = 0.19.
Find the following conditional distributions.
- ● The conditional distribution for Θ” given the other parameters and the observations.
- ● The conditional distribution for Θ# given the other parameters and the observations.
- ● The conditional distribution for 𝛲 given the other parameters and the observations. “
- ● The conditional distribution for 𝛲 given the other parameters and the observations. #
2. Using the distributions you found in Part 1, draw 10,000 Gibbs samples of (Θ , Θ , 𝛲 , 𝛲 ). Estimate
“#”#
3. Do a traceplot of Θ# – Θ”. Find the autocorrelation function of Θ# – Θ” and the effective sample size
90%credibleintervalsforΘ ,Θ ,Σ =𝛲*,/+,Σ =𝛲*,/+,andΘ -Θ . “#””###”
for your Monte Carlo sample for Θ# – Θ”.
4. Comment on your results. Compare with Assignment 6.