我正在尝试使用OpenBugs
在R中运行以下模型model {
# Likelihood.
for ( i in 1 : N ) {
Y[i] ~ dnorm( mu[i], tau )
mu[i] <- alpha+beta*x[i]}
# Prior.
sigma <- 1/sqrt(tau)
alpha ~ dnorm(0, 1.0E-6)
beta ~ dnorm(0, 1.0E-6)
tau ~ dgamma(1.0E-3,1.0E-3)
}
datalist = list(
x = c( 1,2,3,4,5 ) ,
Y = c( 1,3,3,3,5) ,
N = 5
)
它遇到以下错误: “模型在语法上是正确的 数据加载 模型编译 无法为UpdaterGamma.Updater类型的节点生成初始值 samplesSet中的错误(c(“alpha”,“beta”,“tau”)): 必须在监视器使用“
之前初始化模型但是,在为tau分配一个值后,即tau <-1.0,它可以正常运行。
任何人都知道什么是错的?