R - mc2d Монте-Карло пакет, уровень неопределенности

Question

У меня есть следующий вопрос относительно пакета mc2d для моделирования в Монте-Карло.

Учитывая узел mc, т. Е. Объект mc. Как мы можем получить неопределенность для значений распределения?

Например, в качестве распределения ввода я использую равномерное распределение, где min равно, например, равным 2, а максимальное значение равно 8. Учитывая это, мы создаем объект mc, применяем его к mc.

Резюме функция выдает значения, такие как медиана, среднее, 97,5% и т.д. и т.п.

Но, как я сказал, как можно получить оценку неопределенности для заданного значения?

Заранее благодарен!

Answer 1

Ну, вы должны собрать второй импульс Тогда

v = <x^2> - <x>^2 
u = sqrt(v)/sqrt(N-1) 

a = <x> +-u 

Чтобы сделать вещи более ясно, вы выборочные события

x = 2 + (8-2)*U(0,1) 

где-то в функции суммарного вы вычислить сумму событий

m = m + x 

поэтому после запуска N событий вы сообщаете mean=m/N

Вы должны добавить код, чтобы собрать второй импульс, что-то вроде

m2 = m2 + x*x 

Таким образом, после запуска можно вычислить

v = m2/N - mean*mean 
u = sqrt(v)/sqrt(N-1) 

и сообщить в виду с неопределенностью mean +-u