Я пытаюсь запустить WinBUGS из R, используя пример из вводного текста. На самом деле, я использую тот же пример и имею очень похожую проблему с этим парнем:
R2WinBUGS — предупреждающие сообщения
но я получаю другое сообщение об ошибке, и предложения в этой теме (например, запустить R от имени администратора, определить каталог BUGS) не сработали.
Мой компьютер принадлежит моему университету, поэтому я подозреваю, что виновата какая-то проблема с разрешениями. Первоначально у меня был установлен WinBUGS в папке Program Files, но после прочтения этого сайта я попытался установить его в папку 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 и получил такое же сообщение об ошибке.
y1000
, чтобы другие могли воспроизвести ваш пример? - person guyabel   schedule 07.01.2015