古いGoogleグループではなく、ここで質問すべきだったことを取り上げて、PyMCでベータ分散線形モデルをパラメータ化することを検討していました.Chris Fonnesbeckは、次のように提案しました:
「\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 を使用しています
ありがとうアーロン