負二項式，泊松伽馬混合 winbugs

[英]Negative binomial , Poisson-gamma mixture winbugs

Winbugs 陷阱錯誤 這只是數據樣本，model 問題來自 Ezra Hauer 8.3.2，道路安全回歸的藝術，model 提供了一個**錯誤未定義的真實結果。 ** model 的目標是完全貝葉斯和一步 model 而不是使用經驗貝葉斯。 結果應類似於 MLE，其中 beta0 為 ...

WinBUGS 錯誤 - 我的“未定義變量”是什么？

[英]WinBUGS Error - What is my "undefined variable"?

這是我第一次運行 WinBUGS。 我的 model 相當簡單，只有一個截距和兩個隨機效應。 但是，當我運行 model 時，我得到了一個完整的錯誤列表（見下文）。 第一個錯誤是“未定義的變量”。 我無法確定我的哪些變量在語法中未定義。 任何幫助是極大的贊賞。 這是我的 R 代碼。 我正在使用 R ...

需要逗號或右括號的 openBugs 錯誤。 我的錯誤在哪里？

[英]openBugs error with expected a comma or right parenthesis. where is my mistake?

列表（案例=結構（。數據=c（1.10000E+01，1.20000E+01，1.10000E+01，1.20000E+01，1.10000E+01，2.20000E+01，2.30000E+01，2.20000E+01 , 2.30000E+01, 2.20000E+01, 3.30000E+0 ...

OpenBUGS：分層 Model 無法編譯

[英]OpenBUGS: Hierarchical Model unable to compile

我根據以下數據在 OpenBUG 中實現了 model： 數據包括 10 項獨立研究——2 種不同葯物的 5 次試驗（Med = 0 或 1），然后是該試驗的感染參與者總數（感染）——以及試驗的參與者總數。 我原來的 model 可以寫成： 我可以像這樣在 OpenBUGS 中實現： 但是，然后我 ...

“變量名中不允許有空槽”（OpenBUGS，R2OpenBUGS）

[英]“empty slot not allowed in variable name” (OpenBUGS, R2OpenBUGS)

我是 OpenBUGS 的初學者，我通過 R2OpenBUGS R package 使用它。 我嘗試設置 state 空間 model 以識別非常嘈雜的數據中的對數正態信號。 經過多次試驗和錯誤，我設法獲得了此代碼，但我仍然收到以下錯誤消息：“變量名錯誤 pos 664 中不允許空槽”，我不明白。 ...

有沒有辦法在 JAGS 中獲取和存儲矩陣元素的位置？

[英]Is there a way to obtain and store positions of matrix element in JAGS?

我正在 R 中開發貝葉斯分層 model ，並在 JAGS 中使用 BUGS 代碼。 在我的 model 中，我有兩個矩陣，它們在同一個精確矩陣 position 中包含彼此相關的信息。 我的信息按行排列。 我按行對第一個矩陣 Distmat 應用數學運算： 我有興趣將 diffmat 每一行中每 ...

JAGS遞歸model

[英]JAGS recursive model

我是使用 JAGS/WinBUGS/STAN 語言進行貝葉斯統計分析的新手。 我正在嘗試安裝擴展的 SIR model，似乎錯誤位於第 2 行，因為我總是收到警告：“解析 model 文件時出錯：第 2 行靠近“for”的語法錯誤但是我無法解決這個錯誤，數據：y是矩陣形式（T x K），N向量大小 ...

winbugs：數組索引大於 model 數據類型分類（序數）的數組上限

[英]winbugs: array index is greater than array upper bound for model with data type kategoric (ordinal)

我需要幫助來查找我的 WINBUGS 代碼 (v. 1.4.3) 中的錯誤。 在“模型規范”步驟中，model 看起來語法正確。 但是，在我嘗試加載數據時，我收到了這個錯誤： 數組索引大於a 的數組上限有人可以幫我嗎？ 我的 model 提供如下： ...

普通 CAR u[22] 的鄰接矩陣必須是對稱的 - WINBUGS

[英]Adjacency matrix of normal CAR u[22] must be symmetric - WINBUGS

我正在嘗試為我的空間數據構建一個卷積模型，但它一直給我一個錯誤，說“普通 CAR u[22] 的鄰接矩陣必須是對稱的”。 任何人都可以給我一個想法，我做錯了什么？ 我認為這個問題與 obs 42 沒有任何鄰居有關。 謝謝你的時間！ 非常感謝。 ...

JAGS/BUGS 中的 if-else 條件

[英]If-else condition in JAGS/BUGS

我有很大的參數矩陣。 關鍵是在許多參數中，一些帶有任意索引的參數會導致值錯誤，我想修復它們。 玩具示例如下： 我看到了一些例子，比如15414303和46730232 ，但我無法解決如何在我的案例中應用這些技巧（或類似技巧）的問題。 有沒有一種簡單的方法可以在 JAGS/BUGS 中實現這 ...

將 WINBUGS 模型轉換為 PyMC3

[英]Convert WINBUGS model to PyMC3

我目前正在學習貝葉斯統計課程，我們可以使用任何包來計算求解模型，但所有示例都在 WINBUGS 中提供。 我更喜歡使用 Python 和 PyMC3。 我沒有太多使用 PyMC3 的經驗，可以使用一些幫助來了解如何將這個簡單的 WINBUGS 模型轉換為 PyMC3 模型。 示例 WINBUGS ...

隨機截距和斜率 model 在 Winbugs 中的個體變異性內具有相關性和復雜性

[英]Random intercept and slope model with correlation and complex within individual variability in Winbugs

我正在嘗試在 WinBUGS 中實現具有復雜誤差可變性的隨機截距和斜率。 我得到“預期的多變量節點”。 這是帶有錯誤消息的示例代碼。 任何幫助將非常感激。 請讓我知道是否需要進一步澄清。 ...

R，HPD（最高后驗密度）區間基於后驗樣本，WinBUGS

[英]R, HPD (Highest posterior density) interval based on samples from posterior, WinBUGS

如何從后驗樣本中計算 HPD（最高后驗密度）區間？ 我有四個參數，我從后驗參數分布中生成 1000 個樣本。 現在如何在 R 軟件中計算 HPD。 我使用了 package 代碼但我收到一個錯誤 其中“winbugsresult”是一個包含后驗樣本的列表。 我還使用了一個向量，但出現以下錯誤 我只 ...

為先驗定義約束時 JAGS 中的錯誤消息嘗試重新定義節點

[英]Error messgae in JAGS when defining constraints for priors Attempt to redefine node

我試圖從先前的分布中采樣並將某些總和和差異限制為大於 0，但是，當我運行此通常在 OpenBugs/WinBugs 中有效的代碼時，我收到以下錯誤“嘗試重新定義節點”。 我想知道是否有任何可能的方法來定義我的約束以防止發生此錯誤。 任何幫助將不勝感激。 謝謝 ...

數組索引大於 d 的數組上限

[英]array index is greater than array upper bound for d

有人可以幫我解決這個問題嗎？ 我使用了具有不同數量的研究 (NS) 和治療 (NT) 的其他數據集，並且效果很好。 任何幫助將不勝感激。 數據集如下： 在哪里 ...

WinBUGS中計算后驗分布的代碼語法

[英]Code syntax in calculating posterior distribution in WinBUGS

最近，我讀了《 The BUGS Book-貝葉斯分析實用入門》以學習WinBUGS。 WinBUGS描述后驗分布的方式使我感到困惑。 讓我們以本書中的示例4.1.1來說明： 假設我們觀察到死亡y的數量在給定的醫院進行高風險操作。 令n表示所執行的此類操作的總數，並假設我們希 ...

如何修復第4行中發生的“預期右括號”：winbugs中的tau0 *（t [j + 1] -t [j]）

[英]How to fix the 'expected right parenthesis' happening in the 4th line : tau0*(t[j+1]-t[j]) in winbugs

我需要使用winbugs估計四個參數： tau[i], mu1[i],mu2[i]and tau0[i] 。 當我點擊檢查模型時，它顯示第4行中tau0*和(t[j+1]-t[j])之間的expected right parenthesis 。 當我刪除(t[j+1]-t[j]) ，它顯 ...

基本層次貝葉斯分析代碼

[英]Code for Basic Hierarchical Bayesian Analysis

我用winBUGS編寫了以下代碼： 我在如何在R / JAGS中執行此操作時遇到了麻煩。 實際上，我什至不能完全確定此代碼試圖做什么（我想計算后驗嗎？）。 我以前從未使用過winBUGS，我是R的新手。這也是我的第一個貝葉斯課程，一旦引入了代碼，我就很迷路。 另外，如何計算比例 ...

“多元正態Y [1,1]的均值必須與Y [1,1]的分量數相同”-WinBUGS中的錯誤

[英]“mean of multivariate normal Y[1,1] must have the same number of components as Y[1,1]” - Error in WinBUGS

如何處理這樣的編譯錯誤：“多元正態Y [1,1]的均值必須與Y [1,1]的分量數相同”？ 這是我在WinBugs中的代碼和數據： 數據： ...

如何解釋回歸變量的后驗概率？

[英]How to interpret posterior probability of an regression variable?

我在WINBUGS中使用以下模型來運行分層貝葉斯回歸，其中beta是我的協變量： 如果我通過添加以下代碼來修改此模型： 然后，我可以評估β協變量（正或負）關聯的后驗概率（ PP ）。 我的beta值是： 使用此代碼， beta0的PP為0.959033 ，而be ...