RosettaCodeData/Task/Van-der-Corput-sequence/Stata/van-der-corput-sequence-2.stata

16 lines
419 B
Text
Raw Permalink Normal View History

2017-09-23 10:01:46 +02:00
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) ///
2018-06-22 20:57:24 +00:00
|| scatter z x, msize(tiny) color(green) legend(off) xtitle("") ///
title(Distribution: Van der Corput (top) vs pseudorandom) ///
ylabel(, angle(0) format(%3.1f))