clear mata st_addobs(2500) st_addvar("double","x") st_addvar("double","y") st_addvar("double","z") k=1::2500 st_store(k,1,k) st_store(k,2,0.5*runiform(2500,1)) st_store(k,3,0.5:+0.5*halton(2500,1)) end twoway scatter y x, msize(tiny) color(blue) /// || scatter z x, msize(tiny) color(green) legend(off) xtitle("") /// title(Distribution: Van der Corput (top) vs pseudorandom) /// ylabel(, angle(0) format(%3.1f))