如何仅求R中两条曲线之间的正面积

Question

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

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")

​

正如你所看到的，如果你取两条曲线之间的差，只和正数，你就得到了这两条曲线之间的总正面积。