A comparison of multiple imputation with EM algorithm and MCMC method for quality of life missing data

This study investigated the performance of multiple imputations with Expectation-Maximization (EM) algorithm and Monte Carlo Markov chain (MCMC) method in missing data imputation. We compared the accuracy of imputation based on some real data and set up two extreme scenarios and conducted both empirical and simulation studies to examine the effects of missing data rates and number of items used for imputation. In the empirical study, the scenario represented item of highest missing rate from a domain with fewest items. In the simulation study, we selected a domain with most items and the item imputed has lowest missing rate. In the empirical study, the results showed there was no significant difference between EM algorithm and MCMC method for item imputation, and number of items used for imputation has little impact, either. Compared with the actual observed values, the middle responses of 3 and 4 were over-imputed, and the extreme responses of 1, 2 and 5 were under-represented. The similar patterns occurred for domain imputation, and no significant difference between EM algorithm and MCMC method and number of items used for imputation has little impact. In the simulation study, we chose environmental domain to examine the effect of the following variables: EM algorithm and MCMC method, missing data rates, and number of items used for imputation. Again, there was no significant difference between EM algorithm and MCMC method. The accuracy rates did not significantly reduce with increase in the proportions of missing data. Number of items used for imputation has some contribution to accuracy of imputation, but not as much as expected.     

Pair trading based on quantile forecasting of smooth transition GARCH models

Pair trading is a statistical arbitrage strategy used on similar assets with dissimilar valuations. We utilize smooth transition heteroskedastic models with a second-order logistic function to generate trading entry and exit signals and suggest two pair trading strategies: the first uses the upper and lower threshold values in the proposed model as trading entry and exit signals, while the second strategy instead takes one-step-ahead quantile forecasts obtained from the same model. We employ Bayesian Markov chain Monte Carlo sampling methods for updating the estimates and quantile forecasts. As an illustration, we conduct a simulation study and empirical analysis of the daily stock returns of 36 stocks from U.S. stock markets. We use the minimum square distance method to select ten stock pairs, choose additional five pairs consisting of two companies in the same industrial sector, and then finally consider pair trading profits for two out-of-sample periods in 2014 within a six-month time frame as well as for the entire year. The proposed strategies yield average annualized returns of at least 35.5% without a transaction cost and at least 18.4% with a transaction cost.     

Bayesian analysis of a Tobit quantile regression model

This paper develops a Bayesian framework for Tobit quantile regression. Our approach is organized around a likelihood function that is based on the asymmetric Laplace distribution, a choice that turns out to be natural in this context. We discuss families of prior distributions on the quantile regression vector that lead to proper posterior distributions with finite moments. We show how the posterior distribution can be sampled and summarized by Markov chain Monte Carlo methods. A method for comparing alternative quantile regression models is also developed and illustrated. The techniques are illustrated with both simulated and real data. In particular, in an empirical comparison, our approach out-performed two other common classical estimators.     

Closed-skew normality in stochastic frontiers with individual effects and long/short-run efficiency

This paper considers the estimation of Kumbhakar et al. (J Prod Anal. doi:10.1007/s11123-012-0303-1, 2012) (KLH) four random components stochastic frontier (SF) model using MLE techniques. We derive the log-likelihood function of the model using results from the closed-skew normal distribution. Our Monte Carlo analysis shows that MLE is more efficient and less biased than the multi-step KLH estimator. Moreover, we obtain closed-form expressions for the posterior expected values of the random effects, used to estimate short-run and long-run (in)efficiency as well as random-firm effects. The model is general enough to nest most of the currently used panel SF models; hence, its appropriateness can be tested. This is exemplified by analyzing empirical results from three different applications.     

Ordinal-response GARCH models for transaction data: A forecasting exercise

We use numerous high-frequency transaction data sets to evaluate the forecasting performances of several dynamic ordinal-response time series models with generalized autoregressive conditional heteroscedasticity (GARCH). The specifications account for three components: leverage effects, in-mean effects and moving average error terms. We estimate the model parameters by developing Markov chain Monte Carlo algorithms. Our empirical analysis shows that the proposed ordinal-response GARCH models achieve better point and density forecasts than standard benchmarks.     

Nonparametric Bayes inference for concave distribution functions

Bayesian inference for concave distribution functions is investigated. This is made by transforming a mixture of Dirichlet processes on the space of distribution functions to the space of concave distribution functions. We give a method for sampling from the posterior distribution using a Pólya urn scheme in combination with a Markov chain Monte Carlo algorithm. The methods are extended to estimation of concave distribution functions for incompletely observed data.     

A general framework for approximate sampling with an application to generating points on the boundary of bounded convex regions

We consider the problem of generating a sample of points according to some given probability distribution over some region. We give a general framework for constructing approximate sampling algorithms based on the theory of Markov chains. In particular, we show how it can be proven that a Markov chain has a limiting distribution. We apply these results to prove convergence for a class of so-called Shake-and-Bake algorithms, which can be used to approximate any absolutely continuous distribution over the boundary of a full-dimensional convex body.     

Bayesian semiparametric stochastic volatility modeling

This paper extends the existing fully parametric Bayesian literature on stochastic volatility to allow for more general return distributions. Instead of specifying a particular distribution for the return innovation, nonparametric Bayesian methods are used to flexibly model the skewness and kurtosis of the distribution while the dynamics of volatility continue to be modeled with a parametric structure. Our semiparametric Bayesian approach provides a full characterization of parametric and distributional uncertainty. A Markov chain Monte Carlo sampling approach to estimation is presented with theoretical and computational issues for simulation from the posterior predictive distributions. An empirical example compares the new model to standard parametric stochastic volatility models.     

Solidarity properties and a Doeblin decomposition for a class of non-Markovian stochastic processes

We are concerned with solidarity and a Doeblin decomposition for a class of non-Markovian discrete parameter stochastic processes. Since any such process is associated with a certain general Markov chain whose transition probability function has a special form, we use the theory of Markov chains with continuous components to this particular chain in order to get properties of the non-Markovian process. We illustrate our results on a model closely related to learning theory.     

Empirical Bayes analysis of log-linear models for a generalized finite stationary Markov chain

This article presents the empirical Bayes method for estimation of the transition probabilities of a generalized finite stationary Markov chain whose ith state is a multi-way contingency table. We use a log-linear model to describe the relationship between factors in each state. The prior knowledge about the main effects and interactions will be described by a conjugate prior. Following the Bayesian paradigm, the Bayes and empirical Bayes estimators relative to various loss functions are obtained. These procedures are illustrated by a real example. Finally, asymptotic normality of the empirical Bayes estimators are established.     

An analytical approximation formula for European option pricing under a new stochastic volatility model with regime-switching

In this paper, an analytical approximation formula for pricing European options is obtained under a newly proposed hybrid model with the volatility of volatility in the Heston model following a Markov chain, the adoption of which is motivated by the empirical evidence of the existence of regime-switching in real markets. We first derive the coupled PDE (partial differential equation) system that governs the European option price, which is solved with the perturbation method. It should be noted that the newly derived formula is fast and easy to implement with only normal distribution function involved, and numerical experiments confirm that our formula could provide quite accurate option prices, especially for relatively short-tenor ones. Finally, empirical studies are carried out to show the superiority of our model based on S&P 500 returns and options with the time to expiry less than one month.     

A Paradigm for Masking (Camouflaging) Information

This is an expository paper. Here we propose a decision-theoretic framework for addressing aspects of the confidentiality of information problems in publicly released data. Our basic premise is that the problem needs to be conceptualized by looking at the actions of three agents: a data collector, a legitimate data user, and an intruder. Here we aim to prescribe the actions of the first agent who desires to provide useful information to the second agent, but must protect against possible misuse by the third. The first agent is under the constraint that the released data has to be public to all; this in some societies may not be the case.
A novel aspect of our paper is that all utilities—fundamental to decision making—are in terms of Shannon's information entropy. Thus what gets released is a distribution whose entropy maximizes the expected utility of the first agent. This means that the distribution that gets released will be different from that which generates the collected data. The discrepancy between the two distributions can be assessed via the Kullback–Leibler cross-entropy function. Our proposed strategy therefore boils down to the notion that it is the information content of the data, not the actual data, that gets masked. Current practice of "statistical disclosure limitation" masks the observed data via transformations or cell suppression. These transformations are guided by balancing what are known as "disclosure risks" and "data utility". The entropy indexed utility functions we propose are isomorphic to the above two entities. Thus our approach provides a formal link to that which is currently practiced in statistical disclosure limitation.     

Replicating the Results in ‘A New Model of Trend Inflation’ Using Particle Markov Chain Monte Carlo

An article by Chan et al. ( 2013 ) published in the Journal of Business and Economic Statistics introduces a new model for trend inflation. They allow the trend inflation to evolve according to a bounded random walk. In order to draw the latent states from their respective conditional posteriors, they use accept–reject Metropolis–Hastings procedures. We reproduce their results using particle Markov chain Monte Carlo (PMCMC), which approaches drawing the latent states from a different technical point of view by relying on combining Markov chain Monte Carlo and sequential Monte Carlo methods. To conclude: we are able to reproduce the results of Chan et al. ( 2013 ). Copyright © 2015 John Wiley & Sons, Ltd.     

Issue Information

Much research studies US inflation history with a trend‐cycle model with unobserved components, where the trend may be viewed as the Fed's evolving inflation target or long‐horizon expected inflation. We provide a novel way to measure the slowly evolving trend and the cycle (or inflation gap), by combining inflation predictions from the Survey of Professional Forecasters (SPF) with realized inflation. The SPF forecasts may be treated either as rational expectations (RE) or updating according to a sticky information (SI) law of motion. We estimate RE and SI state‐space models with stochastic volatility on samples of consumer price index and gross national product/gross domestic product deflator inflation and the associated SPF inflation predictions using a particle Metropolis–Markov chain Monte Carlo sampler. The trend converges to 2% and its volatility declines over time—two tendencies largely complete by the late 1990s.     

Stock prices and monetary policy shocks: A general equilibrium approach

Empirical literature documents that unexpected changes in the nominal interest rates have a significant effect on real stock prices: a 100-basis point increase in the nominal interest rate is associated with an immediate decrease in broad real stock indices that may range from 2.2 to 9%, followed by a gradual decay as real stock prices revert towards their long-run expected value. We assess the ability of a general equilibrium New Keynesian asset-pricing model to account for these facts. We consider a production economy with elastic labor supply, staggered price and wage setting, as well as time-varying risk aversion through habit formation. We find that the model predicts a stock market response to policy shocks that matches empirical estimates, both qualitatively and quantitatively. Our findings are robust to a range of variations and parametrizations of the model.     

Inter‐state dynamics of invention activities, 1930–2000

We study the dynamics of the cross‐section distribution of patents per capita for the 48 continental US states from 1930 to 2000 using a discrete‐state Markov chain. We test for and find evidence in favor of the (knowledge) convergence hypothesis. The distribution of patents is converging to a limiting distribution that is significantly more concentrated than its initial distribution. States in the extreme are more mobile than states in the middle of the cross‐sectional distribution and are likely to move to the middle. However, the rate of convergence to the limiting distribution is ‘slow’. Copyright © 2006 John Wiley & Sons, Ltd.     

Patterns of rational default

In this paper, we move closer to the interests of empirical research by providing the entire distribution of default’s severity, not just a summary total. The total price of a mortgage’s default option already provides a rough summary statistic of the unconditional expected severity of default, but as will be seen, distributions of severity are both disperse and skewed and so inadequately described by even their conditional means. The results make it clear that the likelihood of default and the likely magnitude of that default must be treated as separate phenomena.     

Latent Markov Modelling of Recidivism Data

This article discusses the application of latent Markov modelling for the analysis of recidivism data. We briefly examine the relations of Markov modelling with log–linear analysis, pointing out pertinent differences as well. We show how the restrictive Markov model may be more easily applicable by adding latent variables to the model, in which case the latent Markov model is a dynamic version of the latent class model. As an illustration, we apply latent Markov analysis on an empirical data set of juvenile prosecution careers, showing how the Markov analyses producing well-fitting and interpretable solutions. We end by comparing the possible contributions of Markov modelling in recidivism research, outlining its drawbacks as well. Recommendations and directions for future research conclude the article.     

Variance-Gamma and Normal-Inverse Gaussian models: Goodness-of-fit to Chinese high-frequency index returns

In this study Variance-Gamma (VG) and Normal-Inverse Gaussian (NIG) distributions are compared with the benchmark of generalized hyperbolic distribution in terms of their fit to the empirical distribution of high-frequency stock market index returns in China. First, we estimate the considered models in a Markov regime switching framework for the identification of different volatility regimes. Second, the goodness-of-fit results are compared at different time scales of log-returns. Third, the goodness-of-fit results are validated through bootstrapping experiments. Our results show that as the time scale of log-returns decrease NIG model outperforms the VG model consistently and the difference between the goodness-of-fit statistics increase. For high-frequency Chinese index returns, NIG model is more robust and provides a better fit to the empirical distributions of returns at different time scales.     

Boekbesprekingen

Book reviewed in this article:
Analysing Qualitative Data, A. E. Maxwell, Methuen's monographs on applied probability and statistics, Methuen
Population Redistribution and Economic Growth, United States 1870–1950, E. S. Lee, A. R. Miller, C. P. Brainerd, R. A. Easterlin, S. Kuznets, D. S. Thomas
The Ownership of Major Consumer Durables, J. S. Cramer
Elementaire Statistiek, L. A. van Wijk, J. B. Wolters
Elementary Statistics, P. G. Hoel
The passage problem for a stationary Markov chain, J. H. B. Kemperman, Statistical Research Monographs, Volume 1
An introduction to Statistical Communication Theory, D. Middleton
An Introduction to infinitely many variates, E. A. Robinson, Grifin's Statistical Monographs and Courses
Queues, D. R. Cox en W. L. Smith, Methuen's Monographs on Applied Probability and Statistics