Efficient MCMC for Binomial Logit Models

Share:

Description

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).

Details

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.

Author(s)

Agnes Fussl

Maintainer: Agnes Fussl <avf@gmx.at>

References

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.

See Also

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.

Examples

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