Generating Halton sequences using Mata

Abstract

This paper discusses the advantages of Halton sequences over pseudo-random uniform numbers when using simulation to approximate integrals numerically. We describe two types of sequences and give Mata examples. Finally, we document the Mata function halten(), currently in release 9.1 of Stata, which computes both a Halton sequence and its Hammersley variant. Options to use these point sets are available in the Stata 9 program asmprobit, a multinomial-probit estimator, and in the Stata 9.1 Mata function ghk(), the Geweke-HajivassiliouKeane multivariate-normal simulator. © 2006 StataCorp LP.