在coxph对象中的脆弱估计

Question

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

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

谢谢！

特奥多尔

Answer 1

frailty页面的默认分布在?frailty页面中可以看作“gamma”。如果你看一下脆弱函数（没有隐藏），你会发现它只是将分布的名称粘贴到“脆弱”上。并使用get()来检索正确的函数。因此，请查看frailty.gamma（也未隐藏）以查找问题的答案。再次回顾帮助页面，您可以看到我应该能够在不查看代码的情况下解决所有问题，因为它位于页面顶部。但是R有许多知识途径（它们是ML，而不是“经验贝叶斯”，估计。）

帮助页面告诉我，作者（Therneau）希望您咨询Therneau和Grambsch，以获取阅读代码时不明显的更多详细信息。如果你正在认真研究R中的生存模型，这是一本非常有用的书。它非常清楚，有助于理解“生存”包的基础。