変数を列挙する（PyMCをPyMC3に移植する）

Question

PyMCからPyMC3にthis codeを翻訳してPyMC3から始めます。

私はこのセグメントを翻訳するかどうかはわかりません。

v = pymc.Beta('v', alpha=1, beta=alpha, size=N_dp) 
@pymc.deterministic 
def p(v=v): 
    """ Calculate Dirichlet probabilities """ 

    # Probabilities from betas 
    # this line creates the error: 
    value = [u*np.prod(1-v[:i]) for i,u in enumerate(v)] 
    # Enforce sum to unity constraint 
    value[-1] = 1-sum(value[:-1]) 

    return value 

z = pymc.Categorical('z', p, size=len(set(counties))) 

私はp(v)と最後の行にpを交換し、@pymc.deterministicを削除する必要がありますが、問題は、私はを列挙することができないということのようだと仮定v：ValueError: length not known: ViewOp [id A] 'v'。

翻訳のやり方や、ドキュメントの関連するビットへのリンク方法を教えていただけますか？ありがとう。

Answer 1

全体のコードブロックがで置き換えることができるように、ディリクレ分布は、実際にbuilt into pymc3次のとおりです。興味のために

with pm.Model(): 
    ... 
    v = pm.Beta('v', alpha=1, beta=alpha, shape=N_dp) 
    p = pm.Dirichlet('p', a=v, shape=N_dp) 
    ... 
    trace = pm.sample(20000)