Question: extension of Firth's logistic regression to cox regression model
0
gravatar for ycding
11 months ago by
ycding0
ycding0 wrote:

for a 2 by 2 table, if one entry has value of zero, we can run Firth's Logistic regression using "logistf" R package. what about the outcome or dependent variable is survival time, is there a corresponding R package for running Cox regression model for rare variants or copy number alterations?

thanks you,

Ding

R • 501 views
ADD COMMENTlink modified 11 months ago • written 11 months ago by ycding0
0
gravatar for Kevin Blighe
11 months ago by
Kevin Blighe35k
Republic of Ireland
Kevin Blighe35k wrote:

I believe that you are looking for the Cox proportional hazards regression model:

By a strange coincidence of statistics, the log likelihoods of a conditional logistic regression model are also equivalent to those of the Cox model. So, you can additionally fit a conditional logistic regression model if you have a matching stratum:

Kevin

ADD COMMENTlink modified 11 months ago • written 11 months ago by Kevin Blighe35k

no, I worry that proportional hazards regression model does not fit well to data set with rare variants as predictor variables. for example, for a marker or gene, 3 of 100 patients carry a mutations (1) and all those three patients show long survival time (> 5 years, with no recurrence) compared to the other 97 patients (< 5 years) who do not carry the specific mutation. if considering survival time as a categorical dependent variable, 0 for < five years and 1 for > five years, one would use Firth's logistic regression model instead of regular logistic regression model for those sparse data. However, If one wants to consider survival time as continuous variable and fit a Cox proportional hazard regression model, the hazard ratio would be large or estimated inaccurately. I would think someone has already develop a modified Cox regression model to deal with those data set with rare variant variable as independent variable.

ADD REPLYlink written 11 months ago by ycding0
0
gravatar for ycding
11 months ago by
ycding0
ycding0 wrote:

OK, "coxphf" R package is what I was looking for, which Implements Firth’s penalized maximum likelihood bias reduction method for Cox regression.It has been shown to provide a solution in case of monotone likelihood (nonconvergence of likelihood function). The program fits profile penalized likelihood confidence intervals which were proved to outperform Wald confidence intervals.

ADD COMMENTlink modified 11 months ago • written 11 months ago by ycding0

Great - I will note this for future myself.

ADD REPLYlink written 11 months ago by Kevin Blighe35k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 639 users visited in the last hour