huff 1.5.3__tar.gz → 1.5.4__tar.gz

{huff-1.5.3 → huff-1.5.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: huff
-Version: 1.5.3
+Version: 1.5.4
 Summary: huff: Huff Model Market Area Analysis
 Author: Thomas Wieland
 Author-email: geowieland@googlemail.com
@@ -18,12 +18,10 @@ Thomas Wieland [ORCID](https://orcid.org/0000-0001-5168-9846) [EMail](mailto:geo
 See the /tests directory for usage examples of most of the included functions.
 
 
-## Updates v1.5.3
+## Updates v1.5.4
 - Bugfixes:
-  - InteractionMatrix.summary() and HuffModel.summary(): No KeyError when param is None anymore
-- Extensions:
-  - Confidence intervals (bootstrap) for Huff Model ML estimations
-  - ValueError and KeyError are more precise (function is included)
+  - Use of check_vars() is now optional (default: True)
+  - Correction of args argument when calling minimize() in huff_ml_fit() 
   
 
 ## Features

{huff-1.5.3 → huff-1.5.4}/README.md

@@ -10,12 +10,10 @@ Thomas Wieland [ORCID](https://orcid.org/0000-0001-5168-9846) [EMail](mailto:geo
 See the /tests directory for usage examples of most of the included functions.
 
 
-## Updates v1.5.3
+## Updates v1.5.4
 - Bugfixes:
-  - InteractionMatrix.summary() and HuffModel.summary(): No KeyError when param is None anymore
-- Extensions:
-  - Confidence intervals (bootstrap) for Huff Model ML estimations
-  - ValueError and KeyError are more precise (function is included)
+  - Use of check_vars() is now optional (default: True)
+  - Correction of args argument when calling minimize() in huff_ml_fit() 
   
 
 ## Features

{huff-1.5.3 → huff-1.5.4}/huff/models.py

@@ -4,8 +4,8 @@
 # Author:      Thomas Wieland 
 #              ORCID: 0000-0001-5168-9846
 #              mail: geowieland@googlemail.com              
-# Version:     1.5.3
-# Last update: 2025-07-15 17:22
+# Version:     1.5.4
+# Last update: 2025-07-18 18:07
 # Copyright (c) 2025 Thomas Wieland
 #-----------------------------------------------------------------------
 
@@ -766,7 +766,10 @@ class InteractionMatrix:
         self.supply_locations.metadata = supply_locations_metadata
         self.customer_origins.metadata = customer_origins_metadata
 
-    def utility(self):
+    def utility(
+        self,
+        check_df_vars: bool = True
+        ):
         
         interaction_matrix_df = self.interaction_matrix_df
 
@@ -781,10 +784,11 @@ class InteractionMatrix:
         if interaction_matrix_df["A_j"].isna().all():
             raise ValueError ("Error in utility calculation: Attraction variable is not defined")
 
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["A_j", "t_ij"]
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["A_j", "t_ij"]
+                )
         
         customer_origins = self.customer_origins
         customer_origins_metadata = customer_origins.get_metadata()
@@ -793,9 +797,9 @@ class InteractionMatrix:
         if tc_weighting["func"] == "power":
             interaction_matrix_df["t_ij_weighted"] = interaction_matrix_df["t_ij"] ** tc_weighting["param"]
         elif tc_weighting["func"] == "exponential":
-            interaction_matrix_df["t_ij_weighted"] = np.exp(tc_weighting["param"] * interaction_matrix_df['t_ij'])
+            interaction_matrix_df["t_ij_weighted"] = np.exp(tc_weighting["param"] * interaction_matrix_df["t_ij"])
         elif tc_weighting["func"] == "logistic":
-            interaction_matrix_df["t_ij_weighted"] = 1+np.exp(tc_weighting["param"][0] + tc_weighting["param"][1] * interaction_matrix_df['t_ij'])
+            interaction_matrix_df["t_ij_weighted"] = 1+np.exp(tc_weighting["param"][0] + tc_weighting["param"][1] * interaction_matrix_df["t_ij"])
         else:
             raise ValueError ("Error in utility calculation: Transport costs weighting is not defined.")
                            
@@ -818,8 +822,8 @@ class InteractionMatrix:
             
             for key, attrac_var in enumerate(attrac_vars):
                 
-                attrac_var_key = key #+1
-                if attrac_var_key == 0: #1:
+                attrac_var_key = key
+                if attrac_var_key == 0:
                     continue
                 
                 name = supply_locations_metadata["weighting"][attrac_var_key]["name"]
@@ -850,7 +854,10 @@ class InteractionMatrix:
         
         return self
     
-    def probabilities (self):
+    def probabilities (
+        self,
+        check_df_vars: bool = True
+        ):
 
         interaction_matrix_df = self.interaction_matrix_df
 
@@ -858,6 +865,12 @@ class InteractionMatrix:
             self.utility()
             interaction_matrix_df = self.interaction_matrix_df
 
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["U_ij"]
+                )
+
         utility_i = pd.DataFrame(interaction_matrix_df.groupby("i")["U_ij"].sum())
         utility_i = utility_i.rename(columns = {"U_ij": "U_i"})
 
@@ -876,7 +889,10 @@ class InteractionMatrix:
 
         return self
         
-    def flows (self):
+    def flows (
+        self,
+        check_df_vars: bool = True
+        ):
 
         interaction_matrix_df = self.interaction_matrix_df
 
@@ -885,10 +901,11 @@ class InteractionMatrix:
         if interaction_matrix_df["C_i"].isna().all():
             raise ValueError ("Error in flows calculation: Market size column in customer origins not defined. Use CustomerOrigins.define_marketsize()")
 
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["C_i"]
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["C_i"]
+                )
 
         if interaction_matrix_df["p_ij"].isna().all():
             self.probabilities()
@@ -900,15 +917,19 @@ class InteractionMatrix:
 
         return self
 
-    def marketareas(self):
+    def marketareas(
+        self,
+        check_df_vars: bool = True
+        ):
 
         interaction_matrix_df = self.interaction_matrix_df
         
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["E_ij"],
-            check_zero = False
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["E_ij"],
+                check_zero = False
+                )
         
         market_areas_df = pd.DataFrame(interaction_matrix_df.groupby("j")["E_ij"].sum())
         market_areas_df = market_areas_df.reset_index(drop=False)
@@ -947,9 +968,9 @@ class InteractionMatrix:
             if tc_weighting["func"] == "power":
                 interaction_matrix_df["t_ij_weighted"] = interaction_matrix_df["t_ij"] ** tc_weighting["param"]
             elif tc_weighting["func"] == "exponential":
-                interaction_matrix_df["t_ij_weighted"] = np.exp(tc_weighting["param"] * interaction_matrix_df['t_ij'])
+                interaction_matrix_df["t_ij_weighted"] = np.exp(tc_weighting["param"] * interaction_matrix_df["t_ij"])
             elif tc_weighting["func"] == "logistic":
-                interaction_matrix_df["t_ij_weighted"] = 1+np.exp(tc_weighting["param"][0] + tc_weighting["param"][1] * interaction_matrix_df['t_ij'])
+                interaction_matrix_df["t_ij_weighted"] = 1+np.exp(tc_weighting["param"][0] + tc_weighting["param"][1] * interaction_matrix_df["t_ij"])
             else:
                 raise ValueError ("Error in hansen accessibility calculation: Transport costs weighting is not defined.")
                         
@@ -1073,7 +1094,8 @@ class InteractionMatrix:
     def loglik(
         self,
         params,
-        fit_by = "probabilities"
+        fit_by = "probabilities",
+        check_df_vars: bool = True
         ):
         
         if fit_by not in ["probabilities", "flows"]:
@@ -1160,9 +1182,9 @@ class InteractionMatrix:
         
         interaction_matrix_copy = copy.deepcopy(self)
 
-        interaction_matrix_copy.utility()
-        interaction_matrix_copy.probabilities()
-        interaction_matrix_copy.flows()
+        interaction_matrix_copy.utility(check_df_vars = check_df_vars)
+        interaction_matrix_copy.probabilities(check_df_vars = check_df_vars)
+        interaction_matrix_copy.flows(check_df_vars = check_df_vars)
 
         interaction_matrix_df_copy = interaction_matrix_copy.get_interaction_matrix_df()
         
@@ -1196,7 +1218,8 @@ class InteractionMatrix:
         bounds: list = [(0.5, 1), (-3, -1)],        
         constraints: list = [],
         fit_by = "probabilities",
-        update_estimates: bool = True
+        update_estimates: bool = True,
+        check_df_vars: bool = True
         ):       
 
         supply_locations = self.supply_locations
@@ -1232,7 +1255,7 @@ class InteractionMatrix:
         ml_result = minimize(
             self.loglik,
             initial_params,
-            args=fit_by,
+            args=(fit_by, check_df_vars),
             method = method,
             bounds = bounds,
             constraints = constraints,
@@ -1745,7 +1768,8 @@ class HuffModel:
 
     def loglik(
         self,
-        params
+        params,
+        check_df_vars: bool = True
         ):
 
         if not isinstance(params, list):
@@ -1827,9 +1851,9 @@ class HuffModel:
 
         interaction_matrix_copy = copy.deepcopy(huff_model_copy.interaction_matrix)
 
-        interaction_matrix_copy = interaction_matrix_copy.utility()
-        interaction_matrix_copy = interaction_matrix_copy.probabilities()
-        interaction_matrix_copy = interaction_matrix_copy.flows()
+        interaction_matrix_copy = interaction_matrix_copy.utility(check_df_vars = check_df_vars)
+        interaction_matrix_copy = interaction_matrix_copy.probabilities(check_df_vars = check_df_vars)
+        interaction_matrix_copy = interaction_matrix_copy.flows(check_df_vars = check_df_vars)
 
         huff_model_copy = interaction_matrix_copy.marketareas()
 
@@ -1855,7 +1879,8 @@ class HuffModel:
         constraints: list = [],
         fit_by = "probabilities",
         update_estimates: bool = True,
-        check_numbers: bool = True
+        check_numbers: bool = True,
+        check_df_vars: bool = True
         ):
         
         if fit_by in ["probabilities", "flows"]:
@@ -1866,7 +1891,8 @@ class HuffModel:
                 bounds = bounds,        
                 constraints = constraints,
                 fit_by = fit_by,
-                update_estimates = update_estimates
+                update_estimates = update_estimates,
+                check_df_vars = check_df_vars
                 )
         
         elif fit_by == "totals":
@@ -1914,6 +1940,7 @@ class HuffModel:
             ml_result = minimize(
                 self.loglik,
                 initial_params,
+                args=check_df_vars,
                 method = method,
                 bounds = bounds,
                 constraints = constraints,
@@ -1999,11 +2026,11 @@ class HuffModel:
 
                 else:
                                        
-                    self.interaction_matrix.utility()
-                    self.interaction_matrix.probabilities()
-                    self.interaction_matrix.flows()                    
+                    self.interaction_matrix.utility(check_df_vars = check_df_vars)
+                    self.interaction_matrix.probabilities(check_df_vars = check_df_vars)
+                    self.interaction_matrix.flows(check_df_vars = check_df_vars) 
                     
-                    huff_model_new_marketareas = self.interaction_matrix.marketareas()
+                    huff_model_new_marketareas = self.interaction_matrix.marketareas(check_df_vars = check_df_vars)
                     self.market_areas_df["T_j"] = huff_model_new_marketareas.get_market_areas_df()["T_j"]
             
             self.interaction_matrix.metadata["fit"] = {
@@ -2372,7 +2399,8 @@ class MCIModel:
                   
     def utility(
         self,
-        transformation = "LCT"
+        transformation = "LCT",
+        check_df_vars: bool = True
         ):
         
         interaction_matrix = self.interaction_matrix        
@@ -2384,10 +2412,11 @@ class MCIModel:
         if interaction_matrix_df["A_j"].isna().all():
             raise ValueError ("Error in utility calculation: Attraction variable is not defined")
 
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["A_j", "t_ij"]
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["A_j", "t_ij"]
+                )
 
         customer_origins = interaction_matrix.get_customer_origins()
         customer_origins_metadata = customer_origins.get_metadata()
@@ -2476,7 +2505,8 @@ class MCIModel:
         
     def flows (
         self,
-        transformation = "LCT"
+        transformation = "LCT",
+        check_df_vars: bool = True
         ):
 
         interaction_matrix = self.interaction_matrix
@@ -2488,10 +2518,11 @@ class MCIModel:
         if interaction_matrix_df["C_i"].isna().all():
             raise ValueError ("Error in flows calculation: Market size column in customer origins not defined. Use CustomerOrigins.define_marketsize()")
 
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["C_i"]
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["C_i"]
+                )
 
         if interaction_matrix_df["p_ij"].isna().all():
             self.probabilities(transformation=transformation)
@@ -2504,15 +2535,19 @@ class MCIModel:
 
         return self
 
-    def marketareas (self):
+    def marketareas (
+        self,
+        check_df_vars: bool = True
+        ):
 
         interaction_matrix = self.interaction_matrix
         interaction_matrix_df = interaction_matrix.get_interaction_matrix_df()
         
-        check_vars(
-            df = interaction_matrix_df,
-            cols = ["E_ij"]
-            )
+        if check_df_vars:
+            check_vars(
+                df = interaction_matrix_df,
+                cols = ["E_ij"]
+                )
         
         market_areas_df = pd.DataFrame(interaction_matrix_df.groupby("j")["E_ij"].sum())
         market_areas_df = market_areas_df.reset_index(drop=False)
@@ -2735,7 +2770,7 @@ def load_interaction_matrix(
     xlsx_sheet: str = None,
     crs_input = "EPSG:4326",
     crs_output = "EPSG:4326",
-    check_df_vars = True
+    check_df_vars: bool = True
     ):    
 
     if isinstance(data, pd.DataFrame):
@@ -2976,7 +3011,7 @@ def load_marketareas(
     csv_decimal = ",", 
     csv_encoding="unicode_escape",
     xlsx_sheet: str = None,
-    check_df_vars = True
+    check_df_vars: bool = True
     ):    
 
     if isinstance(data, pd.DataFrame):

{huff-1.5.3 → huff-1.5.4}/huff/tests/tests_huff.py

@@ -4,8 +4,8 @@
 # Author:      Thomas Wieland 
 #              ORCID: 0000-0001-5168-9846
 #              mail: geowieland@googlemail.com              
-# Version:     1.5.3
-# Last update: 2025-07-15 17:22
+# Version:     1.5.4
+# Last update: 2025-07-18 18:06
 # Copyright (c) 2025 Thomas Wieland
 #-----------------------------------------------------------------------
 

{huff-1.5.3 → huff-1.5.4}/huff.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: huff
-Version: 1.5.3
+Version: 1.5.4
 Summary: huff: Huff Model Market Area Analysis
 Author: Thomas Wieland
 Author-email: geowieland@googlemail.com
@@ -18,12 +18,10 @@ Thomas Wieland [ORCID](https://orcid.org/0000-0001-5168-9846) [EMail](mailto:geo
 See the /tests directory for usage examples of most of the included functions.
 
 
-## Updates v1.5.3
+## Updates v1.5.4
 - Bugfixes:
-  - InteractionMatrix.summary() and HuffModel.summary(): No KeyError when param is None anymore
-- Extensions:
-  - Confidence intervals (bootstrap) for Huff Model ML estimations
-  - ValueError and KeyError are more precise (function is included)
+  - Use of check_vars() is now optional (default: True)
+  - Correction of args argument when calling minimize() in huff_ml_fit() 
   
 
 ## Features

{huff-1.5.3 → huff-1.5.4}/setup.py

@@ -7,7 +7,7 @@ def read_README():
     
 setup(
     name='huff',
-    version='1.5.3',
+    version='1.5.4',
     description='huff: Huff Model Market Area Analysis',
     packages=find_packages(include=["huff", "huff.tests"]),
     include_package_data=True,