построение симуляции Монте-Карло в 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()

Ошибка: эстетика должна быть действительными столбцами данных. Проблемная эстетика (ы): у = монеты. Вы неправильно набрали имя столбца данных или забыли добавить 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()