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?