【发布时间】：2015-04-07 07:11:18
【问题描述】：
我有一个包含 135 个物种、6 个生活史特征（响应）和 2 个环境变量（预测变量）的数据集。我有兴趣使用 Caper 包进行 pgls 回归分析，以查看环境变量对每个生活史特征的影响。由于其中一个特征是女性体重，我计算了每个性状对女性体重的 pgls 回归，并在针对环境变量的多元回归中使用了随后的残差。为了检查回归诊断在针对体重的单变量特征回归中是否正常，我使用了 plot.pgls()。这些图显示了一些异常值的存在。在 Caper 中使用 pgls 的一些论文中，学生化残差 >3 的数据点已被排除为异常值。
- 谁能帮我弄清楚如何从 Caper 中的 pgls 回归计算学生化残差？
- pgls 回归给出系统发育和非系统发育残差。其中哪一个是学生化残差计算的？
感谢您的时间和帮助， pb
【问题讨论】：
-
你能举出一个可重现的小例子吗？
标签： r regression phylogeny