【发布时间】：2014-07-25 13:10:25
【问题描述】：
SAS proc phreg 的class 语句中的itprint 选项导致显示迭代历史记录。这包括一个 Ridge 值，以及每次迭代的 beta 值和对数可能性。 Ridge 通常为零，但当对数似然比前一次迭代的对数似然更负时，它就是非零。我需要知道 SAS 如何计算该脊值，但我在该过程的详细信息部分或其他任何地方都找不到任何内容。
看来，默认情况下，Ridge 值始终为 0.0001 * 2^n，并且 SAS 以 n=0 开头并递增 n，直到对数似然度低于上一次迭代时的负数。但我已经测试了至少一个例子，其中 SAS 使用 Ridge=0.4096 而 Ridge=0.2048 就足够了。
更新：我现在认为 SAS 正在迭代 4^n，而不是 2^n。这解释了跳过 2048 并且与我迄今为止的测试一致。
所以我想我已经回答了我自己的问题，现在希望这种方法得到学术支持。我很可能会在Cross Validated 上寻求Robert Penridge 和Joe 的建议。
【问题讨论】：
-
如果您在这里没有成功，您应该考虑将问题标记为将其迁移到交叉验证 (stats.stackexchange.com)。
-
我认为这不是一个合适的简历问题；这是询问特定软件如何执行计算，而不是通常如何执行计算。当然不要交叉发布它，如果您正在寻找更通用的统计建议，则应该标记迁移。
-
也就是说，如果您要将其发布到其他地方，我建议您使用SAS Community forums。
-
谢谢你的建议，乔。我最终追求的是支持 SAS 使用的任何方法的学术论文或文本。我认为将 SAS 的显示值与我自己的 VB 例程中的计算进行比较是一个不错的起点。
-
我认为，对这种方法的学术支持对 CV 来说很好。最初的问题（SAS 是如何...）不会是，但是一旦您进入不再真正依赖特定软件来回答的领域，就可以了。
标签： sas survival-analysis