Commit b5573910 authored by Médéric Boquien's avatar Médéric Boquien
Browse files

The percentage of models with a low reduced χ² was too small by a factor 100.

parent a61694ec
......@@ -42,6 +42,7 @@
- Improve the sanitation of input data so that upper limits of extensive properties are not eliminated from the quantities to be fitted when upper limits are activated. (Médéric Boquien)
- Ensure that the radio module still works with recent versions of numpy. (Médéric Boquien & Laure Ciesla, reported by Wenjia Zhou)
- The unphysical wavy structure of the nebular continuum in the mid-infrared has been eliminated. (Médéric Boquien)
- The percentage of models with a low reduced χ² was too small by a factor 100. (Médéric Boquien, reported by Guang Yang)
### Optimised
- Slight speedup of the computation of the likelihood from the χ² using a multiplication rather than a division. (Médéric Boquien)
- Speedup of the computation of the χ² by ~10% taking the opposite of a scalar rather than of an array. (Médéric Boquien)
......@@ -322,9 +322,9 @@ class BestResultsManager(object):
chi2_red = self.chi2 / (nobs - 1)
# If low values of reduced chi^2, it means that the data are overfitted
# Errors might be under-estimated or not enough valid data.
print(f"{np.round((chi2_red < 1e-12).sum() / chi2_red.size, 1)}% of "
f"the objects have χ²_red~0 and "
f"{np.round((chi2_red < 0.5).sum() / chi2_red.size, 1)}% "
print(f"{np.round(100. * (chi2_red < 1e-12).sum() / chi2_red.size, 1)}%"
f" of the objects have χ²_red~0 and "
f"{np.round(100. * (chi2_red < 0.5).sum() / chi2_red.size, 1)}% "
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment