This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.
Messages - Jan Stanstrup
MZMine is know to use a lot of memory. I imagine that is where your bottleneck is. But you should check that.
XCMS is much more memory efficient. Be aware that each core will use a certain amount of memory. So on a system like yours not using all cores will use less memory and might save you if memory is your bottleneck. Also don't use 80 cores on processes that are bottlenecked by HDD reads (like reading the raw data).
That said, with 10,000 samples you really need to be careful about how greedy you need to be in terms of how low in intensity you want to pick.
An update to this: Latest version of msconvert should support the vendor centroiding. So the above method going through masslynx should no longer be necessary.
What values are you comparing? How do you get a single m/z value from the profile mode data to compare to?
So there is the profile mode data, Waters centroided m/z and the msconvert centroided m/z. The last two will be different due to different centroiding algorithms. The documentation says the CWT method is not very good. You could use Waters centroiding (if that is the one that is good?) if you centroid in masslynx first (to new raw file) and then convert without any additional processing.
Alternatively the R package MSnbase might have more advanced alternatives: https://github.com/lgatto/MSnbase/blob/11c336ebdc3e78cfa404803eb907346b046cd38b/vignettes/v03-MSnbase-centroiding.Rmd
Please do not post in unrelated topics. I Have split your topic here.
That said the XCMS online people don't seem to frequent this forum so you might have better luck contacting them directly.
Some example data would probably make it a lot easier to help.
What is the error?
If you post the output of
it might give a hint about the problem.
One comment: With XCMS3 you make an object with the raw data before the XCMS object. Instead of re-reading the raw files in your function you could reuse this object.
Have a look at these places:
I am not aware of any R focused course for metabolomics, but there are many general online R courses.
Here are some more options. I have only tried the one from Wehrens that seems to work just fine for me.
center is the correct argument for retcor.obiwarp.
XCMS is transitioning to a new interface (new functions with different arguments, but that basically does the same) where centerSample is an argument to ObiwarpParam. That it doesn't work with center=NULL but does when you don't specify might be a bug; it should work I think.
I think the reason you don't get a plot is because you have specified plottype = c("none","deviation"). You should only use one of the two. Otherwise the first option is used i.e. none.
Yes, default values are used if you don't specify explicitly.
family is a parameter for retcor.peakgroups not for retcor.obiwarp (e.g. method="obiwarp"). See the help of the individual methods.
It means the "files" variable is empty (or technically that it contains one string of length zero).
There are several possible errors in your code.
The double forward slash might be causing trouble if you are sure the path is otherwise correct. Try this (I think this might be a non-issue though):
path <- "G:/mzxml/np/neg"
"system.file" is meant to find files included with some package. Since you are using your own files you should be using "list.files" instead. Also keep in mind that the pattern is case sensitive.
If you write "files" (without quotes) in the console you should see the list of your files.
Variables are kept in memory. They are not written to files. Anything you want written to a file you need to do it explicitly.
An external drive should not be a problem, no.
Btw.: Best practice is to create a project in Rstudio and have a dedicated folder for each project. You then keep your scripts (and smaller support files if you have any) here.
Check what is in "files". I am guessing something strange got in there or it is empty.
XCMS online faq would suggest yes https://xcmsonline.scripps.edu/docs/fileformats.html.
In XCMS in R you cannot use vendor formats.
I would however be careful with waters files. Whether or not the lockmass correction is applied seems to depends on the version of masslynx and the msconvert version... So since it is unknown how XCMS online converts I would make sure my files are OK converted to mzML and upload that.
Your problem is probably the polarity switching.
I suggest you try splitting to two files when you do the conversion of your raw data. You should be able to do that with the polarity filter in msconvert.