假设我有一个“肾导管”数据集。数据是关于使用便携式透析设备的肾病患者在插入导管时感染的复发时间。导管可能会因感染以外的原因被移除，在这种情况下，观察会受到审查。每个患者恰好有 2 个观察值。

如果我想使用 EM 算法拟合具有随机效应（此处为 Gamma 脆弱）的 Cox PH 模型。通过使用内置的生存包 R 代码coxph()，我可以轻松地做到这一点

library(survival)
data(kidney)
fit<-coxph(Surv(time, status) ~ age + sex + 
           frailty(id, dist='gamma', method='em'), kidney)

但是如果我想为 EM 算法编写一个逐步函数，我该如何进行呢？