r - 在 R 中找到拟合模型中的缺失

Question

考虑coxph使用 100 个数据点拟合模型。分析中仅包括 95 个，而由于存在NA（即缺失）而排除了 5 个。我提取拟合数据的残差，因此我有一个包含 95 个观测值的残差向量。我想将残差包含回原始数据框中，但我不能这样做，因为长度不同。

如何识别原始数据框中的哪些观察未包含在模型中，以便我可以排除/删除它们以使两个长度相同？

（原始数据要大得多，因此很难找到数据丢失的地方......）

Answer 1

重新拟合您的模型，将na.action参数设置为na.exclude. NA这用s填充作为拟合对象一部分的残差和拟合值。如果您的原始模型是zn50：

zn50_na <- update(zn50, na.action=na.exclude)

这应该给你residuals(zn50_na)和fitted(zn50_na)适当的长度。有关?na.omit更多信息，请参阅。