安装程序包( " 古巴迪 " )
图书馆( 开发工具 ) 开发工具 :: 安装github( " brpetrucci/古巴迪 " )
设置种子( 1 ) 编号0 <- 1 # 物种初始数量
λ <- 0.1 # 物种形成率
亩 <- 0.05 # 消光率
t最大 <- 30 # 最大模拟时间
# 运行模拟
模拟 <- bd.sim(bd-sim)( 编号0 , λ , 亩 , t最大 )
设置种子( 1 ) ρ <- 1 # 采样率
箱子 <- 序列( t最大 , 0 , - 1 ) # 模拟地质层段的东西
# 获取化石出现时间的数据帧
化石 <- 样品.叶片( 模拟 = 模拟 , ρ = ρ , t最大 = t最大 , 箱子 = 箱子 ) # 可视化模拟与化石赋存 图纸.sim( 模拟 , 化石 )
物理(phy) <- make.phylo品牌( 模拟 ) # 用模拟组制作系统发育树
猿 :: phylo绘图( 物理(phy) , root.edge(根.边缘) = 真的 ) # 用树干绘制(需要APE)
猿 :: axisPhylo() # 添加轴