lmer와 p- 값에 대한 혼란 : memisc 패키지의 p- 값과 MCMC의 p- 값은 어떻게 비교됩니까?

패키지 의 함수 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()?

— 텀블 다운
소스

getSummary.mer

p

$p$

p

$p$ memisc

@JasonMorgan 이것은 질문에 대한 합리적인 대답으로 보입니다.

— Glen_b-복지 주 모니카

@JasonMorgan 나는 당신이 말한 것에 동의하지만 mcmcsamp()많은 문제로 인해 현재 사용할 수 없습니다 ( Status of mcmcsamp자세한 내용은 glmm.wikidot.com/faq 에서 섹션을 확인할 수 있습니다 ). 현재로서는 아마도 (파라 메트릭?) 부트 스트랩이 실행 가능하고 구현하기 어렵지 않다고 생각합니다. 기능 bootMer()은 서비스 될 수 있습니다.

— usεr11852

@JasonMorgan 내 의견으로는 귀중한 의견이 실제로 질문에 대답하기 때문에 복사 붙여 넣기가 좋을 것이라고 생각합니다. (단, p- 값이 p- 값이 아닌 빠른 검사로만 사용되거나 근사치가 좋지 않을 가능성이 높은 경우에 대한 설명을 조금만 확대해도됩니다. 또는 MCMC가 안전한 내기 인 이유는 무엇입니까?)

— Glen_b -Reinstate Monica

OP에 대한 다른 기술 / 용어 설명. 방법을 사용하여 얻은 memiscp- 값은 관측 된 테스트 통계를 Wald 통계로 처리하여 얻은 p- 값입니다 ( 이 경우 t 를 Wald z 로 처리 ). 이러한 테스트는 "큰 샘플"가정에 의존하므로 샘플 크기가 커짐에 따라 점점 더 신뢰할 수 있습니다. 내가 아는 한 MCMC 기반 가치는 그러한 가정에 의존하지 않습니다. 어쨌든 Wald 테스트와 그에 대한 대안에 대해 조금만 읽으면 질문에 대한 추가 정보를 얻을 수 있습니다.

— Jake Westfall

getSummary.mer $p$ $p$ memisc

— 제이슨 모건
소스