1

我正在尝试使用介绍性文本中的示例从 R 运行 WinBUGS。实际上，我使用的是相同的示例，并且与这个人有一个非常相似的问题：

R2WinBUGS - 警告消息

但我收到一条不同的错误消息，并且此线程中的建议（例如，以管理员身份运行 R，定义 BUGS 目录）没有奏效。

我的电脑归我的大学所有，所以我怀疑权限问题是有问题的。我最初在 Program Files 文件夹中安装了 WinBUGS，但在阅读此站点后，我尝试将其安装在我的 Users 文件夹中。无论如何，我都收到了相同的错误消息。

这是我收到的警告消息：

警告消息：运行命令 '"C:/Users/crowe106/WinBUGS14/WinBUGS14.exe" /par "C:/Users/crowe106/Desktop/School Stuff/R/WinBUGS/script.txt"' 的状态为 28462

WinBUGS 运行，但 R 崩溃并且未创建对象 out。有谁知道状态 28462 是什么意思？

谢谢您的帮助。

克里斯

这是代码：

library(R2WinBUGS)          # Load the R2WinBUGS library

# Save BUGS description of the model to the working directory
sink("model.txt")

cat("
model {

#Priors
population.mean~dunif(0,5000)          # Normal parameterized by precision
precision <- 1/population.variance     # Precision = 1/variance
population.variance <- population.sd * population.sd
population.sd~dunif(0,100)

# Likelihood
for(i in 1:nobs) {
mass[i]~dnorm(population.mean, precision)
}
}
", fill=TRUE)
sink()

# Package all the stuff to be handed over to WinBUGS
# Bundle data
win.data <- list(mass=y1000, nobs=length(y1000))

# Function to generate starting values
inits <- function()
list(population.mean=rnorm(1,600), population.sd=runif(1,1,30))

# Parameters to be monitored (= to estimate)
params <- c("population.mean", "population.sd", "population.variance")

# MCMC settings
nc <- 3          # Number of chains
ni <- 1000       # Number of draws from posterior (for each chain)
nb <- 1          # Number of draws to discard as burn-in
nt <- 1          # Thinning rate

# Start Gibbs sampler: Run model in WinBUGS and save results in object called out
out <- bugs(data=win.data, inits=inits, parameters.to.save=params, model.file="model.txt",
        n.thin=nt, n.chains=nc, n.burnin=nb, n.iter=ni, debug=TRUE, DIC=TRUE, working.directory=getwd(),
        bugs.directory = "C:/Users/crowe106/WinBUGS14")

我最初在没有 bugs.directory 行的情况下运行它并得到了相同的错误消息。

4

0 回答 0