r - 在 R 中绘制蒙特卡洛模拟

Question

我正在尝试创建一个图，该图绘制了掷硬币的蒙特卡洛模拟，其中 x 上的迭代和 y 上的概率。但是，我很难创建必要的功能。

到目前为止我有

set.seed(2738952)
coins<-function(n){
    sample(c(0,1),n,replace=TRUE)
}
mean(coins(10000)==1)

ggplot(mapping = aes(x=iteration,y=coins))+geom_line()

错误：美学必须是有效的数据列。有问题的美学：y = 硬币。您是否输入错误数据列的名称或忘记添加 stat()？

我不确定如何创建迭代以及概率。

Answer 1

这是你要找的东西吗？

iterations <- 100
flip_results <- coins(iterations)
# As @teunbrand suggested, use cumsum to avoid looping
pct_head <- pct_head <- cumsum(flip_results) / seq_len(iterations)
ggplot(mapping = aes(x=1:iterations,y=pct_head))+geom_line()