【问题标题】:This seems to work fine but keep having "command #Bugs:gen.inits cannot be executed (is greyed out)"这似乎工作正常,但一直有“命令 #Bugs:gen.inits 无法执行(灰显)”
【发布时间】:2015-05-24 09:34:25
【问题描述】:

我正在通过 R2WinBUGS 使用 WinBUGS 做一个简单的回归模型。 结果看起来没问题,但是 WinBUGS 一直显示这条线,我不知道为什么:

command #Bugs:gen.inits cannot be executed (is greyed out)

有人有想法吗? 谢谢!

调用WinBUGS的R代码:

# Prepare data
ni = 12 #Nb of tagging year
prL <- c(0.13,0.23,0.28,0.22,0.29,0.25,0.31,0.32,0.32,0.40,0.33,0.35) #Proportion of tag lost
input.tlr <- list(ni=ni,prL=prL)

# Generates intial values for two chains
nbchains = 2
init1 <- list(beta1=runif(1,-1,1),beta2=runif(1,-1,1),tauG=rgamma(1,1,1))
init2 <- list(beta1=runif(1,-1,1),beta2=runif(1,-1,1),tauG=rgamma(1,1,1))
inits <- list(init1,init2)

# Call WinBUGS
parameters <- c("gamma","beta1","beta2","tauG")
out.tlr <- bugs(data=input.tlr, inits=inits, parameters=parameters, "tlr.txt", n.chains=2, n.iter=20000, n.burnin=5000,n.thin=1,debug=TRUE)

WinBUGS 代码是:

model {
###### Tag lost rate model ######
for (i in 1:ni) {
  prL[i] ~ dnorm(gamma[i],tauG)
  gamma[i] <- beta1 * log(i-0.5) + beta2
  }
### Priors
    beta1 ~ dnorm(0,0.001)
    beta2 ~ dnorm(0,0.001)
    tauG ~ dgamma(0.001,0.001)
} # End model

WinBUGS 的输出是:

display(log)
check(C:/Users/AppData/Local/Temp/RtmpG8dxh1/tlr.txt)
model is syntactically correct
data(C:/Users/alebris/AppData/Local/Temp/RtmpG8dxh1/data.txt)
data loaded
compile(2)
model compiled
inits(1,C:/Users/AppData/Local/Temp/RtmpG8dxh1/inits1.txt)
chain initialized but other chain(s) contain uninitialized variables
inits(2,C:/Users/AppData/Local/Temp/RtmpG8dxh1/inits2.txt)
model is initialized
gen.inits()
command #Bugs:gen.inits cannot be executed (is greyed out)
thin.updater(1)
update(5000)
set(gamma)
set(beta1)
set(beta2)
set(tauG)
set(deviance)
dic.set()
update(15000)
coda(*,C:/Users/AppData/Local/Temp/RtmpG8dxh1/coda)
stats(*)

【问题讨论】:

    标签: r winbugs r2winbugs


    【解决方案1】:

    我认为这只是因为您自己提供了inits。 winBUGS 命令 gen.inits 用于自动生成 inits,但如果您手动提供 inits,它会“变灰”。所以我觉得没什么好担心的。

    【讨论】:

      猜你喜欢
      • 2011-05-05
      • 2023-03-29
      • 2017-09-04
      • 2012-05-06
      • 2013-02-17
      • 2017-10-09
      • 2017-08-03
      • 2014-05-02
      • 2011-09-11
      相关资源
      最近更新 更多