R2WinBUGS 错误 - 陷阱 - 不兼容的副本

Question

我正在尝试从 R 调用 winBUGS 来估计逻辑回归。我正在使用以下代码:

# Directorio de trabajo
setwd("~/3 Diplomado/7 Bayesiana/8t1")

# paquete para hablarse con WinBUGS desde R
library(R2WinBUGS)

# cargamos datos
reg <- read.table("enf.csv", header = TRUE, sep = ",")
edad <- reg$edad
enfer <- reg$efer
n <- length(reg$edad)

# Primeras filas de los datos
head(reg)

# Nombres de los datos para alimentar al modelo en WinBUGS
datos <- list("edad","enfer", "n")

# Construimos el modelo

modelo <- function(){
  for (i in 1:n) {
    enfer[i] ~ dbin(theta[i], 1)
    logit(theta[i]) < - beta0 + beta1 * edad[i] 
    phi[i]<-theta[i]/(1-theta[i])
      }
  beta0 ~ dnorm(0, 0.0001)
  beta1 ~ dnorm(0, 0.0001)

}

## escribimos el modelo a un archivo .bug file:
write.model(modelo, "modelo.bug")
## Para echarle un vistazo al modelo:
file.show("modelo.bug")

# Distribuciones iniciales
iniciales <- function(){
 list(beta0=1,beta1=1)
}

# Corremos esto en WinBUGS

   regresionlogistica.sim = bugs(datos,iniciales,
                 model.file="modelo.bug",
                 parameters=c("beta0","beta1"),
                 n.chains=1,n.iter=20000,n.burnin=5000,n.thin=1,
                 bugs.directory="C:/Users/pythonxy/Documents/3 Diplomado/7 Bayesiana/winbugs14",
                 working.directory="~/3 Diplomado/7 Bayesiana/8t1",          
                 codaPkg=FALSE)

WinBUGS 打开但出现错误。 “不兼容的副本”。我一直在寻找这个错误，发现很多人都遇到过这个问题。我已经尝试了建议的解决方案:model.bug 的较短文件路径，确保工作目录中不再有“数据”文件，使用数据变量名称的名称列表等。

但我仍然无法让它工作。我做错了什么？

Answer 1

如果我没记错的话，如果输入数据或输出参数有问题，就会出现这个错误。我在您的代码中发现了几个潜在的问题:

1. 我认为这 datos <- list("edad","enfer", "n")错了:你应该用 like datos <- list(edad = edad, enfer = enfer, n = n) ，因为你想传递实际数据，而不仅仅是 R 变量名!不要将它与输出参数混淆!

2. phi节点被计算但不被监控。我不确定这是否会成为问题，但 WinBUGS 有时会以这种方式表现得很奇怪。试试吧!

3. 初始值应该是随机生成的，不是固定的(比如rnorm(1, 0, 0.01)等)