Change from -twoway function- to -twoway scatteri- for fit lines plotting. twoway function generates tempvars and unnecessarily sorts the entire dataset.
Bring in fastxtile version 1.22, which fixes a typo that causes crash when nquantiles > 100.  Patch to fastxtile contributed by Laszlo.
Remove percentiles() option