lliu1871 / phybase

Phybase provides functions to read, write, manipulate, simulate, estimate, and summarize phylogenetic trees
GNU General Public License v3.0
4 stars 1 forks source link

what does the “#0.1” in example trees mean? #3

Closed Dawei-Xia closed 2 years ago

Dawei-Xia commented 2 years ago

Hi,professor liu: The function “sim.coaltree.sp.mu” generates a random gene tree from the species tree under the non-clock species tree model. I'm confused “#0.1” in the example tree “ (((A:0.5,B:0.5):1#0.1,C:1.5):1#0.1,D:2.5)#0.1;)”,and when I get a certain theta value (μT/τ ), how to set it in the species tree used for simulation analysis.

Thanks for all your help! Best wishes

lliu1871 commented 2 years ago

In tree (((A:0.5,B:0.5):1#0.1,C:1.5):1#0.2,D:2.5)#0.01, #0.01 indicates that theta = 0.01 for the root and #0.2 indicates that theta = 0.2 for the internal branch leading to the group ((A,B),C)

Liang

Liang Liu Graduate Coordinator Department of Statistics Institute of Bioinformatics University of Georgia Phylogenetics Lab | Phylogenetics Lab (uga.edu)https://faculty.franklin.uga.edu/lliu/

[cid:AF1C81C8-35FF-48FE-9242-011CC2CC4CA0]

On Aug 30, 2022, at 9:19 AM, Dawei Xia @.**@.>> wrote:

[EXTERNAL SENDER - PROCEED CAUTIOUSLY]

Hi,professor liu: The function “sim.coaltree.sp.mu” generates a random gene tree from the species tree under the non-clock species tree model. I'm confused “#0.1” in the example tree “ (((A:0.5,B:0.5):1#0.1,C:1.5):1#0.1,D:2.5)#0.1;)”,and when I get a certain theta value (μT/τ ), how to set it in the species tree used for simulation analysis.

Thanks for all your help! Best wishes

— Reply to this email directly, view it on GitHubhttps://github.com/lliu1871/phybase/issues/3, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ACZITQWPGC73TBR5TJHTVJDV3YC7PANCNFSM6AAAAAAQAMZMPM. You are receiving this because you are subscribed to this thread.Message ID: @.***>

Dawei-Xia commented 2 years ago

In tree (((A:0.5,B:0.5):1#0.1,C:1.5):1#0.2,D:2.5)#0.01, #0.01 indicates that theta = 0.01 for the root and #0.2 indicates that theta = 0.2 for the internal branch leading to the group ((A,B),C) Liang ------------------------------------- Liang Liu Graduate Coordinator Department of Statistics Institute of Bioinformatics University of Georgia Phylogenetics Lab | Phylogenetics Lab (uga.edu)https://faculty.franklin.uga.edu/lliu/ [cid:AF1C81C8-35FF-48FE-9242-011CC2CC4CA0] On Aug 30, 2022, at 9:19 AM, Dawei Xia @.**@.>> wrote: [EXTERNAL SENDER - PROCEED CAUTIOUSLY] Hi,professor liu: The function “sim.coaltree.sp.mu” generates a random gene tree from the species tree under the non-clock species tree model. I'm confused “#0.1” in the example tree “ (((A:0.5,B:0.5):1#0.1,C:1.5):1#0.1,D:2.5)#0.1;)”,and when I get a certain theta value (μT/τ ), how to set it in the species tree used for simulation analysis. Thanks for all your help! Best wishes — Reply to this email directly, view it on GitHub<#3>, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ACZITQWPGC73TBR5TJHTVJDV3YC7PANCNFSM6AAAAAAQAMZMPM. You are receiving this because you are subscribed to this thread.Message ID: @.***>

I was hesitant but not sure. thank you very much for your prompt reply

best wishes