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

如果我记得很清楚，如果输入数据或输出参数存在问题，就会出现此错误。我在您的代码中发现了几个潜在的问题：

我认为这个datos <- list("edad","enfer", "n")是错的：你应该使用像 datos <- list(edad = edad, enfer = enfer, n = n)，因为你想要传递实际数据，而不仅仅是R变量名！不要将此与输出参数混淆！
phi节点已计算但未受监控。我不确定这可能是一个问题，但WinBUGS有时会以这种方式表现得很奇怪。试一试！
初始值应随机生成，而不是像这样固定（如rnorm(1, 0, 0.01)等）

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

1 个答案: