Menu

#3 No such file or directory at RaFAH.pl line 310

1.0
open
nobody
None
2022-10-31
2022-06-09
No

Hello, thank you for your tool. I installed the dependency you mentioned with CONDA, but the following problems occurred when running the test data.

Do you have any suggestions? Thank you very much for your help.

1 Attachments

Discussion

  • Felipe Hernandes Coutinho

    Hi there! Which dependencies are you referring to specifically? There is no conda enviroment for RaFAH. Also, which version of RaFAH are you running and with which command?

     
  • Julia Gustavsen

    Julia Gustavsen - 2022-09-30

    Hello, I think the comment refers to the Docker image as I also get the same error when running the Docker image "--predict --merged_cds_file_name". Do you have any suggestions? Thanks, Julia

     
  • Felipe Hernandes Coutinho

    Hi Julia, can you please specify the full command line you are using?

     
  • Julia Gustavsen

    Julia Gustavsen - 2022-10-01

    Thanks for the response!
    docker run --mount type=bind,source=/Users/juliagustavsen/Documents/viral_analysis/results,target=/results/ fhcoutinho/rafah perl /usr/bin/RaFAH.pl --predict --merged_cds_file_name /results/Total_CP.fasta --output_dir /results

     
  • Felipe Hernandes Coutinho

    Is "/results/Total_CP.fasta" empty by any chance? If there is something in it, are they protein sequences in fasta format?

     
  • Julia Gustavsen

    Julia Gustavsen - 2022-10-03

    yes they are protein sequences in fasta format.

    We can also see an issue (as in the original comment) using the test data - TS3_Toy_Genomes.fasta (downloaded from here https://sourceforge.net/projects/rafah/files/Data/Toy_Set.tgz/download). If I download the data to my "docker" folder and run the command as described in the readme I get this:

    docker run --mount type=bind,source=/Users/juliagustavsen/Documents/thesis_projects/Kuwait_gp20_analysis/docker,target=/results/ fhcoutinho/rafah perl /usr/bin/RaFAH.pl --predict --genomes_dir /results/ --extension .fasta --file_prefix RaFAH_1

    Running host prediction mode
    Warning! ./ already exists. Will overwrite previous files.
    Indexing sequences from /results/
    Processing TS3_Toy_Genomes.fasta
    Processed 10 Genomic Sequences
    Running Prodigal
    Indexing sequences from .//RaFAH_1_CDS_Prediction.faa
    Running hmmsearch. Query: .//RaFAH_1_CDS_Prediction.faa DB: /home/RaFAH/Data/HP_Ranger_Model_3_Filtered_0.9_Valids.hmm
    Obtained 43644 ids from /home/RaFAH/Data/HP_Ranger_Model_3_Valid_Cols.txt
    Parsing .//RaFAH_1_CDSxClusters_Prediction
    Detected 705 OGs across 10 genomic sequences
    Performing host prediction
    [1] "Loading Model from /home/RaFAH/Data/MMSeqs_Clusters_Ranger_Model_1+2+3_Clean.RData"
    No such file or directory at /usr/bin/RaFAH.pl line 311.
    Parsing output of host prediction .//RaFAH_1_Host_Predictions.tsv

     
  • Felipe Hernandes Coutinho

    Ok. And there is no "RaFAH_1_Host_Predictions.tsv" file generated, correct? It seems R is dying after loading the models into memory. This could happen if there is not enough RAM available. RaFAH needs at least 80 Mb of RAM to run properly, so it will most likely not work on a personal computer. Can you try running it on your server and let me know if you still get the same error?

     
  • Julia Gustavsen

    Julia Gustavsen - 2022-10-03

    Thanks for the message. Tracking the memory usage of the Docker container appears to go up to 8Gigs and then is closer to 3 Gigs memory use when it has the file not found error.

     
  • Felipe Hernandes Coutinho

    My bad, I meant 80GB (considering a dataset with thousands of viral sequences and using the predict mode). It only uses this much for a very short time, but enough to get killed if there is not enough much memory available.

     
  • Julia Gustavsen

    Julia Gustavsen - 2022-10-07

    Ok thanks. I ran it on AWS and it was fine. I think it would be helpful for users to know min requirements (or to improve the error handling). I cannot make a PR here on Sourceforge or else I would have submitted some propositions. Thanks for your help, Julia

     
  • Felipe Hernandes Coutinho

    Thanks Julia, I will add min requirements to the readme file. Also, I just created a RaFAH repo on github: https://github.com/felipehcoutinho/RaFAH Its about time I center everything there anyway. Feel free to submit PRs through there. Best, F

     

Log in to post a comment.