Update fit2.py

fixed the weighted fit
2020-10-13 18:07:47 +02:00
parent 6edc9c2ec4
commit fc70bdcb13
1 changed files with 1 additions and 1 deletions
--- a/pyzebra/fit2.py
+++ b/pyzebra/fit2.py
@ -157,7 +157,7 @@ def fitccl(
        ),
    )
    # the weighted fit
-    result = mod.fit(y, params, weights=y_err, x=x, calc_covar=True)
+    result = mod.fit(y, params, weights=[1/y_err[i] for i in range(len(y_err))], x=x, calc_covar=True)
    # u.ufloat to work with uncertanities
    fit_area = u.ufloat(result.params["g_amp"].value, result.params["g_amp"].stderr)
    comps = result.eval_components()