Displaying 1 result from an estimated 1 matches for "below35".
Did you mean:
below
2009 Aug 10
0
survival:: plotting survfit with two predictors
...n able to find any other
solutions in the achieve and I would greatly appreciate any suggestions.
Thanks in advance,
Marshall Knoderbane
Below is my code for the two models and the error message:
a<-read.delim("FDT_weath_dttrans.txt")
library(survival)
model1<-coxph(Surv(a$deg.below35,a$censor)~a$morder);
plot(survfit(model1,conf.type="none",newdata=data.frame(morder=c("Brachycera","Hymenoptera","Nematocera","OtherCol","Zygoptera"))),lty=1:5)
#produces a fine plot with a separate curve for each level of morder
model2<...