Background

The Hubble constant is a fundamental cosmological parameter used to describe the size and age of the universe.  It relates the distance and velocity of objects relative to earth via the relationship v = ηd. Astrophysicists estimate the Hubble constant η from observed ratios H = for a sample of galaxies. However, observations of H for any given galaxy can vary considerably due to measurement error.

Though seemingly simple, measurements of velocity and distance are challenging to obtain with accuracy, and continue to be refined (see a recent article). Velocity is calculated from spectroscopic measurements, which are less complicated to obtain compared with distance measurements.  For the latter, “[r]eliable ‘distance indicators’, such as variable stars and supernovae, must be found in ... a sample of galaxies that are far enough away that motions due to local gravitational influences are negligibly small” (quote from Space Science Short).

The primary objective of the Hubble Space Telescope (HST), launched in 1990, is to collect data that will lead to more precise estimates of the Hubble constant. In 2001, Freedman et al. published a set of measurements based on novel distance indicators and HST data; distance-velocity ratios from measurements for 36 type 1a supernovae using their method of calculation are shown in Figure 1.

Figure 1:  Histogram of H values from 36 type 1a supernovae, with a smooth overlay to aid in visualizing the shape of the distribution. The exact values in the sample are shown as ticks below the histogram. The sample mean and variance are = 72.1861 and SH(2)  = 23.7858.

Objective

Your goal in this assignment is to estimate the age of the universe based on an estimate of the Hubble constant η and to quantify the quality of the estimate.  The age of the universe can be

expressed in years as

c

θ =

where c is the conversion factor = 978, 440, 076, 094.

Assumptions

Assume throughout that H1 , . . . , H36 gamma(α, β), so that the common density is:

f (h; α, β) = hα − 1 e − h > 0

Note that this exam works with the gamma density as parameterized in the textbook – that is, with shape parameter α and scale parameter β .

Assume further that the common mean is in fact the Hubble constant; denote this by η = E[Hi]. Consider the estimators ηˆ = = and θˆ = .

Problems

1.  The following parts guide you through finding the distribution of θˆ.              (a)  Show that c − 1 has a gamma ←nα, _ distribution. Write its density.

(b)  Show that θˆ = c − 1  has an inverse gamma distribution with density

g(u; α, β) = u −nα − 1 e −

For short, U ~ IG ←nα, _.  (Hint: for x > 0, the function is monotone decreasing.)

2.  The following parts guide you through assessing the bias of θˆ.

(a) Write the relationship between the parameters α, β and the parameter of interest θ .

(b)  Review the Wikipedia page on the inverse gamma. Pay close attention to how (their) β

appears in the density function (we have used a ’shape’ parameter and they have used a ’rate’ parameter; one can convert directly between the two by substituting the reciprocal of one for the other in the density, CDF, moments, etc.).  Find an expression for the expected value of θˆ in terms of c, α, β, and n.

(c)  Find an expression for bias ← _θˆ in terms of n, c, α, and β .

3.  The following parts have you apply and interpret your results.

(a)  Compute an estimate for the age of the universe from the data provided.

(b)  Do you expect this is an underestimate, an overestimate, or neither? Explain why in 1-2

sentences.