data(iraqVote)
## probit model
glm1 <- glm(y ~ gorevote + rep,
data=iraqVote,
family=binomial(link=probit))
require(MCMCpack)
albertChib <- MCMCprobit(y ~ gorevote + rep,
data=iraqVote,
burnin=0,
B0=1e-3,
mcmc=1e4,
verbose=0)Run the code above in your browser using DataLab