12.12.2012 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!