伯努利散文的模拟

Question

考虑一个公平的硬币。如果你观察尾巴，让 Hn=0，如果你观察脸，让 Hn=1。那么P(Hn=0)=P(Hn=1)=1/2。

Sn=H1+...+Hn 的分布是参数为 n 和 p=1/2 的二项式分布。

我们使用此分布来找到 Yn=Sn/n 的概率，即人脸的比例。

在 r 中模拟这个。

这是我在R中实现该算法的方法（根据文字不确定是否正确，请检查）：

n=1000
h=rbinom(n,1,1/2)
y=cumsum(h)/n
plot(n,y)

当我运行程序时，标记一个错误

绘图（n，y）

xy.coords(x, y, xlabel, ylabel, log) 中的错误： “x”和“y”长度不同

我不明白为什么自从n=1000 之后我输入y 就知道有多少元素返回了1000 值。

那为什么我会收到错误消息？

Answer 1

n = 1000，作为数字。所以：

1:n => c(1, 2, 3, ..., 998, 999, 1000)

plot(1:n,y)