-
Notifications
You must be signed in to change notification settings - Fork 4
max() arg is an empty sequence #3
Description
Thanks for developing a very much needed tool for SynCom construction, I had problem running the mimic at the first step:
here are the info
(/vol/local/conda_envs/mimic2) [xuxm@blis JGI_com140]$ python /vol/local/xuxm/MiMiC2/MiMiC2.py
-g /vol/local/xuxm/JGI_com140/JGI_com140_bins_refiltered-profile.txt
-t /vol/local/xuxm/JGI_com140/com140_R0001/gtdbtk_out/gtdbtk.bac120.summary.tsv
--taxonomiclevel s
-s /vol/local/xuxm/JGI_com140/JGI_com140_sample_profile.tsv
-p /vol/local/xuxm/MiMiC2/datasets/core/Pfam-A.clans.tsv
--models /vol/local/xuxm/JGI_com140/bins/good_bins/RDS/
-c 4
-o 140SynCom
Reading in pfam vector files.
Number of pfams studied; 17929
Number of genomes studied; 6
Number of GEMs; 6
:: Reading in pfam and vector files
Total samples = 1
All samples will be studied as a single group, without comparative weighting.
Core proteins in Group 1 = 4660
::: Determined functions weight.
:::: Step 1: Iterative sample-specific selection of strains. ::::
0%| | 0/1 [00:01<?, ?it/s]
Traceback (most recent call last):
File "/vol/local/xuxm/MiMiC2/MiMiC2.py", line 367, in
max_match = max(genome_scores, key=genome_scores.get) # Identify the best consortia member this round
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
ValueError: max() arg is an empty sequence
I also wonder if I only have 6 MAGs in the sample, does it mean the SynCom would be the 6 members? No need to further running the MiMiC?
Thanks in advance for your help!