1

我正在解决一个问题，我想使用其他两个线性模型的残差来构建线性模型。我使用 UN3 数据集来显示我的问题，因为它比使用我的实际数据集更容易将问题放在这里。

这是我的 R 代码：

head(UN3)
m1.lgFert.purban <- lm(log(Fertility) ~ Purban, data=UN3)
m2.lgPPgdp.purban <- lm(log(PPgdp) ~ Purban,  data=UN3)
m3 <- lm(residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban))

这是我得到的错误：

> m3 <- lm(residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban))
Error in model.frame.default(formula = residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban),  : 
  variable lengths differ (found for 'residuals(m2.lgPPgdp.purban)')

我真的不明白为什么会发生这个错误。如果是与日志相关的问题，那么我在构建前两个模型时应该会收到错误消息。

4

1 回答 1

6

na.action很可能是您的默认设置na.omit（用 勾选options("na.action")）。这意味着NA值会被无声地删除，从而导致残差向量的长度不同。您可能想要使用，它用sna.action="na.exclude"填充残差。NA

library(alr3)
options("na.action")
#$na.action
#[1] "na.omit"

m1.lgFert.purban <- lm(log(Fertility) ~ Purban, data=UN3,na.action="na.exclude")
m2.lgPPgdp.purban <- lm(log(PPgdp) ~ Purban,  data=UN3,na.action="na.exclude")

m3 <- lm(residuals(m1.lgFert.purban) ~ residuals(m2.lgPPgdp.purban))
#Coefficients:
#                 (Intercept)  residuals(m2.lgPPgdp.purban)  
#                    -0.01245                      -0.18127
于 2013-02-17T18:40:22.753 回答