패키지 의 함수 lmer()
가 lme4
p- 값을 생성하지 않았다는 인상을 받았습니다 ( lmer
p- 값 및 그 밖의 모든 것을 참조하십시오 ).
이 질문에 따라 대신 MCMC에서 생성 된 p 값을 사용 했습니다. lme4
혼합 모델의 중요한 효과 및이 질문 : 의 패키지에서 출력에서 p- 값을 찾을 수 없습니다lmer()
lm4
R
.
최근에 memisc 라는 패키지 getSummary.mer()
를 사용하여 모델의 고정 효과를 csv 파일로 가져 왔습니다 . 마술처럼 마치 p
MCMC p- 값과 매우 일치 하는 열이 나타나고 (을 사용하여 처리 시간이 걸리지 않습니다 pvals.fnc()
).
나는 코드를 잠깐 살펴 보았고 getSummary.mer
p- 값을 생성하는 줄을 발견했습니다.
p <- (1 - pnorm(abs(smry@coefs[, 3]))) * 2
이것은 p 값 lmer
이 달리기보다는 출력 에서 직접 생성 될 수 있다는 것을 의미합니까 pvals.fnc
? 이것이 의심 할 여지없이 ‘p- 값 페티쉬주의’토론을 시작한다는 것을 알고 있지만 관심이 있습니다. memisc
전에 언급 한 적이 없습니다 lmer
.
좀 더 간결 해지려면 : getSummary.mer()
?
답변
getSummary.mer
memisc