diff --git a/riskmapjnr/benchmark/defrate_per_class.py b/riskmapjnr/benchmark/defrate_per_class.py
index 3548e44..78d2bf1 100644
--- a/riskmapjnr/benchmark/defrate_per_class.py
+++ b/riskmapjnr/benchmark/defrate_per_class.py
@@ -133,7 +133,9 @@ class is created (see ``tab_file_defrate``).
     # Relative deforestation rate from model (not annual)
     if period in ["validation", "forecast"]:
         df_mod = pd.read_csv(deforate_model)
-        df["rate_mod"] = (df_mod["ndefor"] / df_mod["nfor"]).values
+        # Some defor class might be missing in the future
+        # so we need to join the tables to get rate_mod.
+        df = df.merge(right=df_mod, on="cat", how="left", suffixes=(None, "_mod"))
     else:
         df["rate_mod"] = df["ndefor"] / df["nfor"]