Skip to main content

Topics

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

Topics - jcmartin

1
XCMS - FAQ / xcms does not execute
I changed my computer and re-installed the whole packages, with the latest R version (3.5.1). Since then, xcms does not work anymore, remaining stacked at the first line

xset<-xcmsSet(method="centWave", peakwidth=c(3,20), snthresh=5, mzdiff=0.01, ppm=3, prefilter=c(5,100000), noise=20000, integrate=1)
xset

a warning message indicate :
Scanning files in directory C:/users/.... (blablabla until the correct file directory).. found 99 files

I also tried with the old R version that I installed and that used to work, but no success either

I'm not a beginner with xcms, and I tried to process my files using another computer in the lab and everything is alright.

I'm using windows 10 with an dell computer equipped with a Xeon processor

any idea of what'sgoing wrong?

many thanks
jc
2
Data Analysis Interest Group / data filtering
Hi everybody,
I\'m performing large-scale LCMS metabolomic profiling, and I\'m looking for a R script that performs quality control-based robust LOESS signal correction to correct for analytical drift within and across analytical batches (such as in  Dunn, W. B., Broadhurst, D., Begley, P., Zelena, E., Francis-McIntyre, S., Anderson, N., Brown, M., Knowles, J. D., Halsall, A. et al. (2011) Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc 6: 1060-1083).
I would be extremely pleased if somebody could help,
thanks and see you in washington

jc martin
biomet platform
marseille, france
3
XCMS / unreliable peak intensity in report table
Hi there,
I'm facing an unresolved problem; we are using a qexactive + in dual mode (+/-) @30000 mass resolution. After conversion into mzXML, we discovered that some features intensity are eratically reported in the diffreport table .tsv, whereas both the EIC caption and also the raw EIC are excellent. See figures below as an illustration
I would really appreciate help since we tested many different scripts without any subtantial improvements
thank you

raw xcalibur EIC
[attachment=2:8i42r00r]xcalibur raw EIC.JPG[/attachment:8i42r00r]

XCMS EIC
[attachment=1:8i42r00r]xcms EIX.JPG[/attachment:8i42r00r]

diffreport xcms Table values (.tsv)
[attachment=0:8i42r00r]diffreport table EIC values.JPG[/attachment:8i42r00r]

Hereafter is the script we tested:
library(xcms)                           
library(Biobase)                           
library(multtest)                           
xset<-xcmsSet(method="centWave"    peakwidth=c(3   15)    snthresh=2.5    mzdiff=-0.005    ppm=2    prefilter=c(3   1500000)    noise=10000    integrate=1
xset                           
xset1<-retcor(xset    method="obiwarp"    profStep=0.1    plottype="deviation")                   
xset1                           
xset2<-group(xset1    method="density"    bw=10    mzwid=0.01    minfrac=0.75    minsamp=1)            
xset2                           
xset3<-retcor(xset2    method="obiwarp"    profStep=0.1    plottype="deviation")                   
xset3                           
xset4<-group(xset3    method="density"    bw=2    mzwid=0.01    minfrac=0.75    minsamp=1)            
xset4                           
xset5<-fillPeaks(xset4)                           
xset5                           
#Pour obtenir le rapport du résultat d'XCMS#                           
reporttab<-diffreport(xset5   ctrl   pool   StdMixSwitchPosNeg_20150414   5000   metlin=0   h=480   w=640)      
library(CAMERA)                           
diffreport <- annotateDiffreport(xset5)                           
write.csv(diffreport    file="StdMixSwitchPosNeg_20150414_diffreport.csv")

[attachment deleted by admin]
4
XCMS / fillpeaks problem
I'm using a q-exactive plus connected to a uplc through an ESI source. After using xcms, the differport matrix contains a lot (too much) of variables with zero intensities, mainly in the samples, less so in the QC's.
Here is the script I  use

xset<-xcmsSet(method="centWave", peakwidth=c(3,20), snthresh=5, mzdiff=0.01, ppm=3, prefilter=c(5,100000), noise=100000)
xset
xset2<-retcor(xset, method="obiwarp", profStep=0.1, plottype="deviation")
xset2
xset2<- group(xset2, bw=5, mzwid=0.008)
xset3<-fillPeaks(xset2)
xset3

 I shall say that we previoulsy used a QTOF (Bruker microTOF Q) with of course a different xcms tuning, but never got such trouble.

Any ideas of what's going wrong with fillpeaks using  he q-exactive data?
5
XCMS / how to discard fillpeaks function
Hi,
I would like to process my data but not using fillpeaks. Hereafter is the return I got:

> reporttab<-diffreport(xset2,"échantillons","pools","urines_par_ech",10,metlin=0.15,h=480,w=640)
Erreur dans .local(object, ...) : NA values in xcmsSet. Use fillPeaks()
any idea how to circumvent this problem and to get to the reporttab step?

many thanks

jc