[R] Plot survival curves after coxph() with frailty() random effects terms
fromnorden at gmail.com
Fri Jun 23 16:17:39 CEST 2017
I would like to plot a survival curves of a group with different categories
after running a Cox model with frailty() random effects terms.
I just could display a survival plot of the covariable’s mean.
Here an example:
fit<-coxph(Surv(time, status) ~ sex+ frailty(litter, dist='gamma',
plot(suf, xscale=1, xlab = "Days", ylab="Survival", conf.int=FALSE)
In `[.survfit`(suf, 1) : survfit object has only a single survival curve
#But if I use the next code I get the 2 groups.
surv_group <- survfit(Surv(time, status) ~ sex, rats)
lines(surv_group[1:2], lwd=2, conf.int=FALSE)
However, I am not sure about these 2 curves are well done, appropriate.
If any of you could help…
Andreu Ferrero Gregori
[[alternative HTML version deleted]]
More information about the R-help