Błąd komunikatu WinBUGS: Oczekiwana struktura słowa kluczowego

Question

Jestem nowym użytkownikiem WinBUGS i nie mogłem uzyskać kodu poniżej. Model jest syntactically correct (jest to hierarchiczny model logitowy z efektami losowymi), ale po kliknięciu na load data pojawia się komunikat o błędzie expected key word structure. Co to znaczy? Czy ktoś może mi pomóc w pracy z poniższym kodem? Mój zestaw danych jest większy, ale aby uprościć pytanie, pracuję tutaj tylko z n = 2 (liczba grup) i k = 5 (liczba osób w każdej grupie).

model{ 
for(i in 1:n){ 
    for(j in 1:k){ 
    yij[i,j] ~ dbern(p[i,j]) 
    logit(p[i,j]) <- alpha + beta*xij[i,j] + ui[i] 
    } 
    ui[i] ~ dnorm(0,tau) 
} 
alpha ~ dnorm(0,0.001) 
beta ~ dnorm(0,0.001) 
tau ~ dunif(0,100) 
} 

z danymi:

list(n=2, k=5, yij=structure(.Data=c(1, 1, 1, 1, 1, 0, 0, 0, 0, 0), .Dim=c(2,5)), 
xij = c(0.0494063719, -0.078101264, 0.2748560749, 0.1245743393, -2.531242809, .6849338859, 0.5302062384, 0.7828560148, -0., 0.5128471157), 
ui = c(0.5031197054, 0.5031197054, 0.5031197054, 0.5031197054, 0.5031197054, -2.13785637, -2.13785637, 2.13785637, -2.13785637, -2.13785637)) 

Answer 1

xij zapowiada się matryca w modelu błędów, ale trzeba się jako wektor w danych. Powinno to zadziałać:

list(n=2, k=5, yij=structure(.Data=c(1, 1, 1, 1, 1, 0, 0, 0, 0, 0), .Dim=c(2,5)), 
xij = structure(.Data=c(0.0494063719, -0.078101264, 0.2748560749, 0.1245743393, -2.531242809, .6849338859, 0.5302062384, 0.7828560148, -0., 0.5128471157), .Dim=c(2,5)), 
ui = c(0.5031197054, 0.5031197054, 0.5031197054, 0.5031197054, 0.5031197054, -2.13785637, -2.13785637, 2.13785637, -2.13785637, -2.13785637))