Hi!
In answer to:
--------
I noticed that you were using what might be called an "externally
created Surv object". I have a memory that Terry Therneau has
criticized that practice. I cannot remember if it was in exactly this
situation but I might ask if setting up the model as:
cox = coxph(Surv(stime, event) ~ bucket*(today + accor + both) +
activity, data = data)
... might give the survival machinery a better handle on where
everything might be found.
------------
I tried to create the Surv object "internally" but I face the same issue:
> (cox.s = coxph(Surv(time=absence, event=(censored==FALSE)) ~
> bucket*(today) + strata(activity), data = small))
Call:
coxph(formula = Surv(time = absence, event = (censored == FALSE)) ~
bucket * (today) + strata(activity), data = small)
coef exp(coef) se(coef) z p
bucket575 0.4526 1.572 0.740 0.612 0.54
todayTRUE -0.0886 0.915 0.676 -0.131 0.90
bucket575:todayTRUE -0.1670 0.846 0.794 -0.210 0.83
Likelihood ratio test=2.32 on 3 df, p=0.509 n= 100, number of events= 100
> fit = survfit(cox.s, newdata=small[1:50,])
Error in model.frame.default(data = small[1:50, ], formula = ~bucket + :
number of variables != number of variable names
Best, and thank you for the suggestion.
ge
--
View this message in context:
http://r.789695.n4.nabble.com/survfit-number-of-variables-number-of-variable-names-tp4649834p4650080.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org 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.
