Efficient and flexible model-based clustering of jumps in diffusion processes

Jump–diffusion processes involving diffusion processes with discontinuous movements, called jumps, are widely used to model time-series data that commonly exhibit discontinuity in their sample paths. The existing jump–diffusion models have been recently extended to multivariate time-series data. The models are, however, still limited by a single parametric jump-size distribution that is common across different subjects. Such strong parametric assumptions for the shape and structure of a jump-size distribution may be too restrictive and unrealistic for multiple subjects with different characteristics. This paper thus proposes an efficient Bayesian nonparametric method to flexibly model a jump-size distribution while borrowing information across subjects in a clustering procedure using a nested Dirichlet process. For efficient posterior computation, a partially collapsed Gibbs sampler is devised to fit the proposed model. The proposed methodology is illustrated through a simulation study and an application to daily stock price data for companies in the S&P 100 index from June 2007 to June 2017.     

Diffusion approximations for insurance risk processes

We consider a risk-reserve process for an insurance company where premium income and the claim sum process are modeled as a renewal reward processes. Moreover, dividends are paid out according to a barrier rule. The aim of the article is to establish a diffusion approximation of this model and to compute ruin probabilities (in finite and in infinite time) and other relevant statistics approximately using the limiting diffusion process. We also demonstrate that, under special circumstances, there exists a stationary distribution for the limiting diffusion.     

Testing for Threshold Diffusion

The threshold diffusion model assumes a piecewise linear drift term and a piecewise smooth diffusion term, which constitutes a rich model for analyzing nonlinear continuous-time processes. We consider the problem of testing for threshold nonlinearity in the drift term. We do this by developing a quasi-likelihood test derived under the working assumption of a constant diffusion term, which circumvents the problem of generally unknown functional form for the diffusion term. The test is first developed for testing for one threshold at which the drift term breaks into two linear functions. We show that under some mild regularity conditions, the asymptotic null distribution of the proposed test statistic is given by the distribution of certain functional of some centered Gaussian process. We develop a computationally efficient method for calibrating the p-value of the test statistic by bootstrapping its asymptotic null distribution. The local power function is also derived, which establishes the consistency of the proposed test. The test is then extended to testing for multiple thresholds. We demonstrate the efficacy of the proposed test by simulations. Using the proposed test, we examine the evidence of nonlinearity in the term structure of a long time series of U.S. interest rates.     

Ancestral Inference in Population Genetics Models with Selection (with Discussion)

A new algorithm is presented for exact simulation from the conditional distribution of the genealogical history of a sample, given the composition of the sample, for population genetics models with general diploid selection. The method applies to the usual diffusion approximation of evolution at a single locus, in a randomly mating population of constant size, for mutation models in which the distribution of the type of a mutant does not depend on the type of the progenitor allele; this includes any model with only two alleles. The new method is applied to ancestral inference for the two‐allele case, both with genic selection and heterozygote advantage and disadvantage, where one of the alleles is assumed to have resulted from a unique mutation event. The paper describes how the method could be used for inference when data are also available at neutral markers linked to the locus under selection. It also informally describes and constructs the non‐neutral Fleming–Viot measure‐valued diffusion.     

A marginalized diffusion model for estimating age at first lower endoscopy use from current status data

Summary.  We propose an approach for estimating the age at first lower endoscopy examination from current status data that were collected via two series of cross-sectional surveys. To model the national probability of ever having a lower endoscopy examination, we incorporate birth cohort effects into a mixed influence diffusion model. We link a state-specific model to the national level diffusion model by using a marginalized modelling approach. In future research, results from our model will be used as microsimulation model inputs to estimate the contribution of endoscopy examinations to observed changes in colorectal cancer incidence and mortality.     

Markov-Switching and Stochastic Volatility Diffusion Models of Short-Term Interest Rates

This article empirically compares the Markov-switching and stochastic volatility diffusion models of the short rate. The evidence supports the Markov-switching diffusion model. Estimates of the elasticity of volatility parameter for single-regime models unanimously indicate an explosive volatility process, whereas the Markov-switching models estimates are reasonable. Itis found that either Markov switching or stochastic volatility, but not both, is needed to adequately fit the data. A robust conclusion is that volatility depends on the level of the short rate. Finally, the Markov-switching model is the best for forecasting. A technical contribution of this article is a presentation of quasi-maximum likelihood estimation techniques for the Markov-switching stochastic-volatility model.     

高频数据驱动的连续时间模型估计

现代金融经济学中连续时间模型能够更方便地描述重要经济变量的动态过程如股价、汇率和利率等。为连续时间模型提出了一种高频数据驱动的二阶段估计方法,增强了连续时间扩展模型的弹性和可操作性。以Vasicek模型为例给出了该方法的应用实例,首先在第一阶段使用实现波动率方法估计出模型的扩散项参数,然后使用实际数据的稳态分布的前向方程估计漂移项参数。此方法对模型初始设定和优化算法依赖程度低,结果较为稳定可靠。     

Likelihood-based inference for a class of multivariate diffusions with unobserved paths

This paper presents a Markov chain Monte Carlo algorithm for a class of multivariate diffusion models with unobserved paths. This class is of high practical interest as it includes most diffusion driven stochastic volatility models. The algorithm is based on a data augmentation scheme where the paths are treated as missing data. However, unless these paths are transformed so that the dominating measure is independent of any parameters, the algorithm becomes reducible. The methodology developed in Roberts and Stramer [2001a. On inference for partial observed nonlinear diffusion models using the metropolis-hastings algorithm. Biometrika 88(3); 603–621] circumvents the problem for scalar diffusions. We extend this framework to the class of models of this paper by introducing an appropriate reparametrisation of the likelihood that can be used to construct an irreducible data augmentation scheme. Practical implementation issues are considered and the methodology is applied to simulated data from the Heston model.     

Estimation of a change point in the hazard rate of Lindley model under right censoring

In this article, we consider a single change point model for a sudden change in the hazard rate of Lindley distribution to model right-censored survival data. We derive the quantile function to generate random numbers from the proposed distribution by using the Lambert function. The maximum likelihood estimation method is used to estimate parameters of the change point model. A simulation study is also carried out to analyze the performance of the estimators. To validate our findings, a dataset on bone marrow transplant for patients of acute lymphoblastic leukemia is analyzed using the proposed model and is compared with the existing exponential single change point model.     

Spectral Estimation of Covolatility from Noisy Observations Using Local Weights

We propose localized spectral estimators for the quadratic covariation and the spot covolatility of diffusion processes, which are observed discretely with additive observation noise. The appropriate estimation for time‐varying volatilities is based on an asymptotic equivalence of the underlying statistical model to a white‐noise model with correlation and volatility processes being constant over small time intervals. The asymptotic equivalence of the continuous‐time and discrete‐time experiments is proved by a construction with linear interpolation in one direction and local means for the other. The new estimator outperforms earlier non‐parametric methods in the literature for the considered model. We investigate its finite sample size characteristics in simulations and draw a comparison between various proposed methods.     

Pricing options with credit risk in a reduced form model

This article investigates the valuation of European option with credit risk in a reduced form model. We assume that the interest rate follows the Vasicek model and the intensity of default is driven by a jump diffusion process. We obtain the closed form formula for the price of the option and provide some numerical illustrations of the results obtained.     

Estimation of the multi-criteria worths of the alternatives in a hierarchical structure of comparisons

The aggregated worths of the alternatives, when compared with respect to several criteria, are estimated in a hierarchical comparisons model introduced by Saaty (1980). A multiplicative model is used for the paired comparisons data which are collected in a ratio scale in this set-up in any level of this hierarchy. An iterative scheme is found for the maximum likelihood estimation of the worth parameters in this multiplicative model. The iterative values are shown to be convergent monotonically to the estimates. We also obtain the asymptotic dispersion matrix of the maximum likelihood estimates of the relative worths of the alternatives according to a single criterion as well as those according to the over-all suitability when compared under several criteria. A numerical example is presented to illustrate the method developed in this paper. Simulation techniques are employed to find the average number of iterations required for the convergence of the above iterative scheme.     

A solvable singular control problem driven by a jump diffusion process with applications

We consider a class of singular control problems driven by a double exponential jump diffusion process, which come from the reversible investment problem. In some interesting cases (e.g., the running cost function is given by the so-called Cobb-Douglas production function), we give the explicit solutions to the singular control problem by using the connection between singular control and optimal switching. We solve a collection of consistent optimal switching problems and yield the explicit solution for the singular control problem. We then give an application to a particular inventory control problem in a single random period.     

Bayesian diffusion process models with time-varying parameters

The diffusion process is a widely used statistical model for many natural dynamic phenomena but its inference is very complicated because complete data describing the diffusion sample path is not necessarily available. In addition, data is often collected with substantial uncertainty and it is not uncommon to have missing observations. Thus, the observed process will be discrete over a finite time period and the marginal likelihood given by this discrete data is not always available. In this paper, we consider a class of nonstationary diffusion process models with not only the measurement error but also discretely time-varying parameters which are modeled via a state space model. Hierarchical Bayesian inference for such a diffusion process model with time-varying parameters is applied to financial data.     

A spatial Bayesian semiparametric mixture model for positive definite matrices with applications in diffusion tensor imaging

Studies on diffusion tensor imaging (DTI) quantify the diffusion of water molecules in a brain voxel using an estimated 3 × 3 symmetric positive definite (p.d.) diffusion tensor matrix. Due to the challenges associated with modelling matrix‐variate responses, the voxel‐level DTI data are usually summarized by univariate quantities, such as fractional anisotropy. This approach leads to evident loss of information. Furthermore, DTI analyses often ignore the spatial association among neighbouring voxels, leading to imprecise estimates. Although the spatial modelling literature is rich, modelling spatially dependent p.d. matrices is challenging. To mitigate these issues, we propose a matrix‐variate Bayesian semiparametric mixture model, where the p.d. matrices are distributed as a mixture of inverse Wishart distributions, with the spatial dependence captured by a Markov model for the mixture component labels. Related Bayesian computing is facilitated by conjugacy results and use of the double Metropolis–Hastings algorithm. Our simulation study shows that the proposed method is more powerful than competing non‐spatial methods. We also apply our method to investigate the effect of cocaine use on brain microstructure. By extending spatial statistics to matrix‐variate data, we contribute to providing a novel and computationally tractable inferential tool for DTI analysis.     

Rate estimation in partially observed Markov jump processes with measurement errors

We present a simulation methodology for Bayesian estimation of rate parameters in Markov jump processes arising for example in stochastic kinetic models. To handle the problem of missing components and measurement errors in observed data, we embed the Markov jump process into the framework of a general state space model. We do not use diffusion approximations. Markov chain Monte Carlo and particle filter type algorithms are introduced which allow sampling from the posterior distribution of the rate parameters and the Markov jump process also in data-poor scenarios. The algorithms are illustrated by applying them to rate estimation in a model for prokaryotic auto-regulation and the stochastic Oregonator, respectively.     

Optimal financing and dividend policy with Markovian switching regimes

This work investigates an optimal financing and dividend problem for an insurer whose surplus process is modulated by an observable continuous-time and finite-state Markov chain. We assume that the insurer should never go bankrupt by issuing new equity. The goal of the insurer is to maximize the expected present value of the dividends payout minus the discounted cost of equity issuance. We obtain the optimal policies and explicit expressions for the value functions when the risk reserve process is modeled by both upward jump model and its diffusion approximation. Numerical illustrations of the sensitivities of the model parameters are provided.     

An empirical likelihood goodness-of-fit test for time series

Summary. Standard goodness-of-fit tests for a parametric regression model against a series of nonparametric alternatives are based on residuals arising from a fitted model. When a parametric regression model is compared with a nonparametric model, goodness-of-fit testing can be naturally approached by evaluating the likelihood of the parametric model within a nonparametric framework. We employ the empirical likelihood for an α -mixing process to formulate a test statistic that measures the goodness of fit of a parametric regression model. The technique is based on a comparison with kernel smoothing estimators. The empirical likelihood formulation of the test has two attractive features. One is its automatic consideration of the variation that is associated with the nonparametric fit due to empirical likelihood's ability to Studentize internally. The other is that the asymptotic distribution of the test statistic is free of unknown parameters, avoiding plug-in estimation. We apply the test to a discretized diffusion model which has recently been considered in financial market analysis.     

Optimum Experimental Designs for Choosing Between Competitive and Non Competitive Models of Enzyme Inhibition

In the study of the inhibition of enzyme reactions the Michaelis–Menten model is extended to include two experimental variables and three or more parameters. We combine the three-parameter models for competitive and non competitive inhibition in a four-parameter model and use optimum design theory to find D- and Ds-optimum designs for discriminating between the models. These designs are compared with compound T-optimum designs which provide the most powerful tests for discrimination between models. A single design is found with high discrimination efficiency whichever model is true.     

Optimal Strategies for a Long-Term Static Investor

The optimal strategies for a long-term static investor are studied. Given a portfolio of a stock and a bond, we derive the optimal allocation of the capitals to maximize the expected long-term growth rate of a utility function of the wealth. When the bond has a constant interest rate, three models for the underlying stock price processes are studied: Heston model, 3/2 model, and jump diffusion model. We also study the optimal strategies for a portfolio in which the stock price process follows a Black-Scholes model and the bond process has a Vasicek interest rate that is correlated to the stock price.