2 Right Censoring and Kaplan-Meier Estimator - NCSU Statistics
2 Right Censoring and Kaplan-Meier Estimator - NCSU Statistics
2 Right Censoring and Kaplan-Meier Estimator - NCSU Statistics
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 2 ST 745, Daowen Zhang<br />
which can be estimated by<br />
= E[var(�m(i−1)|H)] + var[E(�m(i − 1)|H)]<br />
� �<br />
m(i − 1)[1 − m(i − 1)]<br />
= E<br />
+var[m(i−1)] n(i − 1)<br />
� �<br />
1<br />
= m(i − 1)[1 − m(i − 1)]E ,<br />
n(i − 1)<br />
�m(i − 1)[1 − �m(i − 1)]<br />
.<br />
n(i − 1)<br />
Hence the variance of � φi = log(�qi) can be approximately estimated by<br />
� �2 1 �m(i − 1)[1 − �m(i − 1)]<br />
=<br />
n(i − 1)<br />
�qi<br />
�m(i − 1)<br />
[1 − �m(i − 1)]n(i − 1) =<br />
d<br />
(n − d)n .<br />
Now let us look at the covariances among � φi <strong>and</strong> � φj (i �= i). It is very amazing that they<br />
are all approximately equal to zero!<br />
For example, let us consider the covariance between � φ1 <strong>and</strong> � φ2. Since � φ1 = log(�q1) <strong>and</strong><br />
�φ2 = log(�q2), using the same argument for the delta method, we know that we only need<br />
to find out the covariance between �q1 <strong>and</strong> �q2, or equivalently, the covariance between �m(0)<br />
<strong>and</strong> �m(1). This can be seen from the following:<br />
E[�m(0)�m(1)] = E[E[�m(0)�m(1)|n(0),d(0),w(0)]]<br />
= E[�m(0)E[�m(1)|n(0),d(0),w(0)]]<br />
= E[�m(0)m(1)]<br />
= m(1)E[�m(0)]<br />
= m(1)m(0) = E[�m(0)]E[�m(1)].<br />
Therefore, the covariance between �m(0) <strong>and</strong> �m(1) is zero. Similarly, we can show other<br />
covariances are zero. Hence,<br />
var(log( � S R (5))) = var( � φ1)+var( � φ2)+var( � φ3)+var( � φ4)+var( � φ5).<br />
Let � θ = log( � S R (5)). Then � S R (5) = e �θ .So<br />
var( � S R (5)) = (e θ ) 2 var(log( � S R (5))) = (S(5)) 2 [var( � φ1)+var( � φ2)+var( � φ3)+var( � φ4)+var( � φ5)],<br />
PAGE 16