R 에서 수치적으로 계산 가능한 최저 P-value ?

[권경석]

예를 들어..

> summary(aov(i.b.con~factor(week)+factor(season)+factor(week)*factor(season)))
                             Df Sum Sq Mean Sq F value Pr(>F)    
factor(week)                  2  430.8   215.4   262.8 <2e-16 ***
factor(season)                2  970.8   485.4   592.3 <2e-16 ***
factor(week):factor(season)   4  567.5   141.9   173.1 <2e-16 ***
Residuals                   807  661.4     0.8                   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

이런 결과를 얻었는데

2e-16이라는 수치가 엄청 작은 수치긴 한데

이러한 동일한 값이 반복적으로 나온다는 것 (다른 실험 case에서도 지금까지 발견한 최저값이 2e-16 이었습니다)은

2e-16 이 수치적으로 산출할 수 있는 최저의 p-value 라는 의미일까요?

[김종연]

네 그렇습니다. 0 이상의 수 중 (64비트?) 컴퓨터가 저장할 수 있는 제일 작은 수일겁니다..

[김종연]

http://stats.stackexchange.com/questions/78839/how-should-tiny-p-values-be-reported-and-why-does-r-put-a-minimum-on-2-22e-1
참고가 될 링크를 찾아 걸어드립니다:-)

[조희준]

2e-16 보다 작은 exact p-value는 -log2(p-value) 아니면 -log10(p-value) 형태로 구할 수 있어요. summary에서 F-stat(t-stat) 값을 뽑아내신 다음 pf function(혹은 pt) 의 log.p option을 이용하시면 구하실 수 있어요.

[권경석]

답변 감사드립니다 !!