Hopefully, this works

Author: jonbrenas

malariagen_data/anoph/fst.py

@@ -539,8 +539,11 @@ def plot_pairwise_average_fst(
             if annotation == "standard error":
                 fig_df.loc[cohort1, cohort2] = se
             elif annotation == "Z score":
-                zs = fst / se
-                fig_df.loc[cohort1, cohort2] = zs
+                try:
+                    zs = fst / se
+                    fig_df.loc[cohort1, cohort2] = zs
+                except ZeroDivisionError:
+                    fig_df.loc[cohort1, cohort2] = np.NaN
             else:
                 fig_df.loc[cohort1, cohort2] = fst