Can Phylo-HMRF do pairwise comparison?
Opened this issue · 3 comments
Hi,
Very delicacy approach Phylo-HMRF has! I would like to try it in my research.
I have a question though, can Phylo-HMRF do comparison between two species? and how can I do it?
cheers.
Hi @llq0325,
Thank you very much for your interest in Phylo-HMRF and thank you for the question. Yes, Phylo-HMRF can perform comparison between two species. In this case, the phylogenetic tree is simplified to be a tree with four nodes: one remote ancestor node (root node), one ancestor node, and two leaf nodes (representing the two species to study). Please edit the input files edge.1.txt, branch_length.1.txt, and species_name.1.txt accordingly. Please see the README about how to prepare the input files and the input data. Please let me know if you have any other questions. Thanks!
Best regards!
Yang Yang
Hi Yang,
Thank you for replying.
In current version of utility.py, you hard code the species name list. Maybe you can consider leave them out for a wider use.
https://github.com/ma-compbio/Phylo-HMRF/blob/b20f4e901d66f08de5a6d65b0d693583147f5236/utility.py#L2476
regards,
Langqing
Hi Langqing @llq0325 ,
Thank you very much for letting me know this problem. I have updated utility.py accordingly. Please let me know if you have any other questions. Many thanks!
Best regards!
Yang Yang