Модель иногда компилируется, иногда нет, хотя ничего не меняется - PullRequest
Новая
       24

Модель иногда компилируется, иногда нет, хотя ничего не меняется

0 голосов
/ 23 февраля 2019

Итак, у меня есть модель rjags, приведенная ниже.Моя проблема в том, что когда я хочу запустить jags.model(), он иногда компилируется, иногда нет.Таким образом, я мог буквально запустить функцию jags.model(), и при немедленном повторном запуске одного и того же jags.model() с теми же параметрами и без изменений ничего не происходит, я внезапно получаю сообщение об ошибке: 

Error in update.jags(object, n.iter, ...) : 
Error in node ((beta[2]/x1[5])^0.9)
Invalid parent values

У меня действительно нетИдея, как это может произойти, что-то не так с моим рабочим пространством?Ниже приведены мои данные и моя модель. 

    bayesnlm = "
model {

for (i in 1:length(y)) {

y[i] ~ dnorm(mu[i], tau)

mu[i] <- beta[3] + beta[4] * (x2[i]) +
beta[1] *(0.9/beta[2])*(beta[2]/x1[i])^1.9*exp((-(beta[2]/x1[i])^0.9))


residual[i] <- y[i] - mu[i]


}


# Define Prior distribution 

beta[1:4] ~ dmnorm.vcov(mu.p[1:4], Sigma[1:4,1:4])


tau ~ dgamma(1.0E-3, 1.0E-3)
} 
"

Мои входные данные выглядят так: 

$`x1`
[1]  14  75 258  26  70

$x2
[1]  14  75 258 815 859

$y
[1]  0.2091645  1.7007476 -0.2880298 -0.2880298  0.7063589

$Sigma
              [,1]          [,2]          [,3]          [,4]
[1,]  2.479069e+04  3.717394e+03 -1.940021e+01  5.853116e-03
[2,]  3.717394e+03  1.245803e+03  3.743637e+00 -1.778376e-03
[3,] -1.940021e+01  3.743637e+00  3.006899e-01 -1.496941e-04
[4,]  5.853116e-03 -1.778376e-03 -1.496941e-04  9.027231e-08

$mu.p
[1]  5.445813e+02  1.343918e+02 -5.689060e-01 -1.070608e-04

Здесь вы можете скопировать и вставить их: 

    list(x1 = c(14, 75, 258, 26, 70), x2 = c(14, 75, 258, 815, 859
), y = c(0.20916452356296, 1.70074760142997, -0.288029835726044, 
-0.288029835726044, 0.706358882851964), Sigma = structure(c(24790.6941908131, 
3717.39378904953, -19.400211288343, 0.00585311596063579, 3717.39378904953, 
1245.80270600402, 3.74363729374668, -0.00177837577140645, -19.400211288343, 
3.74363729374668, 0.300689871287482, -0.000149694071017865, 0.00585311596063579, 
-0.00177837577140645, -0.000149694071017865, 9.02723134020887e-08
), .Dim = c(4L, 4L)), mu.p = c(544.581343592283, 134.391793200685, 
-0.568905977830338, -0.000107060840574308))
...