考虑用例如 100 个数据点来拟合 coxph 模型。分析中仅包含 95 个，而 5 个因 NA(即缺失)而被排除在外。我提取了拟合数据的残差，所以我有一个包含 95 个观测值的残差向量。我想将残差包括回原始数据框中，但我不能这样做，因为长度不同。
如何确定模型中未包含原始数据框中的哪些观测值，以便排除/删除它们以使两个长度相同？
(原始数据大得多，所以很难定位数据丢失的地方......)
最佳答案
重新拟合您的模型，将 na.action 参数设置为 na.exclude。这会用 NA 填充作为拟合对象一部分的残差和拟合值。如果您的原始型号是 zn50:
zn50_na <- update(zn50, na.action=na.exclude)
这应该为您提供适当长度的 residuals(zn50_na) 和 fitted(zn50_na)。有关详细信息，请参阅 ?na.omit。
关于r - 在R中找到拟合模型中的缺失，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/17428950/