在R中定位拟合模型中的缺失值。

Question

3

考虑使用coxph模型进行拟合，例如使用100个数据点。在分析过程中，只有95个数据点被包含在内，而5个数据点由于是NA（即缺失值）而被排除。我从拟合数据中提取残差，因此我得到一个包含95个观测值的残差向量。我想将这些残差添加回原始数据框中，但是由于长度不同，我无法做到这一点。

如何确定原始数据框中未被纳入模型的观测值，以便我可以排除/删除它们，使两个长度相同？

（原始数据很大，因此很难确定数据缺失的位置...）

- user2543095

也许是 na.omit，但不确定你想要什么。通常，您需要提供一个最小工作示例。 - Tyler Rinker

我知道如何查找特定列中缺失值的行，但我不知道coxph方法如何将某些内容分类为“缺失”。 - user2543095

1

你是自己排除了这5个值，还是coxph采取了某些行动？如果是前者，你知道你排除（或包含）的值的索引，因此将观察映射到该索引集合上。 - Carl Witthoft

谢谢，我会选择后者！ - user2543095

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Hong Ooi · Accepted Answer

重新调整您的模型，将 na.action 参数设置为 na.exclude。这将在拟合对象的残差和拟合值中填充 NA。如果您的原始模型是 zn50：

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

这应该给你适当长度的残差（zn50_na）和拟合值（zn50_na）。有关更多信息，请参见?na.omit。