我在 OpenBUGS 中使用以下代码进行分析:
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
我收到以下错误:未定义变量 CD41[]。我不知道如何解决这个问题,所以任何帮助将不胜感激。