Как сделать распространение ошибок, используя emcee

Question

Я использую конферансье , чтобы проанализировать некоторые данные, относящиеся к SN 1a.Прежде чем заняться реальными данными, я начал с симулированных данных.Полученные результаты кажутся хорошими:

Теперь я хочу вычислить исходные значения некоторых величин, таких как

f (q, j) = j - q ^ 2

с использованием полученных мной постеров для q и j .Кто-нибудь может мне помочь?

Answer 1

Если у вас есть необработанные задние выборки j и q , например, два пустых массива j и q, то при поэлементном вычислении этих массивов будут получены соответствующие выборкидля желаемой переменной.В вашем примере f = j - q**2.

Я знаю, что опции для захвата выборок преобразованных переменных, представляющих интерес непосредственно в семплере, существуют в других семплерах (Stan, PyMC3).Может быть, кто-то, кто знает ведущего (не меня), знает, как это сделать.