时间戳

首页 > 解决方案 > R 中回归之间特定 X 值的置信区间重叠百分比

r - R 中回归之间特定 X 值的置信区间重叠百分比

问题描述

我的问题可能有点难以理解，但我现在会尝试进一步解释。但首先，我们以数据集为例mtcars。我们现在ggplot：

ggplot(mtcars) + 
  geom_jitter(aes(disp,mpg), colour="blue") + geom_smooth(aes(disp,mpg), method=lm, formula= y~poly(x,3), level=0.95, alpha=0.2) +
  geom_jitter(aes(hp,mpg), colour="green") + geom_smooth(aes(hp,mpg), method=lm, formula= y~poly(x,3), level=0.95, alpha=0.2) +
  labs(x = "Percentage cover (%)", y = "Number of individuals (N)")

我们得到了这个情节：

在此处输入图像描述

如您所见，两组的置信限的某些部分重叠。但是，我对重叠区域不感兴趣。

我对一条线上的重叠百分比非常感兴趣。在这种情况下，线是 X 的特定值。

让我们来吧x = 250。在下一张图片中，您会看到两个距离（红色和粉红色），每个距离对应于 1 组的 CL。红色和粉红色之间有重叠，我想计算两个组的百分比。如何计算这两个百分比？

在此处输入图像描述

标签： r

解决方案

如果没有一些数据可以举例说明，这有点棘手，所以，如果我理解正确，对于 x 轴上的给定值，在拟合模型后，您应该在某个地方同时具有区间的下限值和上限值。那么，如果 a 是上线，b 是下线，对于每个 x 点，您可以拥有 (La,Ua) 和 (Lb,Ub)，即两个置信区间的下限值和上限值，那么相对于曲线 a 的重叠将是 100* (Ub - La) / (Ua - La)。以类似的方式，相对于 b 的重叠比例将为 100* (Ua - Lb) / (Ub - Lb)。在这两种情况下，“/”都是除号。

推荐阅读