R中的离散时间危害模型（博客）

该survival软件包R似乎专注于连续时间生存模型。我对估计比例风险模型（互补对数-对数模型）的离散时间版本感兴趣。我有一个非常简单的生存模型，并带有简单的权限检查。

我知道，估算此模型的一种方法是创建一个数据集，该数据集在不是“死角”的每个时期的每个观察值中都有单独的一行。然后，可以使用glm带有cloglog链接的模型。

这种方法似乎对内存效率很低。实际上，它可能会产生一个数据集，该数据集对于我的计算机上的内存而言太大。

第二种方法是自己编写MLE。那将足够简单，但是我希望有一个可以固定这种生存模型的软件包。协作更容易，并且避免使用程序包时出现编码错误。

有人知道这样的包裹吗？

每个观察有几行似乎是多余的，但可能并非如此。如果模型中存在任何随时间变化的协变量，则每个观察月肯定会需要其自己的行。时变协变量的一个特定示例是经过的时间。由于该变量几乎肯定应包含在模型中，因此在每个观察期有单独的一行是有意义的。因此，建议的第一种方法可能是最好的方法。

请注意，这不同于具有威布尔分布的连续时间比例风险模型。在这里，如果经过的时间是唯一随时间变化的协变量，则可以将每次观察的生存模型简化为一条直线（例如，参见此处）。Cox比例风险模型也有类似的结果。