如何计算R中频率多边形下的面积？

Question

我正在努力寻找一种计算频率多边形下面积的方法，如下所示：

x1 <- 1:5                         
y1 <- c(0.2, 0.14, 0.7, 0.11, 0.1) 
plot(x1, y1, type = "l", lwd = 3)   
polygon(c(1, x1, 5), c(0, y1, 0), col = "red") 
points(x1, y1,cex = 2,pch = 15)
segments(x1, 0, x1, y1)

所以基本上红色区域的区域从 1 到 5.. 任何建议将不胜感激！ 非常感谢

Answer 1

如果你的多边形的底边总是0 并且左右两边是垂直的，那么它只不过是一组梯形的面积相加：

h = x1[2:5]-x1[1:4] # heights
b1 = y1[1:4] #first bases
b2 = y1[2:5] #second bases
ta = h*(b1+b2)/2 # trapezoids
pa = sum # total area

当然，您可以泛化索引以使其适用于不同长度的数组。