terminate called after throwing an instance of 'std::out_of_range'

78 views
Skip to first unread message

Shanshan Gu

unread,
Jan 4, 2023, 2:32:42 AM1/4/23
to hahnlab-cafe
Hi,
After running cafe5 two times, the lamuda value is converged. However, error always happened during or after Computing pvalues... step

Error:
Inferring processes for Gamma model
Score (-lnL): 262206.43934379
Maximum possible lambda for this topology: 0.0024154589371981
Computing pvalues...
terminate called recursively
terminate called after throwing an instance of 'std::out_of_range'
terminate called recursively
terminate called recursively
/opt/gridview/slurm/spool/slurmd/job849807/slurm_script: line 24: 321186 Aborted                 $cafe5 -c 20 -t $species_tree -i $final_family_count -k 3 -o $result_dir


This is my command:
Command line: cafe5 -c 20 -t xenTro10.multiz14ways.nwk -i Orthogroups.GeneCount.forCafe.delete20Orthogroups.txt -k 3 -o cafe_result 


By the way, tuning the count of core from 20 to 8 didn't work.


Best Regards,

Hahn, Matthew

unread,
Jan 4, 2023, 9:55:19 AM1/4/23
to Shanshan Gu, hahnlab-cafe
Hi,

At the moment, running CAFE with the gamma model does not support calculating p-values as well. But if you want to know the fastest families in your dataset, you can look in the rate-category with the highest lambda value.


Matt

--
You received this message because you are subscribed to the Google Groups "hahnlab-cafe" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hahnlabcafe...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hahnlabcafe/8f9ed967-5e9d-4a7e-b531-5d6175d12a9fn%40googlegroups.com.

Reply all
Reply to author
Forward
Message has been deleted
0 new messages