如何用getWaveletPeaks检测rstudio中的所有重叠峰

Question

​

我已经收到了一个代码，可以在Rstudio中绘制、检测峰值和分析核磁共振谱。它就像我希望的那样工作，但在选择的高峰期，我很挣扎。这段代码包含了getWaveletPeaks函数，我让海量的trubble计算出如何检测接近的峰值，特别是重叠的峰值(所有局部最大值)。我可以用这个函数检测所有的局部最大值点吗?在这种情况下，如何检测？

Answer 1

我现在不能给你所有的信息，似乎我只能添加一张图片？然而，在添加的图像中，您可以看到具有多个峰值的光谱的一部分，还可以看到代表所有光谱检测到的局部最大值的点(总共29个重叠)。我的问题是，无论我似乎能够改变什么，都无法检测到一些峰值(例如2.055 ppm)。代码中我认为不能像我希望的那样工作的部分是：

enter code here

peaks <- speaq::getWaveletPeaks(Y.spec=Spectra,
                                 X.ppm=ppm, 
                                 window.width = "small",
                                 window.split = 64,                                            
                                 baselineThresh = 2,                                        
                                 SNR.Th = 7,                                           
                                 nCPU = 6,                                                 
                                 include_nearbyPeaks = TRUE)                                
save(peaks,file='BiGh_191009_neuro_csf_800_peaks_snr10_bt10.Rdata')