[R 프로그래밍] 모비율 추론하기 - 신뢰구간, 가설검정, 표본크기결정

 R을 이용해서 모비율을 추론해보겠습니다. (신뢰구간, 가설검정, 표본크기결정)

 이론은 여기에서 확인할 수 있습니다.

 

1. 데이터 입력

n <- 73380 # 표본공간
x <- 38100 # 성공횟수
p <- x/n   # 성공확률

 

2. 비율 신뢰구간

# 신뢰구간
prop.test(x,n)

 

 

 비율 95% 신뢰구간 = [0.5155923, 0.5228358]

 

3. 비율 가설검정

# 가설검정
prop.test(x,n,p=0.5,alternative="greater")

 

 

 p(확률)이 기각역에 포함되므로 귀무가설은 기각

 

4. 표본크기 결정

### 표본크기결정(비율)
SizePropCI <- function(Error, Conf.level=0.95)
{
  alpha <- (1-Conf.level)/2
  (qnorm(1-alpha)/Error)^2/4  
}

SizePropCI(0.031,0.95)

SizePropTest <- function(prop1,prop2,Alpha=0.05,Power=0.8)
{
  delta = prop2-prop1
  term <- qnorm(1-Alpha/2)*sqrt(prop1*(1-prop1))+qnorm(Power)*sqrt(prop2*(1-prop2))
  (term/delta)^2
}

SizePropTest(0.4,0.5)

 

 

구간추정을 위한 표본크기 = 1000

가설검정을 위한 표본크기 = 191