R:用于审查观察的Survreg()包对数似然代码?

时间:2013-10-04 15:40:52

标签: r

我想在survreg.fit()中找到代码的特定部分,其中定义了例如左或右删失观察的密度。让我们说我有响应变量:

Surv(c(1,1,2,3),c(1,2,2,3),c(1,0,2,1),type="interval")

第一次和最后一次观察是按时发生的事件,而第二次和第三次观察分别是右侧和左侧观察。我当然知道你可以在书籍和其他地方找到对数似然的公式 - 我有兴趣查看实际的代码。那么,survreg.fit()中代码的哪一部分处理左边的删失数据呢?

我显然试图自己解决这个问题,但是没有编写这种复杂性的对数似然代码(即,许多不同类型的审查,分布,缩放因子等),我只是迷失了码。特别是考虑到我打电话survreg()我打电话survreg.fit()。请注意,我并不批评 - 我只是想学习!

干杯,

尼古拉

0 个答案:

没有答案