我在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[]未定义。我不确定如何解决这个问题,所以任何帮助都将不胜感激。
发布于 2013-04-18 02:11:32
您确定将向量变量CD41正确地传递给您的dataList吗?
如果使用R,请在将其传递给OpenBUGS之前尝试调用CD41<-as.vector(CD41)。记住,对于BUG和R,1列矩阵和向量是不同的东西。
https://stackoverflow.com/questions/16094450
复制相似问题