search

sczesla / PyAstronomy

A collection of astronomy-related routines in Python
152 stars 35 forks source link

Generalized ESD #45

Closed jcdang closed 5 years ago

jcdang commented 5 years ago

If you want to match NIST's R implementation,

https://github.com/sczesla/PyAstronomy/blob/595f1eaf32f99493931d72404f7e0877e88ca9e3/src/pyasl/asl/outlier.py#L62

should be

    xstd = xm.std(ddof=1)
sczesla commented 5 years ago

Thanks. That is true. I added an option to enable consistent behavior. I did not want to change existing default behavior, which is not wrong either as far as I can tell. Stefan