r - 计算直方图或密度函数中的峰

Question

似乎已经有很多“密度函数峰值”线程，但是我看不到专门解决这一点的线程。如果我错过了，很抱歉重复。

我的问题：给定一个包含1000个值的矢量（附有样本），我想识别数据的直方图或密度函数中的峰。从下面的示例数据的图像中，我可以看到直方图中的峰，分别位于〜0、6200和8400。但是，我需要获得这些峰的精确值，最好是通过简单的过程，因为我有数千个这样的矢量处理。



我最初开始使用直方图输出本身，但无法获得任何峰值搜索命令以使其正常工作（例如，根本没有）。我什至不知道它将如何从splus2R包中获取peaks（）命令以在直方图对象或密度对象上工作。这仍然是我的偏爱，因为我想确定每个峰的最大频率的确切数据值（与密度函数值（略有不同）相对，但我也找不到）。

我会自己发布样本数据，但是在这里看不到一种方法（对不起，如果我只是想丢失的话）。

Answer 1

如果您的y值是平滑的（例如在示例图中），则应该可以很容易地找到这些峰：

peakx <- x[which(diff(sign(diff(y)))==-2)]