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

问题描述 投票:1回答:1

enter image description here I已收到用于在Rstudio中绘制,检测峰并分析NMR光谱的代码。它几乎可以达到我想要的效果,但是在采摘高峰时,我很挣扎。该代码包含getWaveletPeaks函数,并且我正在摸索着如何找出紧邻的峰,尤其是重叠的峰(所有局部最大值)。我可以使用此功能检测所有局部最大点吗,在这种情况下如何?

r detect overlapping
1个回答
0
投票

我现在无法为您提供全部信息,好像我只能添加一张图像?但是,在添加的图像中,您可以看到光谱中不错的部分,其中有几个峰,还显示了代表所有光谱中检测到的局部最大值的点(总共覆盖了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')
© www.soinside.com 2019 - 2024. All rights reserved.