Ошибка в узлах [1] Узел не соответствует родителям - PullRequest
0 голосов
/ 09 ноября 2019

Я подгоняю измененные значения расширенной модели Вейбулла, используя 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?

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...