我有一个包含 135 个物种、6 个生活史特征(响应)和 2 个环境变量(预测变量)的数据集。我有兴趣使用 Caper 包进行 pgls 回归分析,以查看环境变量对每个生活史特征的影响。由于其中一个性状是女性体重,我计算了每个性状对女性体重的 pgls 回归,并在针对环境变量的多元回归中使用了随后的残差。为了检查回归诊断在针对体重的单变量特征回归中是否正常,我使用了 plot.pgls()。这些图显示了一些异常值的存在。在 Caper 中使用 pgls 的一些论文中,学生化残差 > 3 的数据点已被排除为异常值。
- 有人可以帮我弄清楚如何从 Caper 的 pgls 回归中计算学生化残差吗?
- pgls 回归给出系统发育和非系统发育残差。其中哪一个是学生化残差计算的?
谢谢你的时间和帮助,pb