SAS 语句中的itprint
选项导致显示迭代历史。这包括一个 Ridge 值,以及每次迭代的 beta 值和对数可能性。Ridge 通常为零,但当对数似然比前一次迭代的对数似然更负时,它就是非零。我需要知道 SAS 如何计算该脊值,并且在该过程的“详细信息”部分或其他任何地方都找不到任何内容。class
proc phreg
看来,默认情况下,Ridge 值始终为0.0001 * 2^n
,并且 SAS 以 开始n=0
并递增n
,直到对数似然比前一次迭代中的负数更小。但我已经测试了至少一个使用 SAS 的示例Ridge=0.4096
何时Ridge=0.2048
就足够了。
更新:我现在认为 SAS 正在迭代4^n
,而不是2^n
. 这解释了跳过2048
并且与我迄今为止的测试一致。
所以我想我已经回答了我自己的问题,现在希望这种方法得到学术支持。正如Robert Penridge和Joe建议的那样,我可能会在Cross Validated中寻求这一点。