Skip to main content

Show Posts

This section allows you to view all Show Posts made by this member. Note that you can only see Show Posts made in areas you currently have access to.

Messages - cpataki

1
XCMS / Re: peak group detected with `peakDensityParam` but not listed in feature matrix
Thank you for the reference.

As a follow up, I have a very hard time separating these peaks and I'm hoping you can shed some light on what I could do differently.

My problem:
The peaks getting grouped together have an m/z range of: 153.03701335740, 153.054790674362 --> which is much larger than the 10ppm tolerance that is acceptable in our system to group peaks together. The peaks in that m/z range overlap quite a bit in the 8second window that is defined for that feature, so I can see how xcms could be grouping them together. However, I need the peaks in that feature to be separated into two groups. one group having the mz range: 153.037 - 153.053 and the other group at: 153.0535 - 153.055

I'm using PeakDensityParam and `groupChromPeaks(xset, param=gparam)` to group peaks, and I've tried a range of values for bw, binsize, and maxFeatures with no success. Am I on the right path? Is there anything else you suggest?
3
XCMS / peak group detected with `peakDensityParam` but not listed in feature matrix
I used the code below on my data, and I detect the following peak group (see screenshot of peaks). However, when I export the feature details, this feature group is not listed (see screenshot of filtered feature definitions). The closest feature group detected is well outside of the retention time window of the peaks plotted. Any help or explanation would be greatly appreciated.

rtr_1 <- c(155, 175)
mzr <- c(153.05351119538767, 153.05504173815234)
gparam <- PeakDensityParam(sampleGroups=fgroups,
                            minFraction = 1/nrow(filegroups),
                            bw = 3,
                            binSize = .025,
                            maxFeatures = 100)
chr <- chromatogram(xset5, rt = rtr_1 , mz = mzr)
plot(chr)#, col = paste0(cols, 80))
highlightChromPeaks(xset5, rt = rtr_1, mz = mzr)
par(mar = c(4.5, 4.5, 1, 0))
plotChromPeakDensity(xset5, mz = mzr, param = gparam, xlim = rtr_1)



[attachment deleted by admin]