diff --git a/netin/models/dh_model.py b/netin/models/dh_model.py
index 8a4425e..83b2c00 100644
--- a/netin/models/dh_model.py
+++ b/netin/models/dh_model.py
@@ -54,7 +54,7 @@ def _initialize_lfms(self):
         super()._initialize_lfms()
         self.h = TwoClassHomophily.from_two_class_homophily(
             node_class_values=self.graph.get_node_class(CLASS_ATTRIBUTE),
-            homophily=(self.h_m, self.h_M))
+            homophily=(self.h_M, self.h_m))
 
     def compute_target_probabilities(self, source: int)\
         -> np.ndarray:
diff --git a/netin/models/dpah_model.py b/netin/models/dpah_model.py
index 1ffd661..15b3611 100644
--- a/netin/models/dpah_model.py
+++ b/netin/models/dpah_model.py
@@ -58,7 +58,7 @@ def _initialize_lfms(self):
         super()._initialize_lfms()
         self.h = TwoClassHomophily.from_two_class_homophily(
             node_class_values=self.graph.get_node_class(CLASS_ATTRIBUTE),
-            homophily=(self.h_m, self.h_M))
+            homophily=(self.h_M, self.h_m))
         self.pa = InDegreePreferentialAttachment(
             graph=self.graph, N=self._n_nodes_total)
 
diff --git a/netin/models/homophily_model.py b/netin/models/homophily_model.py
index 2fec538..1fd3464 100644
--- a/netin/models/homophily_model.py
+++ b/netin/models/homophily_model.py
@@ -58,4 +58,4 @@ def compute_target_probabilities(self, source: int) -> np.ndarray:
     def _initialize_lfms(self):
         self.h = TwoClassHomophily.from_two_class_homophily(
             node_class_values=self.graph.get_node_class(CLASS_ATTRIBUTE),
-            homophily=(self.h_m, self.h_M))
+            homophily=(self.h_M, self.h_m))