groupFeatures-abundance-correlation | R Documentation |
Features from the same originating compound are expected to have similar
intensities across samples. This method thus groups features based on
similarity of abundances (i.e. feature values) across samples in a
data set.
See also AbundanceSimilarityParam()
for additional information and
details.
This help page lists parameters specific for xcms
result objects (i.e.
XcmsExperiment()
and XCMSnExp()
objects). Documentation of the
parameters for the similarity calculation is available in the
AbundanceSimilarityParam()
help page in the MsFeatures
package.
## S4 method for signature 'XcmsResult,AbundanceSimilarityParam'
groupFeatures(
object,
param,
msLevel = 1L,
method = c("medret", "maxint", "sum"),
value = "into",
intensity = "into",
filled = TRUE,
...
)
object |
|
param |
|
msLevel |
|
method |
|
value |
|
intensity |
|
filled |
|
... |
additional parameters passed to the |
input object with feature group definitions added to (or updated
in) a column "feature_group"
in its featureDefinitions
data frame.
Johannes Rainer
feature-grouping for a general overview.
Other feature grouping methods:
groupFeatures-eic-similarity
,
groupFeatures-similar-rtime
library(MsFeatures)
library(MsExperiment)
## Load a test data set with detected peaks
faahko_sub <- loadXcmsData("faahko_sub2")
## Disable parallel processing for this example
register(SerialParam())
## Group chromatographic peaks across samples
xodg <- groupChromPeaks(faahko_sub, param = PeakDensityParam(sampleGroups = rep(1, 3)))
## Group features based on correlation of feature values (integrated
## peak area) across samples. Note that there are many missing values
## in the feature value which influence grouping of features in the present
## data set.
xodg_grp <- groupFeatures(xodg,
param = AbundanceSimilarityParam(threshold = 0.8))
table(featureDefinitions(xodg_grp)$feature_group)
## Group based on the maximal peak intensity per feature
xodg_grp <- groupFeatures(xodg,
param = AbundanceSimilarityParam(threshold = 0.8, value = "maxo"))
table(featureDefinitions(xodg_grp)$feature_group)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.