Displaying 1 result from an estimated 1 matches for "fdt_weath_dttran".
Did you mean:
fdt_weath_dttrans
2009 Aug 10
0
survival:: plotting survfit with two predictors
...suspect it is not a bug but
rather a problem with my code. I have not been 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...