r - 为什么我会得到不同的 lmList 与 lmList[[x]] 的汇总结果

Question

我使用lmListfrom nlmepackage 来运行基于不同组的各种回归 by 函数。鉴于我看到的解释，我对总结结果有点困惑summary.lmList

The `summary.lm` method is applied to each lm component of object to produce
summary information on the individual fits, which is organized into a list of
summary statistics. 

The returned object is suitable for printing with the print.summary.lmList
method.

我的意思是下面

set.seed(123)
t <- data.frame(
        name=sample(c("a","b","c"),size=500,replace=T),
        x=1:500,
        y=1:500+rnorm(100)
      )
ta <- t[t$name=="a",]
lma <- lm(y~x,ta)
lmL <- lmList(y~x | name,t)
r1 <- summary(lmL)  
r2 <- summary(lmL[["a"]])  
r3 <- summary(lma)

有人可以向我解释为什么 r1 中显示的“a”值与 r2 和 r3 中的值不匹配，而 r2 中的值与 r3 中的值匹配吗？

Answer 1

nlme中的版本lmList及其摘要方法有一个参数：pool

一个可选的逻辑值，指示是否应在计算标准差或汇总标准误差时使用残差标准误差的汇总估计。

我猜如果你FALSE在打电话时将其设置为，summary你会得到相同的值。