\documentclass[12pt]{article} 
%\usepackage{amsbsy} % for \boldsymbol and \pmb 
\usepackage{graphicx} % To include pdf files!
\usepackage{amsmath}
\usepackage{amsbsy}
\usepackage{amsfonts}
\usepackage[colorlinks=true, pdfstartview=FitV, linkcolor=blue, citecolor=blue, urlcolor=blue]{hyperref} % For links
\usepackage{fullpage}
%\pagestyle{empty} % No page numbers


\begin{document}
%\enlargethispage*{1000 pt} 

\begin{center}   
{\Large \textbf{STA 312s19 Assignment Five}}\footnote{This assignment was prepared by  \href{http://www.utstat.toronto.edu/~brunner}{Jerry Brunner},
Department of Mathematical and Computational Sciences, University of Toronto. It is licensed under a 
\href{http://creativecommons.org/licenses/by-sa/3.0/deed.en_US}
     {Creative Commons Attribution - ShareAlike 3.0 Unported License}. Use any part of it as you like and share the result freely. The \LaTeX~source code is available from the course website:
\href{http://www.utstat.toronto.edu/~brunner/oldclass/312s19} {\texttt{http://www.utstat.toronto.edu/$^\sim$brunner/oldclass/312s19}}}
\vspace{1 mm}
\end{center}

\noindent
The paper and pencil part of this assignment is not to be handed in. It is practice for Quiz~5 on February 11th. The R part may be handed in as part of the quiz. \textbf{Bring hard copy of your printout to the quiz}. Do not write anything on your printout in advance except possibly your name and student number.

\begin{enumerate} 

\item \label{paperpencil} Consider a random sample of right censored data from an exponential distribution with parameter $\lambda>0$. In addition to the data $T_1, \ldots, T_n$, we have indicators $\delta_1, \ldots, \delta_n$, where $\delta_j=1$ if  $T_j$ is a failure time, and $\delta_j=0$ if $T_j$ is a censoring time.
    \begin{enumerate}
        \item \label{mle} Derive the maximum likelihood estimate of $\lambda$. Show your work. Include the second derivative test. Circle your final answer. It is on the lecture slides. What happens if \emph{all} the observations are censored?
        \item \label{vhat} Give a formula for the estimated asymptotic variance of $\widehat{\lambda}$. Show your work. Circle your answer.
        \item \label{cilambda} Give a 95\% confidence interval for $\lambda$. Your answer is two formulas, one for the lower confidence limit and one for the upper confidence limit.  
        \item \label{Shat} Give a formula for the MLE of the survival function $S(t)$. Your answer is a formula that could be computed from the sample data for any given $t$.  
        \item \label{ciS1} Using the one-variable delta method (see formula sheet), derive a 95\% confidence interval for the survival function $S(t)$. As above, your answer is two formulas, one for the lower confidence limit and one for the upper confidence limit.
        \item \label{ciS2} Now do the same thing a different way. Writing your confidence interval for $\lambda$ as $0.95 \approx P(A < \lambda < B)$, multiply by $-t$ and take the exponential function to obtain a different confidence interval for $S(t)$, one that is also valid. Is it possible for these confidence limits to be outside the range zero to one?
    \end{enumerate}

\item \label{calculator} The file 
\href{http://www.utstat.toronto.edu/~brunner/data/legal/expo.data2.txt}
     {\texttt{http://www.utstat.toronto.edu/$\sim$brunner/data/legal/expo.data2.txt}}
contains a random sample from an exponential distribution with right censoring. Using R as a high-powered calculator (no numerical search yet),
    \begin{enumerate}
        \item \label{handmle} Give the maximum likelihood estimate of $\lambda$. The answer is a number. You are calculating your answer to Question~\ref{mle}.
        \item \label{handvhat} What is the estimated asymptotic variance of $\widehat{\lambda}$? The answer is a number. You are calculating your answer to Question~\ref{vhat}.
        \item Give a 95\% confidence interval for $\lambda$. Your answer is two numbers, the lower confidence limit and the upper confidence limit. You are calculating your answer to Question~\ref{cilambda}.  
        \item Based on your answers to Questions~\ref{Shat} and~\ref{ciS1}, calculate your estimated $S(t)$ and the confidence interval for the following values of $t$: \texttt{seq(from=0,to=3,by=0.1)}. Bind $t$, $\widehat{S}(t)$ the lower confidence limit and the upper confidence limit into a matrix with 4 columns.
        \item Now do the same thing, using your confidence limits from Question~\ref{ciS2}. The result is another $4 \times 31$ matrix. Which one do you like more and why? % Second one stays inside range 0 to 1.
    \end{enumerate}

\item \label{optim} Using the data from Question~\ref{calculator}, 
    \begin{enumerate}
        \item Obtain the maximum likelihood estimate of $\lambda$ by a numerical search. Compare your answer to Question~\ref{handmle}.
        \item Obtain the estimated variance of $\widehat{\lambda}$ from the output of your search. If you think about it, you will realize that $-\ell^{\prime\prime}(\hat{\lambda})$ is a number on your printout. Display the estimated variance on your printout.  Compare your answer to Question~\ref{handvhat}.
    \end{enumerate}
After that, what you'd do would be the same as in Question~\ref{calculator}. You don't have to do it twice. 
 









% \pagebreak %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\end{enumerate} % End of all the questions

\noindent
\textbf{Bring your printout for Questions~\ref{calculator} and~\ref{optim} to the quiz.} All the requested numbers and the code that produced them should appear on your printout. Do not write anything on your printout in advance except possibly your name and student number. 


\end{document}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% See work and quiz files for R work.