R JAGS:节点中的父值无效(组合dcat和dnorm)

时间:2017-07-28 19:19:17

标签: r jags rjags

在使用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()中的所有参数需要使用此数据“验证”?对这种系统有更深入的了解会很棒。

0 个答案:

没有答案