패키지 의 함수 lmer()가 lme4p- 값을 생성하지 않았다는 인상을 받았습니다 ( lmerp- 값 및 그 밖의 모든 것을 참조하십시오 ).

이 질문에 따라 대신 MCMC에서 생성 된 p 값을 사용 했습니다. lme4혼합 모델의 중요한 효과 및이 질문 : 의 패키지에서 출력에서 p- 값을 찾을 수 없습니다lmer()lm4R .

최근에 memisc 라는 패키지 getSummary.mer()를 사용하여 모델의 고정 효과를 csv 파일로 가져 왔습니다 . 마술처럼 마치 pMCMC p- 값과 매우 일치 하는 열이 나타나고 (을 사용하여 처리 시간이 걸리지 않습니다 pvals.fnc()).

나는 코드를 잠깐 살펴 보았고 getSummary.merp- 값을 생성하는 줄을 발견했습니다.

p <- (1 - pnorm(abs(smry@coefs[, 3]))) * 2

이것은 p 값 lmer이 달리기보다는 출력 에서 직접 생성 될 수 있다는 것을 의미합니까 pvals.fnc? 이것이 의심 할 여지없이 ‘p- 값 페티쉬주의’토론을 시작한다는 것을 알고 있지만 관심이 있습니다. memisc전에 언급 한 적이 없습니다 lmer.

좀 더 간결 해지려면 : getSummary.mer()?

답변