On Generating Monte Carlo Samples of Continuous Diffusion Bridges

Diffusion processes are widely used in engineering, finance, physics, and other fields. Usually continuous-time diffusion processes can be observed only at discrete time points. For many applications, it is often useful to impute continuous-time bridge samples that follow the diffusion dynamics and connect each pair of the consecutive observations. The sequential Monte Carlo (SMC) method is a useful tool for generating the intermediate paths of the bridge. The paths often are generated forward from the starting observation and forced in some ways to connect with the end observation. In this article we propose a constrained SMC algorithm with an effective resampling scheme guided by backward pilots carrying the information of the end observation. This resampling scheme can be easily combined with any forward SMC sampler. Two synthetic examples are used to demonstrate the effectiveness of the resampling scheme.

[1]  D. Talay.  Numerical solution of stochastic differential equations , 1994 .

[2]  Alan G. White,et al.  The Pricing of Options on Assets with Stochastic Volatilities , 1987 .

[3]  A. Pedersen,et al.  Consistency and asymptotic normality of an approximate maximum likelihood estimator for discretely observed diffusion processes , 1995 .

[4]  Tohru Ozaki,et al.  Estimation for nonlinear stochastic differential equations by a local linearization method 1 , 1998 .

[5]  Michael W. Brandt,et al.  Simulated Likelihood Estimation of Diffusions with an Application to Exchange Rate Dynamics in Incomplete Markets , 2001 .

[6]  Huu Tue Huynh,et al.  Solution of Stochastic Differential Equations , 2012 .

[7]  R. C. Merton,et al.  Option pricing when underlying stock returns are discontinuous , 1976 .

[8]  Lan Zhang,et al.  A Tale of Two Time Scales , 2003 .

[9]  N. Shephard,et al.  Econometric analysis of realized volatility and its use in estimating stochastic volatility models , 2002 .

[10]  G. Roberts,et al.  On inference for partially observed nonlinear diffusion models using the Metropolis–Hastings algorithm , 2001 .

[11]  Bjørn Eraker.  MCMC Analysis of Diffusion Models With Application to Finance , 2001 .

[12]  S. Ross,et al.  Option pricing: A simplified approach☆ , 1979 .

[13]  W. Greene,et al.  计量经济分析 = Econometric analysis , 2009 .

[14]  Jun S. Liu,et al.  Sequential Monte Carlo methods for dynamic systems , 1997 .

[15]  Jun S. Liu,et al.  Metropolized independent sampling with comparisons to rejection sampling and importance sampling , 1996, Stat. Comput..

[16]  N. Shephard,et al.  Likelihood INference for Discretely Observed Non-linear Diffusions , 2001 .

[17]  G. Roberts,et al.  MCMC methods for diffusion bridges , 2008 .

[18]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[19]  Yacine Aït-Sahalia,et al.  Disentangling diffusion from jumps , 2004 .

[20]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[21]  P. Fearnhead,et al.  Exact and computationally efficient likelihood‐based estimation for discretely observed diffusion processes (with discussion) , 2006 .

[22]  O. Stramer,et al.  On Simulated Likelihood of Discretely Observed Diffusion Processes and Comparison to Closed-Form Approximation , 2007 .

[23]  A. Gallant,et al.  Numerical Techniques for Maximum Likelihood Estimation of Continuous-Time Diffusion Processes , 2002 .

[24]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[25]  Jun S. Liu,et al.  Sequential Imputations and Bayesian Missing Data Problems , 1994 .