How do I fit a mixed effects model with two crossed random effects for grouped
time survival data?
I tried coxme with no luck.
Suppose that y is survival time, uncens is censoring indicator, trt is
treatment below.
ran.eff1 and ran.eff2 below are two crossed random effects. This way, it
wouldn't run. If I remove one of the random effects, it runs perfectly.
coxme(Surv(y,uncens) ~ trt + (1|ran.eff1)+(1|ran.eff2),data=mydataframe)
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
