cmks / DAS_Tool

DAS Tool
Other
140 stars 17 forks source link

Error in aggregate.data.frame #63

Closed nick-youngblut closed 3 years ago

nick-youngblut commented 3 years ago

DAS-tool 1.1.2 generated the following error:

Closing the input file...  [0s]
Closing the output file...  [0s]
Closing the database file...  [0s]
Deallocating taxonomy...  [0s]
Total time = 0.881s
Reported 1 pairwise alignments, 1 HSPs.
1 queries aligned.
The host system is detected to have 540 GB of RAM. It is recommended to increase the block size for better performance using these parameters : -b12 -c1
starting annotations of single copy cogs...
successfully finished
calculating contig lengths.
WARNING: non-unique bin-IDs found
    creating unique bin-IDs

Error in aggregate.data.frame(arc_scg["count"], by = arc_scg[c("Archaeal.SCG",  :
  no rows to aggregate
Calls: cherry_pick -> aggregate -> aggregate.data.frame
Execution halted

The top of the log:

# Running das-tool
Running DAS Tool using 8 threads.
predicting genes using Prodigal V2.6.3: February, 2016
identifying single copy genes using diamond version 2.0.5
database name of all proteins is /ebio/abt3_projects/Georg_animal_feces/bin/llmga/.snakemake/conda/8f4a60b0/share/das_tool-1.1.2-1/db/bac.all.faa
database name of SCGs is /ebio/abt3_projects/Georg_animal_feces/bin/llmga/.snakemake/conda/8f4a60b0/share/das_tool-1.1.2-1/db/bac.scg.faa
database lookup is /ebio/abt3_projects/Georg_animal_feces/bin/llmga/.snakemake/conda/8f4a60b0/share/das_tool-1.1.2-1/db/bac.scg.lookup
diamond v2.0.5.143 (C) Max Planck Society for the Advancement of Science
Documentation, support and updates available at http://www.diamondsearch.org

#CPU threads: 64
Scoring parameters: (Matrix=BLOSUM62 Lambda=0.267 K=0.041 Penalties=11/1)
Database input file: /ebio/abt3_projects/Georg_animal_feces/bin/llmga/.snakemake/conda/8f4a60b0/share/das_tool-1.1.2-1/db/bac.all.faa
Opening the database file...  [0.001s]

It appears that there's no check on the data.frame size prior to the aggregate function

GeoMicroSoares commented 3 years ago

Same thing happening for me - any news on the issue?

cmks commented 3 years ago

This issue should be solved now.