The R package contains different MCMC schemes to estimate the regression coefficients of a binomial (or binary) logit model within a Bayesian framework: a data-augmented independence MH-sampler, an auxiliary mixture sampler and a hybrid auxiliary mixture (HAM) sampler. All sampling procedures are based on algorithms using data augmentation, where the regression coefficients are estimated by rewriting the logit model as a latent variable model called difference random utility model (dRUM).

Package: | binomlogit |

Type: | Package |

Version: | 1.1 |

Date: | 2014-03-12 |

License: | GPL-3 |

The main functions are `dRUMIndMH`

for independence Metropolis-Hastings sampling, `dRUMHAM`

for hybrid auxiliary mixture sampling and `dRUMAuxMix`

for auxiliary mixture sampling in the dRUM representation of a binomial logit model. The function `IndivdRUMIndMH`

is designed to work with binary instead of binomial observations to estimate the regression coefficients of a logit model. All four functions simulate the posterior distribution of the regression coefficients of the logit model and return the MCMC draws. For more details about the models and the estimation procedures see **References**. The results are returned in a list of class "binomlogit"" and can be displayed by using the functions `print`

and `summary`

. The `plot`

method returns a plot of the MCMC draws and the acf.
The Caesarean birth data set (`caesarean`

) is provided in three different formats to serve as exemplary input for the different MCMC schemes.

Agnes Fussl

Maintainer: Agnes Fussl <avf@gmx.at>

Agnes Fussl, Sylvia Fruehwirth-Schnatter and Rudolf Fruehwirth (2013),
"Efficient MCMC for Binomial Logit Models". *ACM Transactions on Modeling and Computer Simulation 23*, 1, Article 3, 21 pages.

Sylvia Fruehwirth-Schnatter and Rudolf Fruehwirth (2010),
"Data augmentation and MCMC for binary and multinomial logit models."
In *Statistical Modelling and Regression Structures - Festschrift in Honour
of Ludwig Fahrmeir*, T. Kneib and G. Tutz, Eds. Physica-Verlag, Heidelberg, pp. 111-132.

`dRUMIndMH`

, `dRUMAuxMix`

, `dRUMHAM`

, `IndivdRUMIndMH`

To evaluate the efficiency of the different MCMC samplers as in the paper by Fussl, Fruehwirth-Schnatter and Fruehwirth (2013), use e.g. the initial sequence estimator `initseq`

by Geyer (package mcmc), the `numEff`

function by Rossi (package bayesm) or the appropriate functions of the coda package by Plummer, Best, Cowles and Vines.

1 2 | ```
# please run the examples in:
# dRUMIndMH, dRUMAuxMix, dRUMHAM and IndivdRUMIndMH
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.