R프로그래밍 기말고사 정리 14- 모비율의 신뢰구간

개념

비교

이항 분포만 가능

 

모비율에 대한 신뢰구간의 신뢰도 비교 시물레이션

입력값 : 모비율(p), 표본의 크기(n), 자료 set의 크기(m), 신뢰도(alpha)

 

입력

p=0.5;n=5;m=100000;alpha=0.95

코드

n=c(5,10,25,50)
p=0.5;m=100000;alpha=0.95
a=alpha+(1-alpha)/2
qn=qnorm(a)
count=0
for(i in 1: length(n)){
for(j in 1: m){
x=rbinom(n[i],1,p)
phat=sum(x)/n[i]
se=sqrt(phat*(1-phat))/sqrt(n[i])
L=phat-qn*se
U=phat+qn*se
if(p>=L && p<=U)count=count+1
}
P=count/m*100
cat("표본의 수 ",n[i],"포함확률 = ",P,"\n")
count=0
}

결과

 

불규칙적? 왜?