Я подгоняю измененные значения расширенной модели Вейбулла, используя JAGS в R.
Вот некоторые данные и код для соответствия измененной расширенной модели Вейбулла в JAGS:
# Data
set.seed(2)
m=50;lambda=0.5;beta=1;alpha=10;
u=runif(m);A = (lambda*alpha)**-1;B = log2(1-u)
t<-alpha*(log2(1- A*B))**(1/beta)
#load jags, specify a jags model.
library(runjags)
j.model="
data{
C<-1000
for(i in 1:n){
ones[i]<- 1
}
}
model{
for(i in 1:n){
# likelihood
MEW[i] <- lambda*beta*(t[i]/alpha)^(beta-1)*exp((t[i]/alpha)^beta + lambda*alpha*(1-exp((t[i]/alpha)^beta)))
ones[i] ~ dbern(MEW[i]/C)
}
#priors
alpha ~ dgamma(1,1)
beta ~ dnorm(0, .0001) I(0, )
lambda ~ dnorm(0, .0001) I(0, )
}
"
#load data as list
data <- list(t = t, n = length(t))
# run jags model
jags.out <- run.jags(j.model,
data=data,
n.chains=1,
monitor=c('alpha','beta','lambda')
)
summary(jags.out)
Этомодель подходит отлично. Однако, когда я запускаю модель JAGS, происходит сбой с ошибкой:
Ошибка в узлах [1] Узел несовместим с родителями
Как я могу реализовать в JAGS?