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。