伯努利论文的模拟

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)