Fix enrichment stat p_pos from >= to > (#257)

angelolab · Oct 2, 2020 · eaf429e · eaf429e
1 parent 0f7a0aa
commit eaf429e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/ark/utils/spatial_analysis_utils.py b/ark/utils/spatial_analysis_utils.py
@@ -307,7 +307,7 @@ def calculate_enrichment_stats(close_num, close_num_rand):
             # Calculate z score based on distribution
             z[j, k] = (close_num[j, k] - muhat[j, k]) / sigmahat[j, k]
             # Calculate both positive and negative enrichment p values
-            p_pos[j, k] = (1 + (np.sum(tmp >= close_num[j, k]))) / (bootstrap_num + 1)
+            p_pos[j, k] = (1 + (np.sum(tmp > close_num[j, k]))) / (bootstrap_num + 1)
             p_neg[j, k] = (1 + (np.sum(tmp < close_num[j, k]))) / (bootstrap_num + 1)
 
     # Get fdh_br adjusted p values