Correct FDR calculation

FredHutch · Mar 13, 2024 · ff6fe04 · ff6fe04
1 parent bf6bd0e
commit ff6fe04
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/bin/plot_metagenomes.py b/bin/plot_metagenomes.py
@@ -66,7 +66,7 @@ def __init__(
             kw: self.adata.var[kw].fillna(val)
             for kw, val in [("estimate", 0), ("std_error", 0), ("p_value", 0.99)]
         })
-        self.adata.var["qvalue"] = multipletests(self.adata.var["p_value"].fillna(0), 0.1, "fdr_bh")[1]
+        self.adata.var["qvalue"] = multipletests(self.adata.var["p_value"].fillna(1), 0.1, "fdr_bh")[1]
         self.adata.var["neg_log10_pvalue"] = -np.log10(self.adata.var["p_value"])
         self.adata.var["neg_log10_qvalue"] = -np.log10(self.adata.var["qvalue"])
         for line in str(self.adata).split("\n"):