Моделирование эссе Бернулли - PullRequest
0 голосов
/ 18 мая 2019

Рассмотрим честную монету. Пусть Hn = 0, если вы наблюдаете за хвостом, и Hn = 1, если вы наблюдаете за лицом. Тогда P (Hn = 0) = P (Hn = 1) = 1 / 2.

Распределение Sn = H1 + ... + Hn является биномиальным с параметрами n и p = 1 / 2.

Мы используем это распределение, чтобы найти вероятности для Yn = Sn / n, который является пропорцией граней.

Имитация этого в т.

Вот мой подход к реализации этого алгоритма в R (не уверен, если он корректен в соответствии с текстом, пожалуйста, проверьте его):

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

Когда я запускаю программу, отмечается ошибка

участок (п, у)

Ошибка в xy.coords (x, y, xlabel, ylabel, log): Длина 'x' и 'y' различаются

Я не понимаю, почему с n=1000 и когда я набираю y, чтобы узнать, сколько элементов вернуло 1000 значений.

Почему я получаю ошибку?

1 Ответ

2 голосов
/ 18 мая 2019

n = 1000, как число. Итак:

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

plot(1:n,y)

enter image description here

...