Simulating Correlated Hull–White Rates and Local-Volatility Equities
Summary
The document considers pricing a path-dependent option linked to a stock index and an interest-rate index. The proposed setup uses a Hull–White short-rate process and a local-volatility stock model, with correlated Brownian shocks. The question is whether rates and discount factors can be handled analytically at selected observation times while the stock is simulated numerically.
The accepted reply recommends deriving the rate process’s future mean and variance analytically, simulating rate paths, and numerically simulating the stock, for example with Euler–Maruyama. It stresses that the simulation must preserve the specified correlation, and that simulated rates can be used to form discount factors and discount option payoffs. The response outlines a practical division between analytical rate calculations and numerical stock-path generation, but supplies no derivation, scheme details, or validation. It also does not explain how to jointly sample correlated rates and stock increments, or address discretization bias and other accuracy considerations, so further implementation work is needed for a production pricer.
Key ideas
- Hull–White rates can be characterized analytically at future observation times.
- A local-volatility stock process can be simulated numerically.
- Correlated Brownian shocks must be reflected in the joint path simulation.
- Rate paths support discount-factor calculation and discounted payoff valuation.
- The reply gives an outline rather than a complete sampling algorithm or accuracy analysis.
Tags
Full text
# Analytic Hull White model with correlated stochastic processes
# Analytic Hull White model with correlated stochastic processes
I am trying to price a path dependent option which uses two underlyings (a stock index and an interest rate index). I am using Hull White model for interest rate modelling and local vol for stock index modelling. \begin{eqnarray} dr(t) &=& (\theta(t) - kr(t))dt + \sigma_r(t)dw_r(t)\\ \end{eqnarray}
\begin{eqnarray} ds(t) &=& r'(t)sdt + \sigma_s(t,s)sdw_s(t)\\ \end{eqnarray} Further, assume that $<dw_s, dw_r> = \rho dt$ and that $r'$ is deterministic. I need the samples for $s$ and $r$(and stochastic discount factors) only at certain times $t_1, t_2...t_n$. In case we did not have correlated processes, I could have calculated mean and variance of $r$ and discount factors analytically (and then creating the corresponding samples). Can we still solve this by calculating the mean and variance for $r$ and discount factors analytically and simulate $s_t$ numerically?
## Answer by Valter (score 2, accepted)
https://quant.stackexchange.com/a/77646
Yes, in the case of path-dependent option with two underlying assets, you can approach the problem by calculating the mean and variance for the interest rate $r(t)$ and corresponding discount factor analytically, while simulating the stock price $S(t)$ numerically. This will still work with correlated process since we have the correlation coefficient between the Brownian motions.
You can solve the Hull-White SDE analytically to find the mean and variance of the interest rate at future times $t_1,t_2,...,t_n$, then simulate the interest rate path.
For the local volatility model for the stock price, $r'(t)$ is deterministic, and $\sigma_s(t,s)$ is the local volatility model which is dependent on time and the stock price, you can typically simulate this process numerically. You can look into Euler–Maruyama method.
Then when simulating the paths for $S(t)$ and $r(t)$ just make sure you account for the correlation. Using the simulated paths for $r(t)$, you can calculate the discount factor for the respective time $t_1,t_2,...,t_n$.
Now that you have the analytical expression for the interest rate and the discount factor, and the numerical simulation setup for the stock price, you can generate the paths and compute the option's payoff at each point in time. Then you just discount these payoffs back to the present value.
## Answer by marie albertini (score 0)
https://quant.stackexchange.com/a/81416
have a look at the section 3.3 in https://arxiv.org/pdf/2411.05425
## Answer by marie albertini (score 0)
https://quant.stackexchange.com/a/81420
this is an article with a "simple" implementation answering the question asked in this post (too long to be posted here, but it is based on monotonic interpolation)...https://arxiv.org/pdf/2411.05425Shown in full with attribution under the source's licence. Licence: CC BY-SA 4.0 (Stack Exchange)
This summary was written by Stratmill's research agent from the original; it is not a copy of the source.