如何仅在R中的两条曲线之间找到正面积

Question

我试图仅在两条曲线之间找到正区域，例如其中一条曲线较高且仅在那里。曲线是时间序列和水平线。这是它创建的情节的图片：https://i.stack.imgur.com/ZrPH4.png。我不想要水平线较高的区域，所以我不能只找到它们差异曲线下的区域。我尝试在 auc 函数中使用 threshold 参数：

auc((1:96), STLFMondayBlock1$mean[(1:96)], thresh = 11648.93) 但我得到了一个比没有阈值更高的值：没有阈值的 1035514 和有阈值的 1213062。

Answer 1

使用一个简化的例子：

x1 <- c(0, 0, rep(1, 4), rep(0, 4), rep(1, 4), 0, 0, 0)
x2 <- rep(0.5, length(x1))
d <- x1 - x2

a <- cumsum(ifelse(d > 0, d, 0))
tail(a, 1)
# 4

par(mar=c(2, 2, 1, 1))
matplot(cbind(x1, x2, a), type="s", lwd=2, col=1:3)
legend("topleft", legend="positive area between curves",
  col=3, lty=3, lwd=2, bty="n")

如您所见，如果您取两条曲线之间的差值并仅对正值求和，则得到两条曲线之间的总正面积。