Skip to content

Commit

Permalink
update compare method to call WeSME utility script
Browse files Browse the repository at this point in the history
  • Loading branch information
ashuaibi7 committed Jan 6, 2025
1 parent 6f519aa commit a2dcdf8
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion src/dialect/utils/compare.py
Original file line number Diff line number Diff line change
Expand Up @@ -4,6 +4,7 @@
from dialect.utils.fishers import run_fishers_exact_analysis
from dialect.utils.discover import run_discover_analysis
from dialect.utils.megsa import run_megsa_analysis
from dialect.utils.wesme import run_wesme_analysis


# ---------------------------------------------------------------------------- #
Expand Down Expand Up @@ -64,11 +65,14 @@ def run_comparison_methods(cnt_mtx, bmr_pmfs, out, k):
logging.info("Running MEGSA...")
megsa_df = run_megsa_analysis(cnt_df, interactions)

# TODO: Implement WeSME/WeSCO Analysis
logging.info("Running WeSME/WeSCO...")
wesme_df = run_wesme_analysis(cnt_df, out, interactions)

# TODO: Implement SELECT Analysis

merged_df = pd.merge(fisher_df, discover_df, on=["Gene A", "Gene B"], how="inner")
merged_df = pd.merge(merged_df, megsa_df, on=["Gene A", "Gene B"], how="inner")
merged_df = pd.merge(merged_df, wesme_df, on=["Gene A", "Gene B"], how="inner")
comparison_interaction_fout = f"{out}/comparison_interaction_results.csv"
merged_df.to_csv(comparison_interaction_fout, index=False)
logging.info(f"Comparison results saved to {comparison_interaction_fout}")
Expand Down

0 comments on commit a2dcdf8

Please sign in to comment.