Hello R users
Would somebody know how to estimate survival from a frailty cox model,
using the function survfit
and the argument newdata ? (or from any other way that could provide
individual expected survival
with standard error); Is the problem related to how the random term is
included in newdata ?
kfitm1 <- coxph(Surv(time,status) ~ age + sex + disease + frailty(id,
dist='gauss'), kidney)
survfit(kfitm1) #ok for mean expected survival
Call: survfit.coxph(object = kfitm1)
n events median 0.95LCL 0.95UCL
76 58 63 39 132
survfit(kfitm1, newdata=kidney[1,]) #return an error message
Erreur dans x2 %*% coef : arguments inadéquats
Thanks in Advance
--
Emmanuel Tillard
Veterinaire
CIRAD-EMVT
Unite de recherche 18
UMR868 Elevage des Ruminants en Regions Chaudes (ERRC)
Campus ENSA-INRA
2 place Viala
34060 Montpellier cedex 1
tel: 0499612265 (fixe)
0633850598 (gsm)
fax: 0467545694
e-mail: [EMAIL PROTECTED]
______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
