Если X = число головок, показывающих, когда брошены три монеты, найдите P (X = 1) и E (X) .
Эксперимент имеет следующее распределение:
Согласно моим расчетам,
P (X = 1) = 0,375 ≈ 0,40
E (X) = 1,50
А, 
, где n = # повторенийэксперимент.
Исходный код
noOfExperiments = 10000;
mySample <- sample(c(0,1,2,3), noOfExperiments, replace = T)
outcomeCount <- length(which(mySample==1))
prob <- outcomeCount / noOfExperiments
eX <- sum(mySample)/noOfExperiments
Согласно приведенному выше коду, я получил
P (X = 1) = 0,2518
E (X) = 1,4917 ≈ 1,50
Кажется, что значение E (X) идет хорошо, но значение P (X = 1) идет неправильно.
Что не так с моим кодом?
Редактировать:
Я написалследующий код на основе комментария Edward Carney .
noOfExperiments = 1000;
mySample <- replicate(noOfExperiments, sum(sample(c(0, 1), 3, replace=T)) )
outcomeCount <- length(which(mySample==1))
prob <- outcomeCount/noOfExperiments
eX <- sum(mySample)/noOfExperiments
Примечание: связанный вопрос .