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
Hi,
the package 'survival' provides the possibility to estimate Cox PH
models with a frailty term.
library(survival)
?frailty
If this is not what you are looking for, there is also a package called
'frailtypack'. But I have no experience with that.
I hope this helps.
Best,
Roland
P.S. Boccon
Dear R-help,
I am trying to estimate a Cox model with nested effects, or better
h(t,v,w)=v*w*h0(t)*exp(B'x)
where h(t,v,w) is the individual hazard function
w and v are both frailty terms (gamma or normal distributed)
I have 12 clusters and for each one of them I would like to associate a
3 matches
Mail list logo
