中的survival
软件包R
似乎专注于连续时间生存模型。我有兴趣估计比例风险模型的离散时间版本,即互补对数模型。我有一个相当简单的生存模型,带有简单的右审查。
我知道估计此模型的一种方法是创建一个数据集,该数据集对于每个未“死”的时期的每个观察值都有一个单独的行。然后,可以使用glm
带有链接的模型。cloglog
这种方法看起来内存效率很低;实际上,它可能会产生一个对于我机器上的内存来说太大的数据集。
第二种方法是自己编写 MLE。这很简单,但我希望有一个包装可以装上这种生存模型。使用包只会更容易协作并避免编码错误。
有人知道这样的包吗?