You want to use coxme, currently found in the kinship package. It uses lme
style input, e.g.
coxme(Surv(time, status) ~ age + weight, data=mydata,
random= ~1 | school/teacher)
One note: your message seems to imply that you want to use group/subject as
the random effect. A Cox model with one random effect per subject is not
stable
-- actually an one where df(random effect) approaches the number of events is
not --- and you will likely get a variance estimate whose confidence interval
is
(0, huge number).
Terry Therneau
______________________________________________
[email protected] 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.
