That's the only point where I have conceptual problems.
wls_results.scale is not a robust estimator, It could be far off if there are large outliers.
I think we should still start with mad, center=0