首页 | 本学科首页   官方微博 | 高级检索  
     


Free Knot Splines with RJMCMC in Survival Data Analysis
Authors:M. Denis  N. Molinari
Affiliation:1. Laboratoire de Biostatistique, Institut Universitaire de Recherche Clinique , Université Montpellier 1 , Montpellier, France marie.denis@inserm.fr;3. Laboratoire de Biostatistique, Institut Universitaire de Recherche Clinique , Université Montpellier 1 , Montpellier, France;4. Service de Biostatistique, Epidémiologie, Santé Publique et Information Médicale (BESPIM) , Hopital Caremeau, CHU N?mes , N?mes, France
Abstract:The B-spline representation is a common tool to improve the fitting of smooth nonlinear functions, it offers a fitting as a piecewise polynomial. The regions that define the pieces are separated by a sequence of knots. The main difficulty in this type of modeling is the choice of the number and the locations of these knots. The Reversible Jump Markov Chain Monte Carlo (RJMCMC) algorithm provides a solution to simultaneously select these two parameters by considering the knots as free parameters. This algorithm belongs to the MCMC techniques that allow simulations from target distributions on spaces of varying dimension. The aim of the present investigation is to use this algorithm in the framework of the analysis of survival time, for the Cox model in particular. In fact, the relation between the hazard ratio function and the covariates being assumed to be log-linear, this assumption is too restrictive. Thus, we propose to use the RJMCMC algorithm to model the log hazard ratio function by a B-spline representation with an unknown number of knots at unknown locations. This method is illustrated with two real data sets: the Stanford heart transplant data and lung cancer survival data. Another application of the RJMCMC is selecting the significant covariates, and a simulation study is performed.
Keywords:B-splines  Cox model  Reversible Jump Markov Chain Monte Carlo method  Splines
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号