10

在 R 中，当使用 时lm()，如果我na.action = na.pass在对 的调用中设置lm()，那么在汇总表中，对于任何无法估计的系数都有一个 NA（因为在这种情况下缺少单元格）。

但是，如果我使用summary(myModel)$coefficients或仅从汇总对象中提取系数coef(summary(myModel))，则省略 NA。

我希望在提取系数时包含 NA，就像在打印摘要时包含它们一样。有没有办法做到这一点？

设置options(na.action = na.pass)似乎没有帮助。

这是一个例子：

> set.seed(534)
> myGroup1 <- factor(c("a","a","a","a","b","b"))
> myGroup2 <- factor(c("first","second","first","second","first","first"))
> myDepVar <- rnorm(6, 0, 1)
> myModel <- lm(myDepVar ~ myGroup1 + myGroup2 + myGroup1:myGroup2)
> summary(myModel)

Call:
lm(formula = myDepVar ~ myGroup1 + myGroup2 + myGroup1:myGroup2)

Residuals:
       1        2        3        4        5        6 
-0.05813  0.55323  0.05813 -0.55323 -0.12192  0.12192 

Coefficients: (1 not defined because of singularities)
                    Estimate Std. Error t value Pr(>|t|)
(Intercept)         -0.15150    0.23249  -0.652    0.561
myGroup11            0.03927    0.23249   0.169    0.877
myGroup21           -0.37273    0.23249  -1.603    0.207
myGroup11:myGroup21       NA         NA      NA       NA

Residual standard error: 0.465 on 3 degrees of freedom
Multiple R-squared: 0.5605,     Adjusted R-squared: 0.2675 
F-statistic: 1.913 on 2 and 3 DF,  p-value: 0.2914 

> coef(summary(myModel))
               Estimate Std. Error    t value  Pr(>|t|)
(Intercept) -0.15149826  0.2324894 -0.6516352 0.5611052
myGroup11    0.03926774  0.2324894  0.1689012 0.8766203
myGroup21   -0.37273117  0.2324894 -1.6032180 0.2072173

> summary(myModel)$coefficients
               Estimate Std. Error    t value  Pr(>|t|)
(Intercept) -0.15149826  0.2324894 -0.6516352 0.5611052
myGroup11    0.03926774  0.2324894  0.1689012 0.8766203
myGroup21   -0.37273117  0.2324894 -1.6032180 0.2072173
4

3 回答 3

3

为什么不从拟合模型中提取系数：

> coef(myModel)
             (Intercept)                myGroup1b 
             -0.48496169              -0.07853547 
          myGroup2second myGroup1b:myGroup2second 
              0.74546233                       NA

这似乎是最简单的选择。

na.action与此无关。请注意，您没有通过na.action = na.pass您的示例。

na.actionNA是用于处理传递给模型拟合的数据的全局选项，通常与公式结合使用；它也是一个函数的名称na.action()。dataR 从参数和公式中表示的模型的符号表示构建所谓的模型框架。此时，NA将检测到任何变量，并且默认选项na.action是用于通过na.omit()删除任何变量NA的样本来从数据中删除。NA有一些替代方案，最有用的是na.exclude()，它们会在拟合期间移除，但会在拟合值、残差等的正确位置NA添加回来。阅读和更多信息，以及更多信息。NA?na.omit?na.action?options

于 2012-06-07T15:24:59.050 回答
1

summary.lm 的文档说“在返回对象中省略了别名系数，但通过print方法恢复了”。似乎没有参数可以控制此省略。除了使用@Gavin Simpson 建议的 coef(summary(myModel)) 之外，还有另一种解决方法。您可以创建一个矩阵

nr <- num_regressors - nrow(summary(myModel)$coefficients) ##num_regressors shall be defined previously
nc <- 4
rnames <- names(which(summary(myModel)$aliased))
cnames <- colnames(summary(myModel)$coefficients)
mat_na <- matrix(data = NA,nrow = nr,ncol = nc,
           dimnames = list(rnames,cnames))

然后 rbind 两个矩阵：

mat_coef <- rbind(summary(myModel)$coefficients,mat_na)
于 2017-01-12T10:21:08.730 回答
-1

您也可以将汇总拟合表转换为数据框（其中 NA 的变量丢失）：

fit <- as.data.frame(summary(fit)$coefficients)

然后按名称提取系数：

fit["age", "Pr(>|z|)"]

如果已删除“年龄”，则在尝试从数据框中提取年龄的 P 值时，您将获得 NA

于 2018-07-04T11:53:28.837 回答