r - coxph 对象中的脆弱性估计

Question

如果使用 obj=coxph(... + frailty(id) )，则该对象还返回每个个体的 (log)frailty 估计值，可以使用 obj$frail 提取。

有谁知道这些估计是如何获得的？它们是经验贝叶斯估计吗？

谢谢！

西奥多

Answer 1

frailty可以在页面中看到默认分布为?frailty“gamma”。如果您查看 frailty 函数（未隐藏），您会发现它只是将分布的名称粘贴到“frailty”上。并用于get()检索正确的功能。因此，请查看frailty.gamma（也未隐藏）以找到您问题的答案。再次查看帮助页面，您可以看到我应该能够在不查看代码的情况下弄清楚所有这些，因为它就在页面的顶部。但是 R 有很多获取知识的途径。（估计是 ML，而不是“经验贝叶斯”。）

帮助页面向我建议，作者 (Therneau) 希望您咨询 Therneau 和 Grambsch，以获取阅读代码时不明显的更多细节。如果您正在认真研究 R 中的生存模型，那么这是一本非常有用的书。这对于理解“生存”包的基础非常清晰和有帮助。