Homework 9.2: MLE of microtubule catastrophe data (40 pts)

Dataset download


Refresh yourself about the microtubule catastrophe data we explored in previous homeworks. We will again work with this data set here.

a) In their paper, Gardner, Zanic, and coworkers modeled microtubule catastrophe times as Gamma distributed. Perform a maximum likelihood estimate for the parameters of the Gamma distribution. Because you showed in a previous homework that there is little difference between labeled and unlabeled tubulin, you only need to work this out for the labeled tubulin now and in part (b). Be sure to include confidence intervals or a confidence region for your MLE and discuss the method you used to get the confidence intervals or confidence region.

b) Obtain a maximum likelihood estimate for the parameters \(\beta_1\) and \(\beta_2\) from the model you derived in homework 5.2. As a reminder, you derived that the PDF for microtubule catastrophe times is

\begin{align} f(t;\beta_1, \beta_2) = \frac{\beta_1\beta_2}{\beta_2 - \beta_1}\left(\mathrm{e}^{-\beta_1 t} - \mathrm{e}^{-\beta_2 t}\right). \end{align}

Be careful; this is a very tricky calculation. It is possible to make significant analytical progress, and this can help you find the MLE. Be sure to think about what happens when \(\beta_1 \approx \beta_2\). You may also need to think about how you will handle the log of sums of exponentials.

c) Compute bootstrap confidence intervals. This calculation has difficult subtleties that make it even trickier than computing the MLE. This part of the problem is optional and your response will not be graded.