확률과 통계 · p값

p값

p값은 귀무가설 H0이 참이라고 가정하고 계산한, 관측 결과만큼 또는 더 극단적인 결과가 나올 확률입니다. 여기서 기계는 평균적으로 한 봉지에 500g을 채워야 합니다. 곡선은 실제로 그렇다면 n봉지의 평균이 어디에 놓일지 보여 줍니다.

이 시각화에서 보는 p값

가설검정은 기본 주장인 귀무가설 H0H_0 아래에서 데이터가 이례적인지 묻습니다. 여기서 H0H_0은 기계가 평균적으로 한 봉지에 μ=500\mu = 500g을 채운다는 주장입니다. 각 봉지의 무게는 서로 독립이며 σ=10\sigma = 10g인 정규분포를 따릅니다. 한 번의 연구에서 nn봉지를 재고 평균 xˉ\bar x를 구합니다. p값은 H0H_0이 참이라고 가정했을 때, 관측 평균만큼 또는 더 극단적인 평균이 나올 확률입니다. 극단적인 방향은 대립가설 H1H_1이 정합니다. 유의수준 α=0.05\alpha = 0.05인 검정은 p<0.05p < 0.05일 때 경보를 울립니다.

p는 정확히 무엇의 확률일까요? H0H_0 아래에서 평균 Xˉ\bar X는 평균 500, 표준오차 10/n10/\sqrt{n}인 정규분포를 따르므로 관측 평균의 위치는 z=(xˉ−500)/(10/n)z = (\bar x - 500)/(10/\sqrt{n})입니다. H1 ⁣:μ<500H_1\!: \mu < 500이면 p값은 zz의 왼쪽 면적 Φ(z)\Phi(z), μ>500\mu > 500이면 오른쪽 면적, μ≠500\mu \neq 500이면 양쪽 꼬리 면적 2Φ(−∣z∣)2\Phi(-|z|)입니다. 반복 연구는 같은 사실을 개수로 보여 줍니다. 실제 평균 충전량이 500g인 기계에서 평균이 이만큼 또는 더 극단적인 연구의 비율은 p 근처에 안정됩니다. 방향은 데이터를 보기 전에 정해야 합니다. 데이터를 본 뒤 두 단측 p값 중 작은 쪽을 고르면 H0H_0이 참일 때 경보가 울릴 확률이 10%로, 명목 비율의 두 배가 됩니다. 방향을 미리 고정하면 H0H_0 아래에서 p는 0부터 1까지 균등하게 분포하므로 오경보 확률은 정확히 α\alpha입니다.

왜 p는 H0H_0이 참일 확률이 아닐까요? p는 전적으로 H0H_0이 성립하는 세계에서 계산됩니다. 즉 P(이만큼 또는 더 극단적인 평균∣H0)P(\text{이만큼 또는 더 극단적인 평균} \mid H_0)이지 P(H0∣데이터)P(H_0 \mid \text{데이터})이 아닙니다. 조건을 뒤집으려면 p가 쓰지 않는 두 양이 필요합니다. 적게 채우는 기계가 원래 얼마나 흔한지, 그리고 그런 기계에서 검정이 경보를 울릴 확률인 검정력입니다. 25봉지에서 3g이 부족하면 검정력은 약 Φ(1.5−1.645)\Phi(1.5 - 1.645)입니다. 베이즈 정리의 오즈 형식에서, 경보는 부족 대 정상의 오즈를 검정력 ÷ α\alpha배로 만듭니다. 여기서는 약 8.85배이므로 사전 오즈 1 : 4가 약 2.2 : 1이 됩니다. 같은 경보라도 적게 채우는 기계가 80%이면 경보 중 정상 기계에서 나온 비율이 약 2.7%이고, 5%이면 약 68%입니다. 이 비율은 경보가 울렸다는 사실만 사용합니다. xˉ\bar x의 구체적인 값에는 더 많은 정보가 있습니다.

왜 아주 작은 차이도 유의할 수 있을까요? 관측 부족량이 고정되어 있으면 ∣z∣=∣xˉ−500∣n/σ|z| = |\bar x - 500|\sqrt{n}/\sigma는 n\sqrt{n}에 따라 커지므로, 더 많은 봉지를 잴수록 왼쪽 꼬리 p값이 작아집니다. 봉지 수가 네 배가 되면 ∣z∣|z|는 두 배가 됩니다. 왼쪽 단측 검정에서 500g보다 0.5g 적은 평균은 1083봉지부터 α=0.05\alpha = 0.05에서 유의합니다. 작은 p는 μ\mu가 정확히 500이라면 이런 데이터가 이례적이라는 뜻이지, 차이가 실질적으로 중요하다는 뜻은 아닙니다. 반대도 성립하지 않습니다. 실제로 부족해도 경보가 울리지 않는 경우가 흔하므로 p>0.05p > 0.05가 기계가 정상임을 보여 주지는 않습니다. 차이의 크기는 신뢰구간과 함께 보고하세요. α=0.05\alpha = 0.05인 양측 검정이 μ=500\mu = 500을 기각하는 것은 95% 구간 xˉ±z∗σ/n\bar x \pm z^*\sigma/\sqrt{n}(z∗≈1.96z^* \approx 1.96)에 500이 들어 있지 않은 것과 동치입니다.

살펴볼 점

  • 위쪽 봉지는 평균이 정확히 xˉ\bar x가 되도록 배치한 예시 데이터입니다. 검정에는 xˉ\bar x, nn, σ\sigma만 쓰입니다.
  • 아래쪽 축은 항상 500의 양쪽으로 표준오차 여섯 배씩을 표시하므로 곡선의 모양이 유지됩니다. 위쪽 축 아래의 틀이 이 범위를 표시하며, nn이 커지면 좁아집니다.

조작 방법

  • 아래쪽에서 드래그해 xˉ\bar x를 옮길 수 있습니다. xˉ\bar x나 방향 레버를 움직이면 같은 반복 결과를 다시 셉니다. nn을 바꾸면 두 모드 모두 처음부터 시작합니다.

관련 항목