Монте-Карло - равномерное распределение - PullRequest
0 голосов
/ 10 мая 2018

У меня есть эта проблема:

Каждую неделю три заправочные станции возле вашего дома продаются по цене за галлон X1, X2 и X3, соответственно, где Xi - IID равномерных случайных величин на [3.80, 4.20]. Вы покупаете газ на станция с самой низкой ценой. Используйте простые методы Монте-Карло для оценки средней цены за галлон газа, который вы заплатите, с ошибкой оценки не более $ 0,05.

Может кто-нибудь выложить R-код для такого рода проблемы или подсказать, как к этому подойти?

1 Ответ

0 голосов
/ 10 мая 2018

Sure

Давайте начнем с выборки

set.seed(12345)

n <- 1000
stationA <- runif(n, min=3.80, max=4.20)
stationB <- runif(n, min=3.80, max=4.20)
stationC <- runif(n, min=3.80, max=4.20)

whatIpay <- pmin(stationA, stationB, stationC)
avgIpay  <- mean(whatIpay)
stdIpay  <- sd(whatIpay)

теперь вы должны выяснить, какое значение n нужно получить в пределах желаемой ошибки

...