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

Question

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

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

Например, в качестве входного распределения я использую равномерное распределение, где min, например, равно 2, а max равно 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