r - 防止 NA 用于 lm 回归

Question

我有一个包含 future 返回的向量 Y 和一个包含当前返回的向量 X。最后一个 Y 元素是 NA，因为最后一个当前返回也是可用系列的最后。

X = { 0.1, 0.3, 0.2, 0.5 }
Y = { 0.3, 0.2, 0.5, NA }
Other = { 5500, 222, 523, 3677 }

lm(Y ~ X + Other)

我想确保每个系列的最后一个元素是 不是 包含在回归中。我阅读了 na.action 文档，但我不清楚这是否是默认行为。

对于 cor()，这是从计算中排除 X[4] 和 Y[4] 的正确解决方案吗？

cor(X, Y, use = "pairwise.complete.obs")

Answer 1

lm 的出厂默认值是忽略包含 NA 的观察结果值。由于这可以使用全局选项覆盖，您可能需要显式设置 na.action至 na.omit :

> summary(lm(Y ~ X + Other, na.action=na.omit))

Call:
lm(formula = Y ~ X + Other, na.action = na.omit)

[snip]

  (1 observation deleted due to missingness)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

关于你的第二个问题 cor(X,Y,use='pairwise.complete.obs')是正确的。由于只有两个变量， cor(X,Y,use='complete.obs')也会产生预期的结果。