Premium
Bayesian proportional hazards model with time‐varying regression coefficients: a penalized Poisson regression approach
Author(s) -
Lambert Philippe,
Eilers Paul H. C.
Publication year - 2005
Publication title -
statistics in medicine
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.996
H-Index - 183
eISSN - 1097-0258
pISSN - 0277-6715
DOI - 10.1002/sim.2396
Subject(s) - mathematics , covariate , poisson distribution , poisson regression , statistics , proportional hazards model , regression analysis , bayesian probability , bayesian inference , linear regression , parametric statistics , econometrics , population , demography , sociology
One can fruitfully approach survival problems without covariates in an actuarial way. In narrow time bins, the number of people at risk is counted together with the number of events. The relationship between time and probability of an event can then be estimated with a parametric or semi‐parametric model. The number of events observed in each bin is described using a Poisson distribution with the log mean specified using a flexible penalized B‐splines model with a large number of equidistant knots. Regression on pertinent covariates can easily be performed using the same log‐linear model, leading to the classical proportional hazard model. We propose to extend that model by allowing the regression coefficients to vary in a smooth way with time. Penalized B‐splines models will be proposed for each of these coefficients. We show how the regression parameters and the penalty weights can be estimated efficiently using Bayesian inference tools based on the Metropolis‐adjusted Langevin algorithm. Copyright © 2005 John Wiley & Sons, Ltd.