Używam następującego kodu w OpenBUGS do przeprowadzenia analizy:

model
{
for(i in 1:467)
{
probit(p[i])<-gamma0+gamma1*drug[i]+gamma2*CD41[i]
R[i]~dbern(p[i])
junk[i]<-ID[i]
}
gamma0~dnorm(0,.0001)
gamma1~dnorm(0,.0001)
gamma2~dnorm(0,.0001)
}
ID[] drug[] CD41[] R[]
1   0   114 NA
2   1   40  NA
3   1   12  0
4   0   15  0
....
END

I otrzymuję następujący błąd: Zmienna CD41 [] nie jest zdefiniowana. Nie wiem, jak rozwiązać ten problem, więc każda pomoc byłaby bardzo mile widziana.