Hello,

I am new to Stata and I am trying to run a Bayesian analysis. I have a model using Jags in R and tried to run this same model in Stata. The results are different. Would you be able to check whether the 2 models below are the same, or if there is anything I missed? I am using Stata 16 and I received the message that dbern was not supported by Stata.

Jags model in R:
model {
for (i in 1:N) {
fail[i] ~ dbern(pi[i])
probit(pi[i]) <-b[1] + b[2]*temp[i]
}
for (j in 1:2) {b[j] ~ dnorm(0, 0.01)}
}

Model in Stata:
bayesmh fail temp, likelihood(probit) prior({fail}, binom(({fail:temp, xb}),1)) prior({fail:temp}, normal(0,100)) prior({fail:_cons}, normal(0,100)) initial({fail:} 0) saving(challm3)

Many thanks,
Giovanna