R语言恐龙祖先状态重建：尾椎数据构树及paleotree包故障问题

特征矩阵构建流程 #

• 先对150份样本的5个尾椎区段测量值做标准化校正，消除个体尺寸差异对特征值的干扰，推荐采用对数转换+系统发育异速生长残差的处理方式，同一物种的多份样本取均值作为该物种的特征代表值
• 矩阵结构设置为：行对应所有目标恐龙物种，列对应5个尾椎区段的测量特征，缺失值可通过missForest包进行插值，或在建树时开启缺失值兼容参数
• 如需引入外群类群校准树的拓扑结构，可直接将外群的公开尾椎测量数据补充到矩阵的对应行中

系统发育树构建方法 #

• 连续形态特征优先选择最大似然法或贝叶斯推断法建树：
  • 最大似然法可配合phangorn、phytools包实现，先基于特征矩阵计算距离矩阵，再拟合布朗运动或OU进化模型，通过自助法抽样计算节点支持率
  • 贝叶斯推断可使用MrBayes软件，设置连续特征进化参数，运行足够代数确认参数收敛后，即可得到带可信度的拓扑结构
• 若需要结合已有分子数据做联合建树，可将形态矩阵与分子序列矩阵合并后，通过BEAST2做联合推断

• 首先检查R版本适配性：paleotree最新版本要求R≥4.2.0，运行R.Version()$version.string即可查看当前版本，版本过低请先升级R到稳定版
• 强制重装完整依赖版本，运行以下命令：

install.packages("paleotree", dependencies = TRUE, type = "source")

• 若加载时出现包命名冲突报错，可直接通过paleotree::目标函数名的形式调用功能，或先执行detach("package:paleotree", unload = TRUE)卸载残留后重新加载
• 上述操作无效的情况下，可通过.libPaths()查看本地包安装路径，手动删除paleotree相关文件夹后重新执行安装命令即可

内容的提问来源于stack exchange，提问作者Thanushan Ravishankar