Подняв вопрос о том, что должно было быть здесь, а не в старых группах Google, я попытался параметризовать линейную модель с бета-распределением в PyMC, которую Крис Фоннесбек предложил как:
«просто репараметризация бета-версии, где \alpha = \mu \phi и \beta = (1-\mu) \phi. Итак, все, что вам нужно, это что-то вроде: new_beta = Lambda('new_beta', lambda x=x , mu=mu, phi=phi: beta_like(x, mu*phi, (1-mu)*phi))"
и это здорово — мой следующий вопрос — куда втыкать компонент линейной модели, который, я думаю, должен быть на \mu, например:
$$\mu = \exp(b_0+b_1x)$$
и с гамма-распределением \phi:
phi = Gamma('phi', alpha=0.001, beta=0.001)
Это верно? NB мы на PyMC2 здесь
Спасибо, Аарон.