2

我试图从 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 的文件路径更短，确保工作目录中没有更多的“数据”文件，使用数据变量名称的名称列表等。

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

4

1 回答 1

5

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

  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)等）

于 2013-01-31T13:47:38.667 回答