我试图理解收缩率，并且我建立了一个非常简单的模型，我认为应该在参数化二项式分布p的两个参数之间显示收缩率：

with pm.Model() as model:
    alpha = pm.HalfNormal('alpha',3)
    beta = pm.HalfNormal('beta',3)
    ps = pm.Beta('ps',alpha=alpha,beta=beta,shape=2)
    pm.Binomial('obs',p=ps,n=[6,125],observed=[1,110])

    traces = pm.sample(3000,cores=2,tune=500)

我假设由于ps [0]和ps [1]都共享较高的级别分布，因此从ps [1]抽取的次数明显多于ps [0]的概率，但是基于后验估计，看起来没有任何作用。