Displaying 1 result from an estimated 1 matches for "brachycera".
2009 Aug 10
0
survival:: plotting survfit with two predictors
...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<-coxph(Surv(a$deg.below35,a$censor)~a$morder+a$tmean);
summary(model2) #is fine
plot(survfit(model2,conf.type=&q...