jags

    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功能相同的呼叫運行正常。我

    1熱度

    1回答

    保存尖齒/輸出WinBUGS軟件我跑在R,工作正常尖齒/ WinBUGS軟件,但我不能輸出保存到文本或CSV文件(或者將工作) TEST.sim<- jags(data=jags.data, parameters.to.save=jags.params, n.iter=200000,jags.seed=123, model.file="~/Docum

    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

    7熱度

    1回答

    我想寫建模多晶話題車型WinBUGS軟件/尖齒模型(正是這個文件 - >http://www.ryanmcd.com/papers/mg_lda.pdf) 在這裏,我想選擇基於特定值的不同分配。 對於例如:我想這樣做 `if (X[i] > 0.5) { Z[i] ~ dcat(theta-gl[D[i], 1:K-gl]) W[i] ~ dcat(phi-gl[z[i], 1:V]) }

    4熱度

    1回答

    jags.model需要包含BUGS模型的文件名作爲其第一個參數。 爲了在一個腳本中包含所有內容,我有時使用writeLines命令將BUGS模型寫入文件。例如, library(rjags) writeLines(" model { for (i in 1:length(Y1)) { Y1[i] ~ dnorm(Beta0, Beta1) } Bet

    3熱度

    1回答

    我是JAGS的新手,並且正在嘗試理解digsval()如何在JAGS中爲刪失的數據工作。我對粗糙數據進行建模,其中每個數據點只有上限和下限(而不是真值)。下面是如何,我認爲它應該工作的一個簡單的例子: 一些上限和下限爲每點: > head(lim) L U [1,] 14.98266 15.68029 [2,] 21.21827 21.91590 [3,] 18.34953 1

    3熱度

    2回答

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