R coxph()中利用 frailty 创建的 混合效应模型,如下
coxph(Surv(Time_range,
Events_Deceased) ~ Avg_Thickness_FOV + IOP +
HbA1c + ratio + Hypertension + Age+
frailty(ID, distribution = "gaussian",
sparse = FALSE, method = "reml"),
ads)结果和coxme 得到的不一致,代码如下
coxme(Surv(Time_range, Events_Deceased) ~ Avg_Thickness_FOV + IOP +
HbA1c + ratio + Hypertension + Age + (1|ID),
ads)
`