[R 프로그래밍] 유의수준과 검정력 함수 그리기

유의수준과 검정력 함수 그리기

 

 R을 이용해서 유의수준과 검정력 함수를 그려보겠습니다.

 정규분포에서 표본을 36개 추출하였고 분산은 4로 가정하겠습니다.

 또 귀무가설 $H_0$ <= 0, 대립가설 $H_1$ > 0 으로 설정하였습니다.

 검정기준(기각역)은 표준평균이 0.5보다 클 때 귀무가설을 reject 하겠습니다.

## 유의수준을 구하는 함수
prob2 <- function(x,n)
{
  z = (0.5-x)/(2/sqrt(n))
  return(1-pnorm(z))
}

n <- 36

#모평균이 0일 때 유의수준
prob2(0,n)

# 검정역과 유의수준 시각화
curve(prob2(x,n),-2,2,xlab=expression(mu),ylab="Prob(Reject H0)")
abline(h=0)
lines(c(0,0),c(0,prob2(0,n)),col="red")