0

在使用 Jags 和 Rjags 时，我是个新手。

modelString = "
model {
for (i in 1:N){
    y[i] ~ dt(mu,tau,nu)

}

tau <-   nuless2/(sig^2) 
nuless2 <- nu / nuless
nuless <- (nu-2)
sig ~ dnorm(0.02045457,10000000)
mu ~ dnorm(0.0013942308,4000000000)
nu ~ dcat(pi[])
pi <- c(0,0,3,3,3,3,3,3,3)
}"
writeLines(modelString,con="model.txt")
line_data = list("y"=ret,"N"=length(ret))

init_value = list("mu"=0.003)
model <- jags.model("model.txt", data=line_data, n.chains=2)

我一直在努力让这段代码正常工作，但我不知道它给了我错误。

错误：节点 nuless2 中的错误
无效的父值

我四处看看问题所在并阅读了大部分讨论，但我看不出任何原因。

如果有人能告诉我我的愚蠢错误在哪里，那就太好了。

另外，我对 JAGS 很陌生，想知道系统的流程是如何工作的。例如，当我定义y ~ dt()时，我还为 提供数据y。通过这样做，我是否告诉系统这就是数据以及dt()需要使用这些数据“验证”的所有参数？对这样的系统有更深入的了解会很棒。

4

0 回答 0