xcms icon indicating copy to clipboard operation
xcms copied to clipboard

How to reach to .. with XCMSnExp data?

Open linlennypinawa opened this issue 3 years ago • 3 comments

I used to load the data files to xcmsSet data container, and get to diffreport().
Could you tell me the workflow with XCMSnExp data to get the same report?

linlennypinawa avatar Aug 25 '22 23:08 linlennypinawa

the diffreport is not available for XCMSnExp objects. With the updated workflow it's now suggested to call quantify on the results XCMSnExp which returns a SummarizedExperiment object (standard data container from Bioconductor) and then use any other tool or method to perform the statistical analysis. There are way better approaches available in R to do the statistical tests then diffreport did.

As alternative, there is still the possibility to convert the XCMSnExp to a xcmsSet object and then call the diffreport on that. To convert:

xset <- as(xdata, "xcmsSet")

where xdata would be your XCMSnExp object.

jorainer avatar Aug 26 '22 13:08 jorainer

Thank you @jorainer. Would you please show me some resources of better approaches?

linlennypinawa avatar Aug 26 '22 20:08 linlennypinawa

One possibility would be Bioconductor's limma package - although it was developed originally for gene expression (microarray) data it works also good for metabolomics data. You can find the userguide here

jorainer avatar Sep 02 '22 12:09 jorainer

One possibility would be Bioconductor's limma package - although it was developed originally for gene expression (microarray) data it works also good for metabolomics data. You can find the userguide here

Fantastic!

linlennypinawa avatar Nov 17 '22 02:11 linlennypinawa