Поиск наиболее вероятного значения с выхода пробоотборника MCMC

Question

Я относительно новичок в PyMC, и у меня есть быстрый вопрос о выходе сэмплера MCMC. Я хотел бы найти наиболее вероятное значение (максимум заднего) моих переменных, найденное сэмплером MCMC. Есть ли быстрый способ сделать это? Предположительно, значения переменных в максимуме задней части, найденные сэмплером MCMC, могут существенно отличаться от значений, обнаруженных методами MAP PyMC.

Благодарим разработчиков за предоставление PyMC. Это очень полезно для моей работы. Этот вопрос также был отправлен в группу Google PyMC.

Answer 1

Редактировать: Как указывает Ахмед Фасих, я неправильно понял вопрос. Его совет в комментарии к вопросу - это способ сделать это =)

---

Эта функциональность встроена в PyMC. Класс MAP в основном пространстве имен pymc принимает массив переменных pymc (или класс Model) и предоставляет метод fit.

map = mc.MAP(model) #or [var1, var2, .. ] 
map.fit() 

тогда все переменные pymc будут установлены как максимальный апостериор.