diff --git a/finemapper.py b/finemapper.py
index 30af1a1..9d450b9 100644
--- a/finemapper.py
+++ b/finemapper.py
@@ -653,7 +653,7 @@ def estimate_h2_hess(self, prop_keep=0.005, R_cutoff=0.99, pvalue_bound=None):
         #take a maximally independent subset
         np.fill_diagonal(R_pot_csnp,0)
         import networkx as nx
-        G = nx.from_numpy_matrix(np.abs(R_pot_csnp)>R_cutoff)
+        G = nx.from_numpy_array(np.abs(R_pot_csnp)>R_cutoff)
         np.fill_diagonal(R_pot_csnp,1)
         inds = np.sort(nx.maximal_independent_set(G))