Description Usage Arguments Details Value See Also Examples
Defines the configuration of an analysis that can be performed using geneSetAnalysis
, and returns it as a wrapper object.
1 2 3 4 5 6 7 8 9 10 11 12 | gsAnalysis(name,
gls = NULL,
glsParameterNames = NULL,
transformation = NULL,
transformationParameterNames = NULL,
gss = NULL,
gssParameterNames = NULL,
globalStat = NULL,
globalStatParameterNames = NULL,
significance = NULL,
significanceParameterNames = NULL,
testAlternative = c("greater", "less"))
|
name |
A character string describing the analysis. |
gls |
The name of the function that calculates the gene-level statistic for a given dataset. If set to |
glsParameterNames |
A character vector of names of the parameters used by the gene-level statistic defined in |
transformation |
The name of the function that transforms the gene-level statistics values. If set to |
transformationParameterNames |
A character vector of names of the parameters used by the transformation defined in |
gss |
The name of the function that calculates the gene set statistics from untransformed or transformed gene-level statistic values. If set to |
gssParameterNames |
A character vector vector of names of the parameters used by the gene set statistic defined in |
globalStat |
If the gene set analysis consists of a global test and cannot be divided into the calculation of a gene-level statistic, a transformation and the calculation of the gene set statistics, this parameter defines the name of the function that performs this global analysis. In this case, the parameters |
globalStatParameterNames |
A character vector names of the parameters used by the global analysis defined in |
significance |
The name of a method that performs a significance assessment for the gene set statistic values. If set to |
significanceParameterNames |
A character vector of names of the parameters used by the significance assessment method defined in |
testAlternative |
Specifies the alternative hypothesis of the significance test for the gene set enrichment, which may be dependent on the chosen gene set statistic. Must be one of |
The function provides a way of flexibly defining the steps of the gene set analysis pipeline. This pipeline consists of a subset of the following steps, each of which may have specific parameters:
Gene-level analysis:
A gene-level statistic scores the relationship between the measurements for a specific gene and the class labels. Typical measures include correlation coefficients, the t statistic or the fold change between the groups (see gls
for gene-level statistics included in the package).
Transformation of gene-level statistics:
Optionally, the gene-level statistic values can be postprocessed, e.g. by taking the absolute value or the square for correlation values or by binarizing or ranking values. See transformation
for transformations included in the package.
Gene set analysis:
Based on the (possibly transformed) gene-level statistics, the gene set(s) of interest is/are scored. Examples are the median, the mean or the enrichment score of the gene-level statistic values in the gene set(s). See gss
for gene set statistics included in the package.
Significance assessment:
To assess the significance of the gene set statistic value(s) with respect to a null distribution, computer-intensive tests are performed. These tests repeatedly sample random label vectors or gene sets and calculate their gene set statistic values. These values can then be compared to the true gene set statistics. See significance
for significance assessment methods included in the package.
Global analysis:
As an alternative to the above pipeline steps, it is possible to define a single, global method that directly calculates an enrichment p-value for a supplied data set and gene set. See global
for the global analysis tests included in the package.
Several state-of-the-art analyses have predefined configuration objects in which the above steps are defined accordingly (see predefinedAnalyses
).
An object of class gsAnalysis
with components corresponding to the above parameters.
predefinedAnalyses
, geneSetAnalysis
, evaluateGeneSetUncertainty
, gls
, transformation
, gss
, global
, significance
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | # defines an analysis that corresponds to gsAna1()
gsa <- gsAnalysis(
name = "averageCorrelation",
gls = "gls.cor",
glsParameterNames = c("labs","method"),
transformation = "transformation.abs",
transformationParameterNames = NULL,
gss = "gss.mean",
gssParameterNames = NULL,
globalStat = NULL,
globalStatParameterNames = NULL,
significance = "significance.sampling",
significanceParameterNames = c("numSamples"),
testAlternative = "greater")
print(gsa)
# load data
require(GlobalAncova)
data(vantVeer)
data(phenodata)
data(pathways)
# apply the previously defined analysis
res <- geneSetAnalysis(
# global parameters
dat = vantVeer,
geneSets = pathways[1],
analysis = gsa,
# parameters for the specific analysis gsAna1
labs = phenodata$metastases,
numSamples = 10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.