R中的间隔审查Cox比例风险模型

给定间隔审查的生存时间，我如何在中执行间隔审查的Cox PH模型R？随便搜索即可找到该软件包intcox，该软件包不再存在于R存储库中。我几乎肯定该包中的coxph功能survival无法处理间隔检查的生存数据。

另外，我不想插入数据然后使用该coxph函数。该方法低估了系数的标准误差，因为您忽略了间隔检查的不确定性。

如上所述，您可以使用survreg函数。不过请注意：严格来说，这不是Cox PH模型，而是位置比例模型。使用默认的日志转换，这是船尾模型。在指数分布的情况下，比例风险和后部模型是等效的，因此，如果将分布设置为指数，则这是具有指数基线的比例风险模型。同样，如果使用基线威布尔分布后部模型，则参数估计值只是具有威布尔基线分布的比例风险模型中所使用的线性估计。但总的来说，survreg不适合Cox PH模型。

如果需要一个半参数模型（如intcox中实现的模型），请当心：在当前版本的intcox中存在几个问题（算法通常会在距MLE较远的地方终止，在未经审查的情况下彻底失败，没有标准错误）自动呈现）。

可以使用的新替代方法是软件包“ icenReg”。

接纳偏见：这是icenReg的作者。

若要在R中进行间隔审查分析，必须创建一个Surv对象，然后使用survfit（）。如果您有多个变量，则intcox包可以解决该问题。