转换为SAtree
SA树
对象可以从化石
对象和关联的树。提示(已灭绝和现存)将是在标签上标明其采样的物种,然后附上索引:给定物种的最古老样本将获得索引1和所有其他样品将按从大到小的顺序订购。这个SA树
格式还包括一个字段完成
这表明树是否应被视为包含处理或仅采样血统。请注意,如果树是完整的给定灭绝物种的最年轻尖端代表灭绝事件而不是化石样本。
t=猿::树(6)f=模拟化石泊松(比率=2,树=t)SAt=SAtree.from.complass(树=t,化石=f)打印(SAt$tree)
## ##具有12个尖端和11个内部节点的系统发育树。## ##提示标签:##t5_1、t6_1、t4_1、t3_1、t1_1、t2_2。。。## ##生根;包括分支长度。
打印(SAt$化石)
##sp边缘hmin hmax h tip.label##1 6 0.5815387 0.5815387 0.5815387 t2_1##3 8 8 1.6669608 1.6669968 1.666968 t8_1##2 8 8 1.5907748 1.59077481.59077480 t8_2##4 10 10 1.3482607 1.3482607,1.3482607t10_1##5 10 10 1.2132023 1.2132022 3 1.213202 3 t10_2##6 10 10 1.2125041 1.212504l 1.2125041t10_3##化石记录,6次出现,代表3个物种##未使用分类法模拟化石记录:假设所有物种形成事件都是对称的
打印(SAt$tree$complete)
##[1]正确
其他有用功能
其他功能位于化石模拟人
修改树。修枝化石
将清除所有中间化石并只保留任何物种的第一次和最后一次出现。
SAt_pruned=修剪化石(tree=SAt$tree)绘图(SAt_pruned)
样本.树.来自组合
删除所有未采样树的谱系,可以选择应用采样概率现有样本,或只保留现有样本的指定列表。
SAt_sampled=sampled.tree.from.combined(树=SAt$tree)绘图(SAt_sampled)
这两个函数可以应用于任何类型的对象门(phylo)
,但设计用于树木,包括化石样本。