# 擲硬幣500次
N <- 500
flipsequence <-sample(x=c(0,1),size = N,replace = TRUE,prob = c(0.5,0.5))
r <- cumsum(flipsequence) ##統計每項前出現1的次數矩陣
n <- 1:N ##1~500的順序
runprop = r/n

#繪圖
plot(n,runprop,type="o",log="x",
xlim=c(1,N), ylim=c(0.0,1.0), cex.axis=1.5, ##cex.axis設計字體粗細
main="Running probability of Heads")

#添加一條直線y=0.5
lines(c(1,N),c(0.5,0.5),lty=3)

#添加標註
flipletters <- paste(c("T","H")[flipsequence[1:10]+1], collapse = "")
displaystring <- paste("FlipSequence= ", flipletters,‘...‘,sep="" )
text(5,0.9, displaystring, adj = c(0,1),cex=1.3)
#顯示模擬500次得到Head的頻率
text(100,0.9, paste("End probability = ",runprop[N]),adj=c(1,0), cex=1.3)

R擲硬幣500次，出現真面朝上的概率