miniDiffReport | R Documentation |
Generates conventional XCMS-like diff report comparing peaks between two sample classes, e.g. "control", and "treatment" from an xcmsSet object configured for X13CMS (i.e. one in which the sample classes have to be demarcated as unlabeled and labeled)
miniDiffReport(xcmsSet, class1sampNames, class2sampNames, varEq = FALSE, intChoice)
xcmsSet |
XCMS object |
class1sampNames |
The subset of rownames(xcmsSet@phenoData) corresponding to names of the unlabeled samples of class 1 (e.g. "control"). |
class2sampNames |
The subset of rownames(xcmsSet@phenoData) corresponding to names of the unlabeled samples of class 2 (e.g. "perturb"). |
varEq |
Boolean indicating whether to assume that absolute ion intensities of each peak in both sample classes are drawn from distributions with equal variance. Defaults to FALSE. |
intChoice |
one of "maxo", "into", or "intb"–the choice of which peak intensity measurement to use from the XCMS object. |
A dataframe whose rows correspond to peak groups identified in the xcmsSet object and whose columns are: - pvalue: p-value of Welch's t-test comparing absolute intensities of the peak in the two conditions - foldChange: change in mean peak intensity in condition 2 compared to conditions 1 - mzmed: median m/z of peak - rtmed: median retention time of peak - means1: mean peak intensity in condition 1 - means2: mean peak intensity in condition 2 - peakIntensities1: all peak intensities in the group in condition 1 - peakIntensities2: all peak intensities in the group in condition 2
Xiaojing Huang
## Not run:
## diffReport = miniDiffReport(xcmsSet, class1sampNames, class2sampNames, varEq = FALSE, intChoice = "intb")
## From the example given in getIsoLabelReport(), the variables class1sampNames and class2sampNames would be c("control_unlabeled1", "control_unlabeled2"), and c("perturb_unlabeled1", "perturb_unlabeled2"), respectively.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.