Qm
Advanced

Brownian Bridge And Stratified Sampling

Simulate the end of a path before its middle, and the biggest, most important random draws land first, exactly where a low-discrepancy sequence works best.

Prerequisites: Monte Carlo Option Pricing, Quasi-Monte Carlo And Sobol Sequences

Ask someone to sketch a stock's path over a year and they usually draw it left to right, day by day. But if all you actually need is where the path ends, which is all a European option cares about, drawing it day by day is backwards. It is like painting a portrait by filling in one square millimetre at a time in reading order: technically it works, but the shape of the face isn't visible until the very last stroke. A better painter blocks in the eyes and the jawline first, then fills in the detail. Simulating a Monte Carlo path works the same way, and the standard approach, draw one random step, add it, draw the next, add it, repeat hundreds of times, is the millimetre-by-millimetre painter. Every one of those steps consumes a fresh random number, and by the time you finally know where the path ends, hundreds of draws have already been spent.

Paint the big strokes first

A Brownian bridge is exactly that reordering. Instead of simulating a path step by step from today to the end, it simulates the endpoint first, then the midpoint conditional on the endpoint, then the quarter-points conditional on their neighbours, and so on, refining the path from coarse structure down to fine jitter. The mathematics behind each step is simple: given the value at two known times, the value at a time in between is normally distributed around the straight line connecting them, with a variance that depends only on how far that midpoint sits from each end,

WtWs,Wu    N ⁣(Ws+tsus(WuWs),  (ts)(ut)us),s<t<u.W_t \mid W_s, W_u \;\sim\; \mathcal{N}\!\left(W_s + \frac{t-s}{u-s}(W_u - W_s),\; \frac{(t-s)(u-t)}{u-s}\right), \qquad s < t < u.

In plain English: the bridge's value at time tt is a weighted average of its two known neighbours (closer to whichever one is nearer), plus a random wiggle whose size is largest exactly in the middle and shrinks to zero at the two ends, since the endpoints are already fixed.

Why this matters for quasi-Monte Carlo

A quasi-Monte Carlo path with 250250 time steps uses 250250 random numbers, one per dimension of a low-discrepancy sequence. The Sobol sequence's low-discrepancy property is strongest in its first handful of dimensions and degrades in later ones, a Sobol sequence's 200th coordinate is barely better spread than an ordinary random number. A step-by-step path spreads the payoff's dependence evenly across all 250250 dimensions, so Sobol gets no advantage there. A Brownian-bridge path concentrates the large moves, the endpoint, then the midpoint, into the first two or three dimensions, and dumps the fine, low-impact jitter into the later, weaker dimensions. The payoff ends up depending mostly on a handful of well-covered dimensions, which is what makes quasi-Monte Carlo actually pay off on multi-step, path-dependent options.

Worked example 1: building one bridge path with 4 points

Simulate a Brownian path with W0=0W_0 = 0, over [0,1][0,1], at times 0,0.25,0.5,0.75,10, 0.25, 0.5, 0.75, 1, using standard normal draws z1,z2,z3=1.0,0.4,0.6z_1, z_2, z_3 = 1.0, 0.4, -0.6 in bridge order.

  1. Endpoint (t=1t=1): W1=z11=1.00W_1 = z_1 \sqrt{1} = 1.00.
  2. Midpoint (t=0.5t=0.5, between W0=0W_0=0 and W1=1.00W_1=1.00): mean =0.5(1.00)=0.50= 0.5(1.00) = 0.50, variance =(0.5)(0.5)/1=0.25= (0.5)(0.5)/1 = 0.25, so std =0.5=0.5. W0.5=0.50+0.5(0.4)=0.70W_{0.5} = 0.50 + 0.5(0.4) = 0.70.
  3. Quarter point (t=0.25t=0.25, between W0=0W_0=0 and W0.5=0.70W_{0.5}=0.70): mean =0.5(0.70)=0.35=0.5(0.70)=0.35, variance =(0.25)(0.25)/0.5=0.125=(0.25)(0.25)/0.5=0.125, std =0.354=0.354. W0.25=0.35+0.354(0.6)=0.138W_{0.25} = 0.35 + 0.354(-0.6) = 0.138.

Notice only three random numbers built four path values, and the first one already pins down where the path ends up, that is the whole point. Compare this with ordinary step-by-step simulation, which would draw four separate increments, ΔW1,ΔW2,ΔW3,ΔW4\Delta W_1, \Delta W_2, \Delta W_3, \Delta W_4, one per quarter-step, and only learn the endpoint after adding up all four. Both approaches produce a path with exactly the same statistical law, the bridge is not a different model, only a different order of construction, but the bridge tells you the single most decision-relevant number, the endpoint, using only its first random draw.

Worked example 2: stratified sampling on the same endpoint

Stratified sampling is the cousin idea applied to the draw itself: instead of drawing the endpoint's normal variate freely, split the [0,1][0,1] probability range into NN equal strata and force exactly one draw into each, so the tails are never accidentally under-sampled. With N=4N=4 strata for the endpoint W1W_1: draw uniforms from [0,0.25),[0.25,0.5),[0.5,0.75),[0.75,1)[0, 0.25), [0.25, 0.5), [0.5, 0.75), [0.75, 1), say 0.10,0.40,0.60,0.900.10, 0.40, 0.60, 0.90, and invert each through the normal CDF to get z=1.28,0.25,0.25,1.28z = -1.28, -0.25, 0.25, 1.28. Compare a plain random draw of 4 uniforms, which might by chance land at 0.48,0.51,0.55,0.600.48, 0.51, 0.55, 0.60 (all clustered near the median, missing both tails entirely, giving z0.05,0.03,0.13,0.25z \approx -0.05, 0.03, 0.13, 0.25). The stratified draw guarantees representation of extreme outcomes, the very paths that determine an out-of-the-money option's value, while the naive draw can miss them by pure chance.

time t=0 1: end 2: mid 3: quarter 4: 3/4 pt
Numbers show the order points are simulated in, not their position on the path. The endpoint and midpoint (large, structural moves) are drawn first with fresh random numbers; the quarter-points (fine detail) are drawn last, conditional on their neighbours.

The paths below are built the ordinary step-by-step way, not as a bridge, but they show the raw material a bridge reorders: watch how much the path wanders between any two points, and imagine simulating the right-hand endpoint of each path first, then filling in the wandering detail afterward.

Path explorer
13055time →
end (bold path) 100.38spread of ends 58.966 independent paths, same settings

What this means in practice

Any desk pricing path-dependent products, barriers, Asians, cliquets, by Monte Carlo builds paths with a Brownian bridge whenever it also wants low-discrepancy sequences, and combines it with other variance reduction like control variates. The gain compounds: bridge construction concentrates the payoff's sensitivity into a few dimensions, and Sobol covers those few dimensions extremely well. The same reordering trick also helps plain Monte Carlo run diagnostics: because the endpoint is known after the very first random draw, a desk can price a rough, first-pass estimate of a book using far fewer total draws than a full step-by-step simulation would need, refining the path detail only for the positions that turn out to matter.

A Brownian bridge changes the order random numbers are consumed, not the distribution of the path, the simulated stock still has the same law, and standard step-by-step simulation and bridge simulation give statistically identical prices under plain Monte Carlo. The mistake is expecting the bridge alone to reduce variance; it only helps once it is paired with a low-discrepancy or stratified sampler that benefits from concentrating structure in early draws. Used with plain pseudo-random numbers, a Brownian bridge buys nothing.

Simulating the endpoint first and the fine jitter last packs a path's important randomness into the first few draws, which is precisely where quasi-Monte Carlo and stratified sampling do their best work.

Discussion

Sign in to join the discussion · reading is open to everyone

💡 Discussion rules

  1. Ask and answer about this concept. Off-topic gets removed.
  2. No homework dumps. Show what you tried first.
  3. Corrections are welcome. Cite a source when you claim an error.

Loading discussion…

Related concepts

Practice in interviews

Further reading

  • Glasserman, Monte Carlo Methods in Financial Engineering (Ch. 3, 5)
  • Caflisch, Morokoff & Owen (1997), Valuation of Mortgage Backed Securities Using Brownian Bridges
ShareTwitterLinkedIn