Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error when running MIDAS on KEGG metabolism gene sets #1

Open
dBenedek opened this issue Nov 23, 2021 · 0 comments
Open

Error when running MIDAS on KEGG metabolism gene sets #1

dBenedek opened this issue Nov 23, 2021 · 0 comments

Comments

@dBenedek
Copy link

Hello,

MIDAS runs perfectly when running on the example gene sets, however, when running on KEGG metabolism pathways, I get the following message:

[Info] Calculate edge activty & determine differentially activated subpaths
[Info] Target Pathway:  hsa04110/Cell_cycle
[Info] Done.
[Info] Target Pathway:  hsa04010/MAPK_sig
[Info] Done.
[Info] Target Pathway:  hsa04151/PI3K_AKT
[Info] Done.
[Info] Target Pathway:  hsa00010/Glycolysis_gluconeog
Traceback (most recent call last):
  File "..///bin//calc_edge_activity.py", line 174, in <module>
    (total_SAS_mean_results, edge_activity_dic) = calc_edge_activity_whole(flow_graph_dic, centrality_dic, edge_info_dic, act_iht_filter_dic,max_expr_dic, edge_weight_case, edge_weight_dic, act_iht_case, sample_num, p_cores)
  File "/san1/BD/softwares/MIDAS/bin/calc_edge_act_parallel.py", line 148, in calc_edge_activity_whole
    total_SAS_mean_results = map(lambda x : (x[0]/total_num, x[1]/total_num), total_SAS_sum_results)
  File "/san1/BD/softwares/MIDAS/bin/calc_edge_act_parallel.py", line 148, in <lambda>
    total_SAS_mean_results = map(lambda x : (x[0]/total_num, x[1]/total_num), total_SAS_sum_results)
ZeroDivisionError: float division by zero
grep: ./test_result//single_pathway_analysis///hsa00010/entry//Edge_act.txt.act_iht_result.txt.MGD: Nincs ilyen fájl vagy könyvtár
Traceback (most recent call last):
  File "..///bin//kruskal_subpath_mining_ver_pac_1.py", line 64, in <module>
    header = next(edge_act_file_reader)
StopIteration
Traceback (most recent call last):
  File "..///bin//kruskal_perm_subpath_parallel.py", line 71, in <module>
    header = next(edge_act_file_reader)
StopIteration
Error in read.table(args[1], sep = "\t", header = T, stringsAsFactors = F) : 
  no lines available in input
Execution halted
[Info] Done.

Pathway set file:

hsa04110	Cell_cycle
hsa04010	MAPK_sig
hsa04151	PI3K_AKT
hsa00010	Glycolysis_gluconeog

Trying with other KEGG metabolism pathways gives the same error messages.
I am attaching the full log file (nohup).

I would appreciate your help,
Benedek

nohup.txt

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant