PyMC中的Beta分佈式線性模型

Question

在我看來應該是一個問題而不是舊的Google Groups， 我一直在尋找參數化PyMC中的Beta分佈式線性模型，Chris Fonnesbeck建議是：

「只是一個測試版的新參數，其中\阿爾法= \畝\披和\測試=（1 \畝）\披 所以，你需要是一樣的東西： new_beta =拉姆達（ 'new_beta'，lambda x = x，mu = mu，phi = phi：beta_like（x，mu * phi，（1-mu）* phi））「

這很棒 - 我的下一個問題是在哪裏線性模型組件，我認爲是HOULD是對\畝，例如：

$$ \畝= \ EXP（B_0 + b_1x）$$

和用γ-分佈式\披：

phi = Gamma('phi', alpha=0.001, beta=0.001) 

這是正確的？ NB我們在這裏PyMC2

感謝 亞倫

Answer 1

這取決於你想怎麼你的線性模型影響測試模型。正如你在這裏所描述的mu，它不是beta的意思，它只是平均值的標準化常數。如果你想阿爾法和描述均值和公測的方差測試版，它是類似如下：

alpha = mu * (mu*(1-mu)/var - 1) 
beta = (1 - mu) * (mu*(1-mu)/var - 1) 

也許一個簡單的方法是在一個簡陋的mu和樣本大小nu方面：

alpha = mu * nu 
beta = (1-mu) * nu