Блог

Я искал в библиотеке scipy какие-либо встроенные модули для подгонки байесовской кривой, и я не могу их найти. Все, что я нашел, это: scipy.optimize.curve_fit

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

Байесовский вывод не является частью библиотеки SciPy — он просто выходит за рамки scipy . Существует ряд отдельных модулей python, которые занимаются этим, и, похоже, вы действительно пропустили довольно много из них — в первую очередь реализации алгоритмов Монте-Карло с цепью Маркова pymc и emcee которые, вероятно, являются наиболее используемыми пакетами MCMC. Они оба относительно просты в настройке, но, на мой взгляд emcee , с ними проще начать.

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