#figure9.2, 在 100 個估計下, 樣本比率一致性的模擬  
#使用法: 將此檔案 (命名為 figure9.2.txt)  置於 R 的工作目錄中, 執行指令 source('figure9.2.txt')

y=c()
cnt=matrix(rep(0, 200000), nrow=100)

p=0.5
epsilon=0.03
lbd=p-epsilon
ubd=p+epsilon

for(n in 1:2000) y[n]=rbinom(1, n, p)/n

for(j in 1:2000){
	if(lbd<y[j] & y[j]<ubd) cnt[1,j]=1
	else cnt[1,j]=0
}


plot(1:2000, y, type='l',  xlim=c(0,2000), ylim=c(0, 1), xlab='投擲數', 
	ylab='100 個 p 的估計')

for(i in 2:100){
	for(n in 1:2000) y[n]=rbinom(1, n, p)/n
	lines(1:2000, y, type='l',  xlim=c(0,2000), ylim=c(0, 1))

      for(j in 1:2000){
		if(lbd<y[j] & y[j]<ubd) cnt[i,j]=1
		else cnt[i,j]=0
	}
}

percnt=epsilon*100

for(i in seq(100, 2000, by=100)) cat('投擲數 n =', i, '時, 在正負',  
		percnt, '個百分點內的比率 =', sum(cnt[, i])/100, '\n')

abline(h=p)
abline(h=lbd, lty=2, col=2)
abline(h=ubd, lty=2, col=2)