r - 将多个峰拟合到数据集并在 R 中提取单个峰信息

Question

我有许多数据集，这些数据集由不同海拔的要素组成。目前有1-30m每1m间隔的数据。绘制时，我的许多数据集都显示 3-4 个峰，表示高度层。

这是一个示例数据集:

Height <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30) Counts <-c(4000,2000,500,300,200,100,0,0,400,700,800,800,500,1000,1500,2000,2500,2200,1700,1100,500,0,0,1000,1500,2000,3000,4000,4000,2000)

我想对这些数据集拟合某种曲线函数，以确定“峰”的总数、峰中心位置(即高度)和峰宽。前段时间我可以通过使用 fityk 软件手动拟合多个高斯函数来执行这种分析，但是我想知道是否可以通过 R 自动执行这样的过程？

我已经探索了许多其他关于将峰值拟合到直方图的帖子，例如通过 mixtools 包，但是我不知道您是否可以提取单个峰值信息。

如果您能提供任何帮助，我们将不胜感激。

Answer 1

“如何将曲线拟合到我的数据”这个问题太宽泛了，因为有无数种方法可以做到这一点。它也可能更适合 https://stats.stackexchange.com/比这里。然而，来自基础 R 的 ksmooth 是基本平滑器的一个很好的起点:

plot(Height,Counts)
smoothCounts<-ksmooth(Height,Counts,kernel="normal",bandwidth=2)
dsmooth<-diff(smoothCounts$y)
locmax<-sign(c(0,dsmooth))>0 & sign(c(dsmooth,0))<0
lines(smoothCounts)
points(smoothCounts$x[locmax],smoothCounts$y[locmax],cex=3,c=2)