Bias-reduced and separation-proof conditional logistic regression with small or sparse data sets.
Conditional logistic regression is used for the analysis of binary outcomes when subjects are stratified into several subsets, e.g. matched pairs or blocks. Log odds ratio estimates are usually found by maximizing the conditional likelihood. This approach eliminates all strata-specific parameters by conditioning on the number of events within each stratum. However, in the analyses of both an animal experiment and a lung cancer case-control study, conditional maximum likelihood (CML) resulted in infinite odds ratio estimates and monotone likelihood. Estimation can be improved by using Cytel Inc.'s well-known LogXact software, which provides a median unbiased estimate and exact or mid-p confidence intervals. Here, we suggest and outline point and interval estimation based on maximization of a penalized conditional likelihood in the spirit of Firth's (Biometrika 1993; 80:27-38) bias correction method (CFL). We present comparative analyses of both studies, demonstrating some advantages of CFL over competitors. We report on a small-sample simulation study where CFL log odds ratio estimates were almost unbiased, whereas LogXact estimates showed some bias and CML estimates exhibited serious bias. Confidence intervals and tests based on the penalized conditional likelihood had close-to-nominal coverage rates and yielded highest power among all methods compared, respectively. Therefore, we propose CFL as an attractive solution to the stratified analysis of binary data, irrespective of the occurrence of monotone likelihood. A SAS program implementing CFL is available at: http://www.muw.ac.at/msi/biometrie/programs.
Heinze
Georg
G
Section of Clinical Biometrics, Core Unit of Medical Statistics and Informatics, Medical University of Vienna, Spitalgasse 23, Vienna A-1090, Austria. georg.heinze@meduniwien.ac.at
Puhr
Rainer
R
D000783:Aneurysm / Q000453:epidemiology; D000818:Animals; D015982:Bias; D056808:Biostatistics; D001943:Breast Neoplasms / Q000532:radiotherapy; D016022:Case-Control Studies; D003198:Computer Simulation / Q000706:statistics & numerical data; D015987:Effect Modifier, Epidemiologic; D005260:Female; D006493:Heparin / Q000009:adverse effects / Q000627:therapeutic use; D006801:Humans; D016013:Likelihood Functions; D016015:Logistic Models; D008175:Lung Neoplasms / Q000453:epidemiology / Q000209:etiology; D009381:Neoplasms, Radiation-Induced / Q000453:epidemiology; D051381:Rats; D012307:Risk Factors; D012907:Smoking / Q000009:adverse effects / Q000453:epidemiology; D012984:Software / Q000706:statistics & numerical data; D014183:Transplantation, Heterologous / Q000009:adverse effects / Q000706:statistics & numerical data
D006493:Heparin