John W. Lau

Bayesian Model-Based Clustering Procedures

This article establishes a general formulation for Bayesian model-based clustering, in which subset labels are exchangeable, and items are also exchangeable, possibly up to covariate effects. The notational framework is rich enough to encompass a variety of existing procedures, including some recently discussed methods involving stochastic search or hierarchical clustering, but more importantly allows the formulation of clustering procedures that are optimal with respect to a specified loss function. Our focus is on loss functions based on pairwise coincidences, that is, whether pairs of items are clustered into the same subset or not. Optimization of the posterior expected loss function can be formulated as a binary integer programming problem, which can be readily solved by standard software when clustering a modest number of items, but quickly becomes impractical as problem scale increases. To combat this, a new heuristic item-swapping algorithm is introduced. This performs well in our numerical experiments, on both simulated and real data examples. The article includes a comparison of the statistical performance of the (approximate) optimal clustering with earlier methods that are model-based but ad hoc in their detailed definition.

Pricing options under a generalized Markov-Modulated jump-diffusion model

Abstract We consider the pricing of options when the dynamics of the risky underlying asset are driven by a Markov-modulated jump-diffusion model. We suppose that the market interest rate, the drift and the volatility of the underlying risky asset switch over time according to the state of an economy, which is modelled by a continuous-time Markov chain. The measure process is defined to be a generalized mixture of Poisson random measure and encompasses a general class of processes, for example, a generalized gamma process, which includes the weighted gamma process and the inverse Gaussian process. Another interesting feature of the measure process is that jump times and jump sizes can be correlated in general. The model considered here can provide market practitioners with flexibility in modelling the dynamics of the underlying risky asset. We employ the generalized regime-switching Esscher transform to determine an equivalent martingale measure in the incomplete market setting. A system of coupled partial-differential-integral equations satisfied by the European option prices is derived. We also derive a decomposition result for an American put option into its European counterpart and early exercise premium. Simulation results of the model have been presented and discussed.

Bayesian mixture of autoregressive models

An infinite mixture of autoregressive models is developed. The unknown parameters in the mixture autoregressive model follow a mixture distribution, which is governed by a Dirichlet process prior. One main feature of our approach is the generalization of a finite mixture model by having the number of components unspecified. A Bayesian sampling scheme based on a weighted Chinese restaurant process is proposed to generate partitions of observations. Using the partitions, Bayesian prediction, while accounting for possible model uncertainty, determining the most probable number of mixture components, clustering of time series and outlier detection in time series, can be done. Numerical results from simulated and real data are presented to illustrate the methodology.

Pricing participating products under a generalized jump-diffusion model

We propose a model for valuing participating life insurance products under a generalized jump-diffusion model with a Markov-switching compensator. It also nests a number of important and popular models in finance, including the classes of jump-diffusion models and Markovian regime-switching models. The Esscher transform is employed to determine an equivalent martingale measure. Simulation experiments are conducted to illustrate the practical implementation of the model and to highlight some features that can be obtained from our model.

Bayesian nonparametric estimation and consistency of mixed multinomial logit choice models

This paper develops nonparametric estimation for discrete choice models based on the Mixed Multinomial Logit (MMNL) model. It has been shown that MMNL models encompass all discrete choice models derived under the assumption of random utility maximization, subject to the identification of an unknown distribution G. Noting the mixture model description of the MMNL, we employ a Bayesian nonparametric approach, using nonparametric priors on the unknown mixing distribution G, to estimate the unknown choice probabilities. Theoretical support for the use of the proposed methodology is provided by establishing strong consistency of a general nonparametric prior on G under simple sufficient conditions. Consistency is defined according to a L1-type distance on the space of choice probabilities and is achieved by extending to a regression model framework a recent approach to strong consistency based on the summability of square roots of prior probabilities. Moving to estimation, slightly different techniques for non-panel and panel data models are discussed. For practical implementation, we describe efficient and relatively easy to use blocked Gibbs sampling procedures. A simulation study is also performed to illustrate the proposed methods and the exibility they achieve with respect to parametric Gaussian MMNL models.

On Bayesian mixture credibility

We introduce a class of Bayesian infinite mixture models first introduced by Lo (1984) to determine the credibility premium for a non-homogeneous insurance portfolio. The Bayesian infinite mixture models provide us with much flexibility in the specification of the claim distribution. We employ the sampling scheme based on a weighted Chinese restaurant process introduced in Lo et al. (1996) to estimate a Bayesian infinite mixture model from the claim data. The Bayesian sampling scheme also provides a systematic way to cluster the claim data. This can provide some insights into the risk characteristics of the policyholders. The estimated credibility premium from the Bayesian infinite mixture model can be written as a linear combination of the prior estimate and the sample mean of the claim data. Estimation results for the Bayesian mixture credibility premiums will be presented.

Modelling long-term investment returns via Bayesian infinite mixture time series models

This paper introduces the class of Bayesian infinite mixture time series models first proposed in Lau & So (2004) for modelling long-term investment returns. It is a flexible class of time series models and provides a flexible way to incorporate full information contained in all autoregressive components with various orders by utilizing the idea of Bayesian averaging or mixing. We adopt a Bayesian sampling scheme based on a weighted Chinese restaurant process for generating partitions of investment returns to estimate the Bayesian infinite mixture time series models. Instead of using the point estimates, as in the classical or non-Bayesian approach, the estimation in this paper is performed by the full Bayesian approach, utilizing the idea of Bayesian averaging to incorporate all information contained in the posterior distributions of the random parameters. This provides a natural way to incorporate model risk or uncertainty. The proposed models can also be used to perform clustering of investment returns and detect outliers of returns. We employ the monthly data from the Toronto Stock Exchange 300 (TSE 300) indices to illustrate the implementation of our models and compare the simulated results from the estimated models with the empirical characteristics of the TSE 300 data. We apply the Bayesian predictive distribution of the logarithmic returns obtained by the Bayesian averaging or mixing to evaluate the quantile-based and conditional tail expectation risk measures for segregated fund contracts via stochastic simulation. We compare the risk measures evaluated from our models with those from some well-known and important models in the literature, and highlight some features that can be obtained from our models.

Viterbi-Based Estimation for Markov Switching GARCH Model

Abstract We outline a two-stage estimation method for a Markov-switching Generalized Autoregressive Conditional Heteroscedastic (GARCH) model modulated by a hidden Markov chain. The first stage involves the estimation of a hidden Markov chain using the Vitberi algorithm given the model parameters. The second stage uses the maximum likelihood method to estimate the model parameters given the estimated hidden Markov chain. Applications to financial risk management are discussed through simulated data.