Байесовская инверсия простая модель - PullRequest
0 голосов
/ 17 января 2020

У меня есть следующий код, который пытается выполнить байесовский вывод над моделью просто поддерживаемого луча с учетом нескольких измерений.

https://github.com/hanzy1110/PyMC3_Examples/blob/master/modelo_barra.py

Уравнение, связывающее отклонение луча и нагрузку, выглядит следующим образом: V = (5/32) * (mu_p * L**4)/(mu_E * b * h**3)

Сначала я хочу проверить влияние неопределенности на модуль Юнга E, а затем учесть неопределенность в распределенной нагрузке p. «Набор данных», который я предоставляю, - это всего лишь 5 измерений, выполненных под одной и той же нагрузкой.

Когда я пытаюсь выполнить это, консоль зависает из-за отсутствия образцов, необходимых для создания объекта трассировки. Я хочу спросить следующее:

- Модель хорошо написана? -Как количество образцов указано в модели? (Консоль продолжает возвращать только 100 образцов в цепочке, даже после того, как я добавлю 1000 случайных чисел к наблюдаемому количеству) -Как я могу заставить его работать!?!

...