Сообщения об ошибках OpenBUGS: ожидается, что оператор коллекции c

Я не мог заставить код ниже работать. это иерархическая односторонняя модель ANOVA, но когда я нажимаю «Загрузить данные», появляется сообщение об ошибке, которое ожидается оператором сбора c. Что это обозначает? Может ли кто-нибудь помочь мне работать с кодом ниже, пожалуйста? Мой набор данных больше, но для упрощения вопроса я работаю здесь только с сезоном = 4 (количество групп) и n = 5 (количество субъектов в каждой группе).

model{
for(i in 1:n){
    Length[i] ~ dnorm(mu.l[i], tau[1]) 
    mu.l[i] <- alpha[j[i]]
}
for (p in 1:J){
    alpha[p]~dnorm(mu, tau[2])
}
mu~dnorm(0,0.0001)
for(k in 1:2){
    tau[k]<-pow(sigma[k],-2)
    sigma[k]~dunif( 0, 10)
}
}

#initials    
list( ort=1.0, alpha=c(NA, 0,0,0), tau=c(1 ,1))  

#Data
list(n=5, j=4)  
j[] length[]  
1   17.00
1   17.50
1   17.50
1   16.20
1   18.00
2   13.70
2   17.4
2   17.70
2   16.40
2   17.70
3   16.4
3   15.00
3   19.60
3   14.70
3   18.00
4   18.20
4   13.60
4   17.30
4   17.3
4   14.5

bayesian winbugs winbugs14
person Sedat    schedule 04.02.2015    source источник

Ответы (1)


arrow_upward
2
arrow_downward

Несколько проблем в вашем коде.

  1. В вашей модели у вас есть J и Length. В ваших данных у вас есть j (дважды) и length, т.е. нет Length или J.

  2. В ваших начальных значениях у вас есть ort, который не используется нигде в модели, и tau, который не является стохастическим (возможно, вы имеете в виду sigma?).

Ваша приведенная выше модель должна работать со следующим (нужна ли вам эта модель или нет, решать вам):

#Initials
list( alpha=c(NA, 0,0,0))  

#Data
list(n=5, J=4, j=c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4), Length=c(17,17.5,17.5,16.2,18,13.7,17.4,17.7,16.4,17.7,16.4,15,19.6,14.7,18,18.2,13.6,17.3,17.3,14.5))
person guyabel    schedule 04.02.2015
comment
Спасибо гжабель. Извините за опечатку кода. Проблемы, о которых вы упомянули, из-за потери перевода с турецкого :) Я сменил модель, как вы рекомендуете. И это хорошо. - person Sedat; 04.02.2015
comment
@gjabel, не могли бы вы увидеть эту проблему stackoverflow.com/questions/58810405/ - person Uddin; 12.11.2019