#figure9.1, 在單一估計下, 針對不同投擲數, 樣本比率的模擬  
#使用法: 將此檔案 (命名為 figure9.1.txt)  置於 R 的工作目錄中, 執行指令 source('figure9.1.txt')

y=c()

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

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

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

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