Description Usage Arguments Details Value References See Also Examples
View source: R/geneSetAnalysis.R
The main function of the package that performs a gene set analysis for a list of gene sets.
1 2 3 4 5 6 7 8 9 | geneSetAnalysis(
...,
dat,
geneSets,
analysis,
signLevel = 0.05,
preprocessGeneSets = FALSE,
adjustmentMethod = p.adjust.methods,
cluster = NULL)
|
... |
Additional parameters for the different steps of the analysis pipeline, depending on the concrete configuration supplied in |
dat |
A numeric matrix of gene expression values for all analyzed genes. Here, each row corresponds to one gene, and each column corresponds to one sample. The rows must be named with the gene names used in the gene sets. |
geneSets |
A list of gene sets, where each gene set is a vector of gene names corresponding to the row names of |
analysis |
An object of type |
signLevel |
The significance level for the significance assessment of the gene sets (defaults to |
preprocessGeneSets |
Specifies whether the gene sets in |
adjustmentMethod |
The method to use for the adjustment for multiple testing (see |
cluster |
If the analyses should be applied in parallel for the gene sets, this parameter must hold an initialized cluster as returned by |
This is the main interface function of the package for gene set enrichment analyses. Analyses usually consist of a pipeline of steps. Often, the first step is the calculation of a summary statistic for the relation of each gene to the class labels. These values or transformations thereof are employed to calculate a gene set statistic for each of the supplied gene sets. The significance of gene set enrichments can be determined according to different methods, and the robustness of gene sets can be evaluated by slightly modifying the gene sets. To provide a flexible mechanism for the plethora of different approaches arising from the different choices, basic pipeline configurations are encapsulated in gsAnalysis
objects which can be created using the gsAnalysis
function. Ready-to-use configuration objects for certain well-known methods are included in the package (see predefinedAnalyses
). Parameters of the chosen analysis pipeline can be set in the ...
parameter.
An object of the type gsaResult
with the following elements:
adjustedPValues |
A vector of p-values, one for each gene set. These values are already adjusted for multiple testing according to the |
rawPValues |
The raw unadjusted p-values, one for each gene set. |
res.all |
A list comprising the detailed results for each gene set. Each element of this list is another list with the following components:
|
signLevel |
The significance level used for this analysis. |
analysis |
The performed analysis (of type |
analysisType |
A character string identifying the analysis as an enrichment analysis pipeline ( |
adjustmentMethod |
The method used to adjust the p-values in |
Ackermann, M., Strimmer, K. (2009) A general modular framework for gene set enrichment analysis. BMC Bioinformatics, 10(1), 47.
gsAnalysis
, gls
, transformation
, gss
, global
, significance
, evaluateGeneSetUncertainty
, hist.gsaResult
, preprocessGs
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | # load data
require(GlobalAncova)
data(vantVeer)
data(phenodata)
data(pathways)
# apply predefined analysis for gene set enrichment analysis
res <- geneSetAnalysis(
# parameters for geneSetAnalysis
dat = vantVeer,
geneSets = pathways[1],
analysis = analysis.averageCorrelation(),
adjustmentMethod = "fdr",
# additional parameters for analysis.averageCorrelation
labs = phenodata$metastases,
method = "pearson",
numSamples = 10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.