2013-02-13 1 views
2

에 대한 요약 결과가 다른 이유는 무엇입니까? by nlme 패키지에서 lmList 패키지를 사용하여 다른 그룹 by by 기능을 기반으로 다양한 회귀 분석을 실행하십시오. 나는 조금 값 "이 표시 왜 내 말하는lmList vs lmList [x]]

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) 

아래 누군가가 나에게 설명해 주시겠습니까이다 내가 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. 

에 대한 참조 설명 주어진 요약 결과에 대한 혼란 스러워요 "은 r2와 r3의 결과와 r2의 결과가 r3의 결과와 일치합니까?

답변

5

lmList의 버전과 NLME 인수 pool이 그것을 요약 방법이다 :

잔류 표준 오류의 풀링 추정 표준 의 계산에 사용할지 여부를 나타내는 선택적인 논리 값 편차 또는 표준 오류를 나타냅니다.

난 당신이 호출 할 때 FALSE 해당 설정하면 summary 당신이 동일한 값을 얻을 것 같은데요.

+0

그 트릭을했다. 고맙습니다! – JPC

+0

내 허락을 받아들이는 것을 잊어 버려! – JPC

관련 문제