A Markov chain Monte Carlo EM algorithm for analyzing interval-censored data under the Cox proportional hazards model.
This paper proposes a Monte Carlo EM (MCEM) algorithm for fitting the proportional hazards model for interval-censored failure-time data. The algorithm generates orderings of the failures from their probability distribution under the model. We maximize the average of the log-likelihoods from these completed data sets to obtain updated parameter estimates. As with the standard Cox model, this algorithm does not require the estimation of the baseline hazard function. The performance of the algorithm is evaluated using simulations, and the method is applied to data from AIDS and cancer studies. Our results indicate that our method produced more precise and unbiased estimates than methods of right and midpoint imputation.
Goggins
W B
WB
Biostatistics Department, Harvard School of Public Health, Boston, Massachussetts 02115, USA.
Finkelstein
D M
DM
Schoenfeld
D A
DA
Zaslavsky
A M
AM
D000163:Acquired Immunodeficiency Syndrome / Q000188:drug therapy; D000465:Algorithms; D000704:Analysis of Variance; D019380:Anti-HIV Agents / Q000627:therapeutic use; D001699:Biometry / Q000379:methods*; D001943:Breast Neoplasms / Q000628:therapy; D003198:Computer Simulation; D004352:Drug Resistance, Microbial; D005260:Female; D015658:HIV Infections / Q000635:transmission; D006467:Hemophilia A / Q000628:therapy; D006801:Humans; D016013:Likelihood Functions; D008390:Markov Chains; D009010:Monte Carlo Method; D016016:Proportional Hazards Models; D065227:Transfusion Reaction; D015215:Zidovudine / Q000627:therapeutic use
D019380:Anti-HIV Agents; D015215:Zidovudine