我在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 []。我不知道如何解决这个问题,所以任何帮助将不胜感激。
答案 0 :(得分:2)
您确定要将向量变量CD41
正确传递给dataList吗?
如果正在使用R,请在将其传递给OpenBUGS之前尝试调用CD41<-as.vector(CD41)
。请记住,1个列矩阵和一个向量对于BUG和R都是不同的。