Displaying 2 results from an estimated 2 matches for "hymenoptera".
2009 Aug 10
0
survival:: plotting survfit with two predictors
...derbane
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="none",newdata=d...
2006 May 17
0
script S-plus -> R
...argumentos adicionais n?o pareados s?o desconsiderados in:
density.default(tdon21, kernel = "gaussian", bandwidth = h[j],
I'm sorry for so lengthy e-mail, so here are my apologies. But I don't
know what I'm doing wrong.
Thanks,
Rog?rio
--
Rog?rio R. Silva
Laborat?rio de Hymenoptera
Museu de Zoologia da USP
Av. Nazar? 481, 04263-000 S?o Paulo SP
GNU/Linux User # 354364
Linux-Debian Etch: 2.6.15