I已收到用于在Rstudio中绘制,检测峰并分析NMR光谱的代码。它几乎可以达到我想要的效果,但是在采摘高峰时,我很挣扎。该代码包含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')