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

如何使用getWaveletPeaks检测rstudio中的所有重叠峰,r,detect,overlapping,R,Detect,Overlapping,我收到了一个代码,可以在Rstudio中绘制、检测峰和分析核磁共振谱。它的工作原理和我想要的差不多,但在采摘高峰期我很挣扎。代码包括getWaveletPeaks函数,我正在研究如何检测非常接近的峰值,尤其是重叠峰值(所有局部最大值)。我可以用这个功能检测所有的局部最大点吗?在这种情况下,如何检测?我现在无法向您提供全部信息,似乎我只能添加一张图像?然而,在添加的图像中,你可以看到光谱的一个很好的部分,有几个峰,还有代表所有光谱检测到的局部最大值的点(总共29个重叠)。我的问题是,无论我似乎能够

我收到了一个代码,可以在Rstudio中绘制、检测峰和分析核磁共振谱。它的工作原理和我想要的差不多,但在采摘高峰期我很挣扎。代码包括getWaveletPeaks函数,我正在研究如何检测非常接近的峰值,尤其是重叠峰值(所有局部最大值)。我可以用这个功能检测所有的局部最大点吗?在这种情况下,如何检测?

我现在无法向您提供全部信息,似乎我只能添加一张图像?然而,在添加的图像中,你可以看到光谱的一个很好的部分,有几个峰,还有代表所有光谱检测到的局部最大值的点(总共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')
在此处输入代码

peaks我现在不能给你完整的信息,好像我只能添加一张图片?然而,在添加的图像中,你可以看到光谱的一个很好的部分,有几个峰,还有代表所有光谱检测到的局部最大值的点(总共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')
在此处输入代码
您能否提供您的问题的最小可重复性示例?也就是说,产生的数据+代码+输出,以及明确说明所需输出和您遇到的任何问题(例如,错误、与所需输出的偏差)。您能否提供一个最小的可重现性问题示例?也就是说,生成的数据+代码+输出,以及明确说明所需输出和遇到的任何问题(例如,错误、与所需输出的偏差)。