openbugs

问题 I am trying to do Bayesian NHST between two groups. Each group consists of many variables, and in order to avoid multiple hypothesis corrections I opted for Bayesian method. However my code in OpenBUGS is giving the following error: "Error in handleRes(res) : Internal "trap" error in OpenBUGS, or non-existent module or procedure called." I am using BRugs R-package with modelCheck(...) - model is syntactically correct and with modelData(...) - data loaded are appearing but when tried to

问题 This question was migrated from Cross Validated because it can be answered on Stack Overflow. Migrated 6 years ago . I'm using the following code in OpenBUGS to perform an analysis: 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 And I'm receiving the following error: Variable CD41[] is not