Начало работы
Пример
Вот простой пример, демонстрирующий пакет в действии:
using Turing
using StatsPlots
# Определим простую нормальную модель с неизвестными средним и дисперсией.
@model function gdemo(x, y)
s² ~ InverseGamma(2, 3)
m ~ Normal(0, sqrt(s²))
x ~ Normal(m, sqrt(s²))
y ~ Normal(m, sqrt(s²))
end
# Запустим сэмплер, соберем результаты.
chn = sample(gdemo(1.5, 2), HMC(0.1, 5), 1000)
# Обобщим результаты.
describe(chn)
# Построим график и сохраним результаты.
p = plot(chn)
savefig("gdemo-plot.png")