r2jags

    3熱度

    2回答

    我有一個關於在JAGS和BUGS中運行模型的細節的快速問題。 說我用n.burnin=5000,n.iter=5000和thin=2運行模型。這是否意味着程序將: 運行5000次迭代,並放棄結果;然後 運行另一個10,000次迭代,只保留每秒的結果? 如果我將這些模擬保存爲CODA對象,那麼所有的對象都保存了10,000次,或者只有稀疏的5000次?我只是想了解使用哪組迭代來製作ACF圖?

    2熱度

    1回答

    我被告知使用貝葉斯推理,而不是僅使用輪詢數據分析工作。但是,我有一個問題;我有一個小數據集,可以猜測各方的先前分佈情況,並且我有來自民意調查的數據。我如何從Gibbs模擬中獲得邊際值? prior <- a <- c(.30, .15, .15, .10, .10, .08, .12) polls <- data.frame(rbind( a <- c(.24, .23, .20, .

    5熱度

    1回答

    當使用jags.parallel,我得到以下錯誤: > out <- jags.parallel(win.data, inits, params, "Poisson.OD.t.test.txt", + nc, ni, nb, nt); Error in get(name, envir = envir) : invalid first argument 使用jags功能相同的呼叫運行正常。我

    0熱度

    1回答

    我正在運行以下模型,其中我爲每個投票周進行了測量。每個派對產生超過100個阿爾法矢量。可能的問題是,我怎樣才能以一種我可以繪製折線圖的方式儲存它們? model{ ## measurement for(i in 1:NPOLLS){ p1[i] ~ dnorm(alpha1[WoY[i]] + pollster1[org[i]], prec1[i]) p2[i]

    4熱度

    1回答

    我正在使用R2jags和CODA爲我的MCMC鏈生成一些診斷統計信息,但我遇到了麻煩。我想運行MCMC如下: modelfit <- jags(data=jags.data, inits=jags.inits, model.params, n.iter = 100000, model.file=jags.model, model.params) 的錯誤是: Error in v