asteroid_spinprops 1.2.2__tar.gz → 1.2.4__tar.gz

{asteroid_spinprops-1.2.2 → asteroid_spinprops-1.2.4}/PKG-INFO

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

{asteroid_spinprops-1.2.2 → asteroid_spinprops-1.2.4}/asteroid_spinprops/ssolib/modelfit.py

@@ -225,14 +225,14 @@ def get_fit_params(
                 )
 
                 for ra, dec in zip(ra_init, dec_init):
-                    for period_in in period_scan:
+                    for period_sc in period_scan:
                         p_in = [
                             shg1g2_params[H_key],
                             0.15,
                             0.15,  # G1,2
                             np.radians(ra),
                             np.radians(dec),
-                            period_in,  # in days
+                            period_sc,  # in days
                             shg1g2_params["a_b"],
                             shg1g2_params["a_c"],
                             0.1,
@@ -268,8 +268,9 @@ def get_fit_params(
                             SOCCA_opt["Nbs"] = Nbs
                         except Exception:
                             SOCCA_opt["Period_class"] = -1 # Classification error
-                    SOCCA_opt["prms"] = p_rms
-                    SOCCA_opt["k_terms"] = k_val
+                    if period_in is None:
+                        SOCCA_opt["prms"] = p_rms
+                        SOCCA_opt["k_terms"] = k_val
                 except Exception:
                     SOCCA_opt = {"Failed at SOCCA inversion": 1}
                 return SOCCA_opt

{asteroid_spinprops-1.2.2 → asteroid_spinprops-1.2.4}/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[

{asteroid_spinprops-1.2.2 → asteroid_spinprops-1.2.4}/pyproject.toml

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