Analysis of marginally specified semi-nonparametric models for clustered binary data

Generalized linear mixed models are widely used for analyzing clustered data. If the primary interest is in regression parameters, one can proceed alternatively, through the marginal mean model approach. In the present study, a joint model consisting of a marginal mean model and a cluster-specific conditional mean model is considered. This model is useful when both time-independent and time-dependent covariates are available. Furthermore our model is semi-parametric, as we assume a flexible, smooth semi-nonparametric density of the cluster-specific effects. This semi-nonparametric density-based approach outperforms the approach based on normality assumption with respect to some important features of 'between-cluster variation'. We employ a full likelihood-based approach and apply the Monte Carlo EM algorithm to analyze the model. A simulation study is carried out to demonstrate the consistency of the approach. Finally, we apply this to a study of long-term illness data.     

What are the advantages of MCMC based inference in latent variable models?

Recent developments in Markov chain Monte Carlo [MCMC] methods have increased the popularity of Bayesian inference in many fields of research in economics, such as marketing research and financial econometrics. Gibbs sampling in combination with data augmentation allows inference in statistical/econometric models with many unobserved variables. The likelihood functions of these models may contain many integrals, which often makes a standard classical analysis difficult or even unfeasible. The advantage of the Bayesian approach using MCMC is that one only has to consider the likelihood function conditional on the unobserved variables. In many cases this implies that Bayesian parameter estimation is faster than classical maximum likelihood estimation. In this paper we illustrate the computational advantages of Bayesian estimation using MCMC in several popular latent variable models.     

Markov switching quantile autoregression

This paper considers the location‐scale quantile autoregression in which the location and scale parameters are subject to regime shifts. The regime changes in lower and upper tails are determined by the outcome of a latent, discrete‐state Markov process. The new method provides direct inference and estimate for different parts of a non‐stationary time series distribution. Bayesian inference for switching regimes within a quantile, via a three‐parameter asymmetric Laplace distribution, is adapted and designed for parameter estimation. Using the Bayesian output, the marginal likelihood is readily available for testing the presence and the number of regimes. The simulation study shows that the predictability of regimes and conditional quantiles by using asymmetric Laplace distribution as the likelihood is fairly comparable with the true model distributions. However, ignoring that autoregressive coefficients might be quantile dependent leads to substantial bias in both regime inference and quantile prediction. The potential of this new approach is illustrated in the empirical applications to the US inflation and real exchange rates for asymmetric dynamics and the S&P 500 index returns of different frequencies for financial market risk assessment.     

Bayesian estimation of log odds ratios from R × C and 2 × 2 × K contingency tables

In this paper, Bayesian estimation of log odds ratios over R × C and 2 × 2 × K contingency tables is considered, which is practically reasonable in the presence of prior information. Likelihood functions for log odds ratios are derived for each table structure. A prior specification strategy is proposed. Posterior inferences are drawn using Gibbs sampling and Metropolis–Hastings algorithm. Two numerical examples are given to illustrate the matters argued.     

Women investors, 'that nasty south sea affair' and the rage to speculate in early eighteenth-century England

The excursions of the five unmarried Hastings sisters and their widowed friend Jane Bonnell into the stock market show how changes in the availability of credit and the services offered by banks in the early eighteenth century had an impact on ordinary citizens. At the time of the South Sea Bubble all six bought South Sea shares through their bank. But their trading activities and investment strategies differed and had different outcomes, showing there are no easy associations between gender and ideas of risk or safe investment.     

Modeling and calculating the effect of treatment at baseline from panel outcomes

We propose and examine a panel data model for isolating the effect of a treatment, taken once at baseline, from outcomes observed over subsequent time periods. In the model, the treatment intake and outcomes are assumed to be correlated, due to unobserved or unmeasured confounders. Intake is partly determined by a set of instrumental variables and the confounding on unobservables is modeled in a flexible way, varying both by time and treatment state. Covariate effects are assumed to be subject-specific and potentially correlated with other covariates. Estimation and inference is by Bayesian methods that are implemented by tuned Markov chain Monte Carlo methods. Because our analysis is based on the framework developed by Chib [2004. Analysis of treatment response data without the joint distribution of counterfactuals. Journal of Econometrics, in press], the modeling and estimation does not involve either the unknowable joint distribution of the potential outcomes or the missing counterfactuals. The problem of model choice through marginal likelihoods and Bayes factors is also considered. The methods are illustrated in simulation experiments and in an application dealing with the effect of participation in high school athletics on future labor market earnings.     

A threshold stochastic volatility model with explanatory variables

In this paper, we introduce a threshold stochastic volatility model with explanatory variables. The Bayesian method is considered in estimating the parameters of the proposed model via the Markov chain Monte Carlo (MCMC) algorithm. Gibbs sampling and Metropolis–Hastings sampling methods are used for drawing the posterior samples of the parameters and the latent variables. In the simulation study, the accuracy of the MCMC algorithm, the sensitivity of the algorithm for model assumptions, and the robustness of the posterior distribution under different priors are considered. Simulation results indicate that our MCMC algorithm converges fast and that the posterior distribution is robust under different priors and model assumptions. A real data example was analyzed to explain the asymmetric behavior of stock markets.     

A Bayesian encompassing test using combined value-at-risk estimates

The Value at Risk (VaR) is a risk measure that is widely used by financial institutions in allocating risk. VaR forecast estimation involves the conditional evaluation of quantiles based on the currently available information. Recent advances in VaR evaluation incorporate conditional variance into the quantile estimation, yielding the Conditional Autoregressive VaR (CAViaR) models. However, the large number of alternative CAViaR models raises the issue of identifying the optimal quantile predictor. To resolve this uncertainty, we propose a Bayesian encompassing test that evaluates various CAViaR models predictions against a combined CAViaR model based on the encompassing principle. This test provides a basis for forecasting combined conditional VaR estimates when there are evidences against the encompassing principle. We illustrate this test using simulated and financial daily return data series. The results demonstrate that there are evidences for using combined conditional VaR estimates when forecasting quantile risk.     

Estimation of Technical Inefficiencies with Heterogeneous Technologies

This paper considers the measurement of firm's specific (in)efficiency while allows for the possible heterogeneous technologies adopted by different firms. A flexible stochastic frontier model with random coefficients is proposed to distinguish technical inefficiency from technological differences across firms. Posterior inference of the model is made possible via the simulation-based approach, namely, Markov chain Monte Carlo method. The model is applied to a real data set which has also been considered in Christensen and Greene (1976), Greene (1990), Tsionas (2002), among others. Empirical results show that the regression coefficients can vary across firms, indicating the adoption of heterogeneous technologies by different firms. More importantly, we find that, without considering this possible heterogeneity, the inefficiency of firms can be over-estimated.     

Tailored randomized block MCMC methods with application to DSGE models

In this paper we develop new Markov chain Monte Carlo schemes for the estimation of Bayesian models. One key feature of our method, which we call the tailored randomized block Metropolis–Hastings (TaRB-MH) method, is the random clustering of the parameters at every iteration into an arbitrary number of blocks. Then each block is sequentially updated through an M–H step. Another feature is that the proposal density for each block is tailored to the location and curvature of the target density based on the output of simulated annealing, following  and  and Chib and Ergashev (in press). We also provide an extended version of our method for sampling multi-modal distributions in which at a pre-specified mode jumping iteration, a single-block proposal is generated from one of the modal regions using a mixture proposal density, and this proposal is then accepted according to an M–H probability of move. At the non-mode jumping iterations, the draws are obtained by applying the TaRB-MH algorithm. We also discuss how the approaches of Chib (1995) and Chib and Jeliazkov (2001) can be adapted to these sampling schemes for estimating the model marginal likelihood. The methods are illustrated in several problems. In the DSGE model of Smets and Wouters (2007), for example, which involves a 36-dimensional posterior distribution, we show that the autocorrelations of the sampled draws from the TaRB-MH algorithm decay to zero within 30–40 lags for most parameters. In contrast, the sampled draws from the random-walk M–H method, the algorithm that has been used to date in the context of DSGE models, exhibit significant autocorrelations even at lags 2500 and beyond. Additionally, the RW-MH does not explore the same high density regions of the posterior distribution as the TaRB-MH algorithm. Another example concerns the model of An and Schorfheide (2007) where the posterior distribution is multi-modal. While the RW-MH algorithm is unable to jump from the low modal region to the high modal region, and vice-versa, we show that the extended TaRB-MH method explores the posterior distribution globally in an efficient manner.