我正在編寫一篇文章中Dirichlet-Multinomial posterior的推導,我已閱讀過該文章,並且遇到了將分佈總結爲1的問題。下面是代碼的非簡化形式: def pcn(X, n, N, c, alpha):
pnc = np.math.factorial(np.sum([n[i] for i in range(len(n))]))/ \
np.product(
我正在使用rJAGS構建一個多層次的貝葉斯模型,我想指定一個Cauchy先前的幾個參數。有沒有辦法在JAGS中執行此操作,還是需要切換到STAN?我的JAGS模型如下。我想用Cauchy替換dnorm分佈,但JAGS找不到標準R Cauchy分佈,例如, dcauchy,pcauchy model_string <- "model{
for (i in 1:n){
y[i] ~ dbi