Description
Activity 06 Today you will work with the same data as last week, and estinmate the average
radius of the sample of Kepler’s detected extrasolar planets.
You will use the same file “Planet_Radii.csv” that can be found in ACTIVITY05 directory on
Canvas and that includes the radii of a large sample of extrasolar planets.
You will write a code to use a Gibbs sampling algorithm to sample from the posterior distribution
of both parameters µ and λ of the LogNormal distribution. You will do this for both the improper
Jeffreys prior and the more informative priors we discussed in class on Monday. Perform a sensitivity
analysis for the proper priors.
Use the posterior to estimate the mean radius of the sample of extrasolar planets, together with
the 95% credibility interval.
Explore how the result changes by varying the chain length: 100, 1000, 10000 when you simulate
with the Gibbs sampler.
[ ]:
1