How to apply quasi-Monte Carlo to path-dependent options?

Following up on my recent question on variance reduction in a Cox-Ingersoll-Ross Monte Carlo simulation, I would like to learn more about using a quasi-random sequence, such as Sobol or Niederreiter, to generate quasi-Monte Carlo interest rate paths.

As before, the objective is to estimate the value of a path-dependent interest rate option. I am evaluating the integral defined by my valuation equation independently for $N$ paths, each of which are $T$ periods long. $T$ may be as high as 40 (quarters, or 10 years) or more in some cases. Should I use the first $N$ $T$-dimensional sequences? The particular Niederreiter implementation I found only allows up to 20 dimensions, which is too small for most securities in my sample. Is this a limitation of the particular implementation or a general limitation of the Niederreiter algorithm? Is it valid to append successive 20-dimensional sequences (with a different random seed each time)?

I've also seen some applications where the first 1000 or so numbers in the sequence are discarded, or where successive innovations to the time-series skip as many as 100 numbers in the sequence between time steps. What is the basis for these modifications to the basic sequences, and does it improve convergence in practical applications?

Finally, is there any research on the pros and cons of the various low-discrepancy sequence algorithms for financial applications? Why does Brian B favor Niederreiter, for example?

Answer

For such high-dimensional path problems you will want to use the Morokov technique (you can find the paper online), which takes QR samples for the "important" dimensions and then reverts to pseudorandom for the less important dimensions in an interest rate problem remarkably similar to yours. (Similar principles apply to using QR sequences in factor model based simulations). Effectively, each sample in $T$-dimensional space has $P$ of its dimensions from the QR sequence and $T-P$ dimensions from a pseudorandom sequence.

Skipping techniques are generally considered a good idea, though the need for them is mitigated by the Morokov ordering. The general principle is to start-skip about as many sample points (in $T$ dimensional space) as you plan to take overall. Intermediate skips are, in my opinion, pointless when Morokov ordering is used.

Do not forget to form any normal samples from the uniform QR samples using Moro inversion rather than Box-Mueller.

As for sequence preferences, my firm did a study on a sample set of equity derivatives and found Niederreiter sequences to perform better than Halton or Sobol, by a slight margin.