Quantitative Simulation American Options RFS (2001)

Longstaff-Schwartz Least Squares Monte Carlo (LSM) Lab

Value American-style derivative contracts via dynamic programming backward induction and cross-sectional OLS regression on in-the-money paths. Examine optimal stopping boundaries, early exercise premiums, and benchmark convergence.

Calibrated Regimes:
Longstaff-Schwartz (2001) Table 1 canonical benchmark: S0=36, K=40, r=6%, σ=20%, T=1.0y, 50 steps, American put.
American (LSM) Price
$4.472
Optimal Stopping
European Price
$3.844
Analytical B-S
Early Exercise Premium
$0.628
American − European
CRR Binomial Benchmark
$4.478
500-Step Tree Target
Monte Carlo SE
±$0.021
95% Confidence

Simulation Engine Controls

$
$
1.00
0.1y (~1 mo) 1.0y 3.0y
%
%
%

Numerical Discretization & Regression

Antithetic variates variance reduction is active across all runs.

Algorithm Diagnostics

Paths Exercised Early: 48.2%
CRR Binomial Tree Discrepancy: -$0.006
95% Confidence Interval: [$4.43, $4.51]
Time Step ($\Delta t$): 0.0200 y

Quantitative Analytics & Stopping Visualizer

Sensitivity Matrix: Strike ($K$) vs. Volatility ($\sigma$)

Early Exercise Premium (EEP in $) across moneyness and diffusion volatility.

Strike \ Vol σ = 10% σ = 20% σ = 30% σ = 40% σ = 50%

Monte Carlo Convergence: Paths ($M$) vs. Steps ($N$)

Simulation Standard Error (SE in ±$) illustrating asymptotic $\mathcal{O}(M^{-1/2})$ precision.

Steps \ Paths M = 1,000 M = 2,000 M = 4,000 M = 8,000 M = 16,000

Mathematical Foundations: Longstaff-Schwartz (2001) Algorithm

1. Optimal Stopping Time Formulation

American option valuation is governed by an optimal stopping problem over a complete probability space $(\Omega, \mathcal{F}, \mathbb{Q})$ with filtration $\{\mathcal{F}_t\}_{t=0}^T$:

$$V_0 = \sup_{\tau \in \mathcal{T}_{[0, T]}} \mathbb{E}^{\mathbb{Q}}\left[e^{-r\tau} h(S_\tau)\right]$$

where $\mathcal{T}_{[0, T]}$ is the set of all stopping times taking values in the discrete exercise grid $\{t_0, t_1, \dots, t_N\}$, and $h(S)$ represents the intrinsic exercise payoff $(K - S)^+$ for puts or $(S - K)^+$ for calls.

2. Cross-Sectional Conditional Expectation via OLS

At time step $t_k$, the continuation value $\Phi(S_{t_k})$ of holding the option until optimal stopping time $\tau > t_k$ satisfies the conditional expectation:

$$\Phi(S_{t_k}) = \mathbb{E}^{\mathbb{Q}}\left[e^{-r(t_{k+1}-t_k)} V_{t_{k+1}}(S_{t_{k+1}}) \;\middle|\; \mathcal{F}_{t_k}\right]$$

By the Markov property, conditional expectation depends solely on $S_{t_k}$. Longstaff and Schwartz approximate this expectation as a linear combination of $P$ orthogonal basis functions $\psi_p(S)$:

$$\Phi(S_{t_k}) \approx \sum_{p=0}^{P-1} \beta_p \psi_p(S_{t_k}), \quad \boldsymbol{\beta} = (\mathbf{X}^\top \mathbf{X})^{-1} \mathbf{X}^\top \mathbf{Y}$$

3. Restriction to In-the-Money (ITM) Trajectories

A key insight of Longstaff and Schwartz is restricting the cross-sectional regression at time $t_k$ strictly to paths that are in-the-money, i.e., $\{m \in \{1, \dots, M\} : h(S_{t_k}^{(m)}) > 0\}$.

Paths that are out-of-the-money have zero immediate exercise value, making early exercise strictly impossible at $t_k$. Omitting OTM paths eliminates noise and dramatically improves the numerical accuracy of $\boldsymbol{\beta}$ in the neighborhood of the critical exercise boundary $S^*(t)$.

4. Backward Induction & Stopping Rule

At expiration $t_N = T$, the cash flow for each path $m$ is initialized to $C_N^{(m)} = h(S_T^{(m)})$ with stopping time $\tau^{(m)} = T$. Moving backward from $t_{N-1}$ to $t_1$:

  1. Identify in-the-money paths at step $k$: $h(S_{t_k}^{(m)}) > 0$.
  2. Regress discounted realized future cash flows $Y_k^{(m)} = e^{-r(\tau^{(m)} - t_k)} C_{\tau^{(m)}}^{(m)}$ onto basis vectors $X_k^{(m)} = [\psi_0(S_{t_k}^{(m)}), \dots, \psi_{P-1}(S_{t_k}^{(m)})]$.
  3. Evaluate estimated continuation value: $\hat{\Phi}(S_{t_k}^{(m)}) = X_k^{(m)} \hat{\boldsymbol{\beta}}$.
  4. Optimal Exercise Decision: If $h(S_{t_k}^{(m)}) > \hat{\Phi}(S_{t_k}^{(m)})$, set $C_k^{(m)} = h(S_{t_k}^{(m)})$ and update $\tau^{(m)} = t_k$. Otherwise, leave $\tau^{(m)}$ unchanged.
$$\hat{V}_0 = \frac{1}{M} \sum_{m=1}^M e^{-r\tau^{(m)}} h(S_{\tau^{(m)}}^{(m)})$$

Interactive Mastery Quiz: Longstaff-Schwartz American Valuation

Test your mastery of dynamic programming backward induction, OLS basis regression, and optimal stopping times.

1. How does the Longstaff-Schwartz (2001) method overcome the traditional curse of dimensionality in American option simulation?

2. What is the exact decision rule for early exercise at time step $t_k$ in the backward induction step?

3. Why is the cross-sectional regression estimated solely on paths that are in-the-money ($h(S_t) > 0$)?

4. Which fundamental relationship always holds between American and European option prices on the same underlying?

Frequently Asked Questions

The Longstaff-Schwartz (2001) Least Squares Monte Carlo (LSM) algorithm is a seminal numerical technique for pricing American-style options by simulation. At each discrete exercise date moving backward from maturity, it estimates the unobservable continuation value of holding the option via cross-sectional ordinary least squares (OLS) regression against basis functions of the current asset price.

Restricting regression estimation to in-the-money (ITM) trajectories focuses numerical degrees of freedom exclusively where early exercise is economically feasible and rational. For out-of-the-money (OTM) paths, the immediate exercise payoff is zero, so early exercise is strictly suboptimal, and omitting them prevents bias in the fitted conditional expectation function.

Backward induction starts at expiration T where option value equals immediate intrinsic payoff. Moving backward step-by-step to t = t_{N-1}, t_{N-2}, ..., t_1, the algorithm compares immediate payoff h(S_t) against estimated continuation value C(S_t). If h(S_t) > C(S_t), stopping is optimal, the stopping time is updated to t, and cash flow is recorded; otherwise, continuation is preferred.

The Early Exercise Premium (EEP) represents the surplus economic value of an American option relative to its European counterpart: EEP = Price_American - Price_European. For puts, EEP is positive even without dividends due to the time value of money on strike cash. For calls, EEP is zero under zero dividends, but becomes strictly positive when continuous dividend yields exceed the risk-free interest rate.

Binomial trees like CRR are computationally ideal for 1D single-asset options but suffer from the curse of dimensionality when scaling to multi-asset baskets or path-dependent payoffs. The LSM Monte Carlo method scales polynomially with dimension, handling multi-factor stochastic processes, jump-diffusions, and cross-asset correlation matrices seamlessly.

Related American Valuation & Option Pricing Labs

Bjerksund-Stensland American Lab

Closed-form analytical approximations and flat early exercise boundary triggers for American calls and puts.

Binomial Tree Option Pricing Lab

Recombining lattice models across Cox-Ross-Rubinstein, Jarrow-Rudd, and Tian early exercise trees.

Fang-Oosterlee COS Pricing Lab

Exponential spectral convergence via Fourier-cosine series expansion of characteristic functions.

Carr-Madan FFT Pricing Lab

Simultaneous full strike spectrum valuation in $\mathcal{O}(N \log N)$ Fast Fourier Transform operations.