Rmarkdown ошибка RJAGS
Каждый раз, когда я пытаюсь создать уценку в R, я получаю следующее сообщение об ошибке:
"Error in jags.model(model.file, data = data, inits = inits.values, n.chains = n.chains, : Nothing to compile calls: <Anonymous>... withVisible -> eval -> jags -> jags.model ->. Call In addition: Warning message: In sink() : no sink to remove Execution halted."
Код работает нормально, проблема как раз в том, когда я пытаюсь его связать. У меня есть кусок кода:
library(R2jags)
setwd("~/[...]")
getwd()
#Model:
sink("model1.txt")
cat("
model
{
for(i in 1:N){
y[i] ~ dnorm(mu[i],tau)
mu[i] <- alpha + beta * (x[i]-x.bar)
}
alpha ~ dnorm(0, 0.0001)
beta ~ dnorm(1,1)
tau ~ dgamma(.25,.25)
sigma <- 1/sqrt(tau)
}
",fill = TRUE)
sink()
#Data:
x = c(1,2,3,4,5)
y = c(1,3,3,3,5)
N = 5
x.bar = 3
jags.data = list("x","y","N","x.bar")
#Parameters:
jags.params = c("alpha", "beta", "tau", "sigma")
#Initial Values:
jags.inits = function(){
list("alpha" = 0, "beta" = 1, "tau" = 1)
}
#Fit Model:
lab1.sim = jags(jags.data, jags.inits, jags.params,
model.file = "model1.txt",
n.chains = 3, n.iter = 11000, n.burnin = 1000)
Я использую Windows 10, Rx64 3.2.3 и RStudio 0.99.903.
1 ответ
Решение
Когда вы пытаетесь связать вашу уценку, он запускает приемник для создания.txt с кодом модели и перезаписывает его. Markdown испортит.txt и сделает его пустым файлом каждый раз, когда вы вяжете его. Так что там нечего компилировать, и вы получите эту ошибку.
Решение: запустите sink(...), чтобы создать свой.txt на консоли и закомментировать его в источнике уценки.