This preview shows pages 1–4. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: Solution to Homework Assignment 4 1 Question 1 a. Let U i = min( T i ,C i ) be the observed time, i = I ( T i C i ) be the failure indicator and n ( U i ) be the number of patients at risk at time U i for the ith individual, where T i and C i are the underlying survival and censoring times, respectively, i = 1 ,..., 12. If combing the two treatments together, we get the ordered time U i listed in the left part of the following table. The KaplanMeier estimator KM ( t ) and its standard error seKM ( t ), for any time t 0, are described in the right part of the following table, where KM ( t ) = Y i : U i t 1 1 n ( U i ) i and se { KM ( t ) } = KM ( t ) " X i : U i t i n ( U i ) { n ( U i ) 1 } # 1 / 2 . Table 1: KaplanMeier estimator of the survival curve and its standard error U i i n ( U i ) t KM ( t ) se { KM ( t ) } 0.73 1 12 [0.00, 0.73) 1.000 0.000 0.95 11 [0.73, 1.43) 0.917 0.080 1.43 1 10 [1.43, 1.58) 0.825 0.113 1.58 1 9 [1.58, 1.98) 0.733 0.132 1.98 1 8 [1.98, 2.25) 0.642 0.144 2.25 1 7 [2.25, 2.30) 0.550 0.150 2.30 1 6 [2.30, 3.24) 0.458 0.150 3.24 1 5 [3.24, 3.30) 0.367 0.146 3.30 1 4 [3.30,+ ) 0.275 0.135 3.89 3 4.03 2 4.20 1 b. Let S 1 ( t ) = P ( T t  trt = A ) and S 2 ( t ) = P ( T t  trt = B ). To test the null hypothesis H : S 1 ( t ) = S 2 ( t ), t 0, we consider the logrank test T = i : i =1 n d 1 ( u i ) d ( u i ) n ( u i ) n 1 ( u i ) o q i : i =1 n 1 ( u i ) n 2 ( u i ) d ( u i ) { n ( u i ) d ( u i ) } n 2 ( u i ) { n ( u i ) 1 } , where d 1 ( u i ) is the number of deaths at time u i from treatment A, d ( u i ) is the total number of deaths from both treatments, n 1 ( u i ), n 2 ( u i ) and n ( u i ) are the number of patients at risk at time u i for treatment A, B and the both treatments, respectively. 1 Figure 1: KaplanMeier estimator of the survival curve 1 2 3 4 0.4 0.6 0.8 1.0 Time in years S(x) From the provided survival data, we see that a total of 8 deaths were observed, but no more than 1 death occurred at the same time. Table 2: Death times u i : d ( u i ) = 1 0.73 1.43 1.58 1.98 2.25 2.30 3.24 3.30 d 1 ( u i ) 1 1 1 1 n 1 ( u i ) 6 4 3 3 3 2 2 1 n ( u i ) 12 10 9 8 7 6 5 4 According to the summary information of the 8 death times as described in the above table 2, we have the estimated logrank test T = i : i =1 n d 1 ( u i ) n 1 ( u i ) n ( u i ) o q i : i =1 n 1 ( u i ) n 2 ( u i ) n 2 ( u i ) = . 9798 1 . 8412 = 0 . 722 , and the corresponding pvalue for a twosided test is 0.4703. 2 2 Question 2 Let and 1 be the hazards of the standard and new treatments, respectively. Since the survival distribution is approximated by an exponential distribution, we have that m m 1 = 1 = exp( A ) , where m and m 1 are the median survival times for these two treatments, respectively, which gives A = log(4 / 5) = . 2231....
View
Full
Document
This note was uploaded on 01/15/2012 for the course PHC 6020 taught by Professor Staff during the Fall '11 term at University of Florida.
 Fall '11
 Staff

Click to edit the document details