
cluster()
in a formula implies that
robust sandwich variance estimators are desired.cluster(x)
x
coxph
, survreg
marginal.model <- coxph(Surv(time, status) ~ rx + cluster(litter), rats,
subset=(sex=='f'))
frailty.model <- coxph(Surv(time, status) ~ rx + frailty(litter), rats,
subset=(sex=='f'))
Run the code above in your browser using DataLab