Skip to content

FUNCE throws error when trying to run more than one terms file #16

@JohnHadish

Description

@JohnHadish

Term List 1 Runs fine:

FUNC-E \
  --background "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_background.txt" \
  --query_list "/home/jah/Documents/05-galaapple_postharvesttreatment/04-Differential_Expression/61-tsv/91-maSigPro_CA.tsv" \
  --terms "/home/jah/Documents/00-reference/FUNC-E/terms/IPR.terms.txt" \
  --terms2features "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_genes2IPR.txt" \
  --ecut 0.01 \
  --outprefix "61-tsv/maSigPro_CA"

Term list 2, runs fine:

FUNC-E \
  --background "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_background.txt" \
  --query_list "/home/jah/Documents/05-galaapple_postharvesttreatment/04-Differential_Expression/61-tsv/91-maSigPro_CA.tsv" \
  --terms "/home/jah/Documents/00-reference/FUNC-E/terms/KEGG.terms.txt" \
  --terms2features "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_genes2KEGG.txt" \
  --ecut 0.01 \
  --outprefix "61-tsv/maSigPro_CA"

Combine the two as described in the documentation, does not run:

FUNC-E \
  --background "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_background.txt" \
  --query_list "/home/jah/Documents/05-galaapple_postharvesttreatment/04-Differential_Expression/61-tsv/91-maSigPro_CA.tsv" \
  --terms "/home/jah/Documents/00-reference/FUNC-E/terms/KEGG.terms.txt" \
  "/home/jah/Documents/00-reference/FUNC-E/terms/IPR.terms.txt" \
  --terms2features "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_genes2KEGG.txt" \
  "/home/jah/Documents/00-reference/07-Malus_x_domestica-genome_GDDH13_v1.1/Malus_x_domestica_GDDH13_v1.1_genes2IPR.txt" \
  --ecut 0.01 \
  --outprefix "61-tsv/maSigPro_CA"

Error message for the combined:

AttributeError: 'float' object has no attribute 'log'

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
  File "/home/jah/.local/bin/FUNC-E", line 8, in <module>
    sys.exit(func_e())
  File "/home/jah/.local/lib/python3.10/site-packages/func_e/cmd.py", line 100, in func_e
    fe.run(modules=args.module, vocabs=args.vocab)
  File "/home/jah/.local/lib/python3.10/site-packages/func_e/FUNC_E.py", line 241, in run
    self.doClustering(modules)
  File "/home/jah/.local/lib/python3.10/site-packages/func_e/FUNC_E.py", line 605, in doClustering
    self.doModuleClustering(module)
  File "/home/jah/.local/lib/python3.10/site-packages/func_e/FUNC_E.py", line 580, in doModuleClustering
    stats = self._calculateClusterStats(final_clusters, module)
  File "/home/jah/.local/lib/python3.10/site-packages/func_e/FUNC_E.py", line 528, in _calculateClusterStats
    gmean = stats.gmean(eterms['Fishers_pvalue'])
  File "/home/jah/.local/lib/python3.10/site-packages/scipy/stats/_axis_nan_policy.py", line 503, in axis_nan_policy_wrapper
    res = hypotest_fun_out(*samples, **kwds)
  File "/home/jah/.local/lib/python3.10/site-packages/scipy/stats/_stats_py.py", line 307, in gmean
    log_a = np.log(a)
TypeError: loop of ufunc does not support argument 0 of type float which has no callable log method

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions