R

我在哪裡聲明 Stan 中的先驗參數?

  • April 7, 2014

我正在使用rstan. 這些模型中的每一個在頂級分佈之前都有一個單元信息。單位先驗的參數最初存儲在 Rlist中。我的問題是:

  1. list我是否從存儲其餘數據的位置傳遞特定模型的單元信息先驗參數?
  2. 我在哪裡聲明這些先驗參數?在參數塊中?

謝謝。

我會區分先驗分佈和先驗分佈的參數。當我開始使用 Stan 時,我會將參數設置為先驗分佈,就像一些值一樣。所以在模型步驟中,我會有類似的東西

model {    
    mu ~ normal(0, 1)
    y ~ normal(mu, s)
}

對於 y 分佈的平均係數的正態先驗。

然而,隨著我更多地使用 Stan,我傾向於將它們包含在數據步驟中(即,您引用的列表中其餘數據的存儲位置)。以這種格式,我將補充數據步驟

data {
   real mu_prior_1
   real<lower=0> mu_prior_2
}

並將模型步驟調整為類似

model {    
    mu ~ normal(mu_prior_1, mu_prior_2)
    y ~ normal(mu, s)
}

我一直這樣做的主要原因是因為它使我更容易更改先驗而不更改 Stan 代碼的其餘部分。

引用自:https://stats.stackexchange.com/questions/92915

comments powered by Disqus