2013-02-13 87 views
2

我使用的lmList包從nlme包中運行基於不同組的功能的各種迴歸。我有點困惑關於彙總結果給我看到summary.lmList爲什麼我會得到lmList vs lmList的不同總結結果[[x]]

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與r2和r3不匹配,而r2與r3匹配?

回答

5

lmList的版本和它的摘要方法在NLME具有一個參數pool

指示是否 殘餘標準誤差的彙總估計應在標準 計算中使用的可選邏輯值摘要的偏差或標準誤差。

我猜如果你設置爲FALSE當你撥打summary你會得到相同的值。

+0

做到了這一點。謝謝! – JPC 2013-02-13 03:33:48

+0

忘記接受這個,我的壞! – JPC 2014-07-15 22:54:47

相關問題