Timeline for Cox baseline hazard
Current License: CC BY-SA 3.0
12 events
when toggle format | what | by | license | comment | |
---|---|---|---|---|---|
Mar 21, 2020 at 13:07 | comment | added | mhy | As \@ocram defined the baseline hazard, you can find this formula in the book Applied Survival Analysis Using R-2016 by Dirk F. Moore (page 64), either. I applied the formula for my dataset and it works, in fact, it estimated the baseline hazrad correctly as I checked it couple of times for some different covariates. And also, the slight difference in the results because of method is used. As @mr.bjerre stated, we will get the same results when we choose the method='breslow' in fitting cox proportional hazard. | |
May 16, 2018 at 7:16 | comment | added | mr.bjerre |
As @tomka mentioned. Replacing the coxph call with fit<-coxph(Surv(time, status)~age, data=kidney, method="breslow") will fix the difference in methods.
|
|
Mar 22, 2017 at 20:45 | history | edited | gung - Reinstate Monica | CC BY-SA 3.0 |
shortened code
|
Feb 28, 2017 at 18:10 | comment | added | tomka |
I need to note some problems with this implementation. Using kidney$time >= y[l] can run into numerical problems when time is numeric due to the tabulation in creating $y$. Furthermore, the way you define your risk set is inaccurate, because if there is a tie of two observations, one with status=0 and one with status=1 , then $d=2$ but your code gives $d=1$ as you exclude all status=0 observations. The latter problem applies higher numbers of ties likewise.
|
|
Feb 27, 2017 at 17:49 | comment | added | tomka | Method of Breslow (1974) | |
Dec 26, 2012 at 11:32 | comment | added | ocram | The discretised hazard is zero, except at event times. This indeed gives the largest contribution to the likelihood if a discrete hazard function is supposed. You might want to interpolate between any two estimates assuming, for example, that the hazard stays constant. | |
Dec 26, 2012 at 11:11 | comment | added | Dihan | Thanks a lot. Yes, there are slight difference for approximation method. But there are 76 time points with ties, if I want to find the baseline hazard for every time point. What can i do? What type of modification in R code is needed? | |
Dec 26, 2012 at 5:38 | vote | accept | Dihan | ||
Dec 25, 2012 at 19:31 | history | edited | ocram | CC BY-SA 3.0 |
added 1 characters in body
|
Dec 25, 2012 at 19:25 | history | edited | ocram | CC BY-SA 3.0 |
added 99 characters in body
|
Dec 25, 2012 at 19:13 | history | edited | ocram | CC BY-SA 3.0 |
added 132 characters in body
|
Dec 25, 2012 at 19:07 | history | answered | ocram | CC BY-SA 3.0 |