diff --git a/nolds/measures.py b/nolds/measures.py
index <HASH>..<HASH> 100644
--- a/nolds/measures.py
+++ b/nolds/measures.py
@@ -1382,7 +1382,7 @@ def hurst_multifractal_dm(data, qvals=[1], max_dists=range(5, 20)):
     ])
   hhcorr = np.array(hhcorr, dtype="float32")
   H = np.array([
-    _aste_line_fit(np.log(range(1, md)), np.log(hhcorr[:md-1, qi]))[1]
+    _aste_line_fit(np.log(range(1, md+1)), np.log(hhcorr[:md, qi]))[1]
     for qi in range(len(qvals))
     for md in max_dists
   ], dtype="float32").reshape(len(qvals), len(max_dists))
bugfix: adjusts limits for line fit (start at 1, include upper limit)