1

我想知道是否可以同时在 MrBayes 中修复它们?

我看到一些关于修复树拓扑的先前帖子是首先定义树然后使用“固定”功能(假设我们已经定义了species_topology)

第一行固定树形拓扑,第二行定义移动到其他拓扑的概率为零。

Prset topologypr=fixed(species_topology);
propset eTBR(Tau)$prob=0;

我想知道我也可以做类似的事情来修复分支长度吗?就转移到分支长度的建议而言,我认为我们应该使用

propset nslider(V)$prob=0

但我不确定要输入什么:

Prset brlenpr= 

提前致谢!我在 MrBayes 的邮件列表中搜索了很多帖子,但他们没有关于同时修复它们的问题。通常,他们只是修复拓扑。

4

1 回答 1

0

在搜索了很多 MrBayes 教程之后,我找到了一种自己实现这一点的方法。在 MrBayes 3.2 中,它允许我们固定树拓扑和分支长度。我们可以将拓扑和分支长度的先验设置为 newick 格式的固定树,然后将从初始状态移动的概率定义为零。请记住,您必须在保存序列信息数据的链接文件中定义具有固定分支长度的固定树。

例如,你可以这样定义你的树:

begin trees;
tree tree_1 = ((t6:0.3279207193,t9:0.9545036491):0.04205953353,t5:0.8895393161,(((t4:0.6557057991,t10:0.7085304682):0.9942697766,t3:0.5440660247):0.6405068138,((t8:0.1471136473,(t2:0.9022990451,t1:0.6907052784):0.9630242325):0.2891597373,t7:0.7954674177):0.5941420204):0.9388911405);
;
end;

然后在您的 MrBayes 模板中,您可以使用以下命令来修复树拓扑和分支长度:

prset topologypr=fixed(tree_1);
prset brlenspr=fixed(tree_1);

Tau 表示拓扑,V 表示分支长度。

于 2015-06-02T18:39:56.533 回答