Run a topological analysis on an expression dataset using topologyGSA.

Share:

Description

This function is deprecated and will be removed in a future release. You can use runTopologyGSA instead.

Usage

1
  runTopologyGSAMulti(pathways, test, exp1, exp2, alpha, maxNodes=150, ...)

Arguments

pathways

a PathwayList object.

test

Either "var" and "mean". Determine the type of test used by topologyGSA.

exp1

Experiment matrix of the first class, genes in columns.

exp2

Experiment matrix of the second class, genes in columns.

alpha

Significance level of the test.

maxNodes

Ignore pathways with more than "maxNodes" nodes. Set to "NULL" to disable the filter.

...

Additional parameters for topologyGSA.

Details

This function produces a warning and whenever the number of genes in common between the expression matrices and a pathway is less than 3.

Value

A list with two elements:

  • results: a list with one entry for each successfully analyzed pathway;

  • errors: a vector containing the error messages of failed analyses.

References

Massa MS, Chiogna M, Romualdi C. Gene set analysis exploiting the topology of a pathway. BMC System Biol. 2010 Sep 1;4:121.

See Also

pathway.var.test pathway.mean.test

Examples

1
2
3
4
5
6
7
8
9
if (require(topologyGSA)) {
  data(examples)

  k <- pathways("hsapiens", "kegg")
  ps <- convertIdentifiers(
    k[c("Acute myeloid leukemia", "Fc epsilon RI signaling pathway")],
    "symbol")
  runTopologyGSAMulti(ps, "var", y1, y2, 0.05)
}