Abstract
Relative risk regression models can be fitted using a log-link binomial GLM, however standard algorithms can suffer convergence problems. Combinatorial EM (CEM) algorithms that provide stable convergence can be computationally intensive, particularly for large models. We present a new approach using an EM algorithm with an overparameterised model that retains the stability of the CEM algorithm but greatly reduces computing time. This is demonstrated with a small example in which modified Fisher scoring fails to converge to the MLE, and a bootstrap analysis of data from a clinical trial in heart attack patients.
Original language | English |
---|---|
Title of host publication | Proceedings of the 31st International Workshop on Statistical Modelling |
Editors | Jean-François Dupuy, Julie Josse |
Place of Publication | Amsterdam, the Netherlands |
Publisher | Statistical Modelling Society |
Pages | 93-98 |
Number of pages | 6 |
Volume | 1 |
Publication status | Published - 2016 |
Event | International Workshop on Statistical Modelling (31st : 2016) - Rennes, France Duration: 4 Jul 2016 → 8 Jul 2016 Conference number: 31st |
Conference
Conference | International Workshop on Statistical Modelling (31st : 2016) |
---|---|
Abbreviated title | IWSM2016 |
Country/Territory | France |
City | Rennes |
Period | 4/07/16 → 8/07/16 |
Keywords
- Binomial regression
- EM algorithm
- Relative risk