asteroid_spinprops 1.2.3__tar.gz → 1.2.5__tar.gz

{asteroid_spinprops-1.2.3 → asteroid_spinprops-1.2.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: asteroid_spinprops
-Version: 1.2.3
+Version: 1.2.5
 Summary: Collection of tools used for fitting sHG1G2 and SOCCA photometric models to sparse asteroid photometry
 License: MIT
 Author: Odysseas

{asteroid_spinprops-1.2.3 → asteroid_spinprops-1.2.5}/asteroid_spinprops/ssolib/periodest.py

@@ -210,7 +210,8 @@ def get_multiterm_period_estimate(
 
             rms = np.sqrt(np.mean(residuals**2))
             n_params = 2 * k + 1 + 3 * len(bands)
-            dof = len(residuals) - n_params
+            dof = n_params
+            N = len(residuals)
 
             results.append((k, f_best, rms, dof, n_params))
 
@@ -219,13 +220,13 @@ def get_multiterm_period_estimate(
         for i in range(len(results) - 1):
             k, f_best, rss, dof, n_params = results[i]
             k_next, f_best_next, rss_next, dof_next, n_params_next = results[i + 1]
-            F = ((rss - rss_next) / (dof - dof_next)) / (rss_next / dof_next)
+            F = ((rss**2/rss_next**2) - 1) * ((N - dof_next) / (dof_next - dof))
 
             # Here crit = Fstat value for which model_2 (more complex) is in fact better than model_1 (less complex)
             crit = ftest.ppf(
                 q=0.99,
                 dfn=n_params_next - n_params,
-                dfd=len(residuals) - n_params_next,
+                dfd=N - n_params_next,
             )
 
             model_comparison.loc[i, "k"] = k
@@ -233,7 +234,7 @@ def get_multiterm_period_estimate(
             model_comparison.loc[i, "f_best"] = f_best
             model_comparison.loc[i, "Fstat"] = F
             model_comparison.loc[i, "alpha_crit"] = crit
-            model_comparison.loc[i, "rms"] = rms
+            model_comparison.loc[i, "rms"] = rss
 
         cond = model_comparison["Fstat"] > model_comparison["alpha_crit"]
         model_comparison = model_comparison[
@@ -241,11 +242,11 @@ def get_multiterm_period_estimate(
         ]  # don't go for the more complex model
         f_chosen = model_comparison[
             model_comparison.k == model_comparison.k.min()
-        ].f_best[0]  # Simplest model
-        k_val = model_comparison[model_comparison.k == model_comparison.k.min()].k[
+        ].f_best.values[0]  # Simplest model
+        k_val = model_comparison[model_comparison.k == model_comparison.k.min()].k.values[
             0
         ]  # Simplest model
-        p_rms = model_comparison[model_comparison.k == model_comparison.k.min()].rms[
+        p_rms = model_comparison[model_comparison.k == model_comparison.k.min()].rms.values[
             0
         ]  # Simplest model
 

{asteroid_spinprops-1.2.3 → asteroid_spinprops-1.2.5}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "asteroid_spinprops"
-version = "1.2.3"
+version = "1.2.5"
 description = "Collection of tools used for fitting sHG1G2 and SOCCA photometric models to sparse asteroid photometry"
 authors = [
     {name = "Odysseas",email = "odysseas.xenos@proton.me"}