库(devtools) 安装github(“fossilsim/fossilsim”)
#使用猿模拟树 提示=8 t=ape::rtree(提示) #使用fossilsim模拟化石 比率=2 f=模拟化石泊松(比率,t) #绘制完整输出 绘图(f,t) #绘制重建的输出 绘图(f,t,重建=真)
#使用fossilsim模拟分类法 β=0.5#对称物种形成概率 lambda.a=0.1#再生率 s=模拟分类法(t,beta,lambda.a) #绘制输出 plot(s,t,legend.position=“bottomright”)