Fluxenrichmentanalysis

FEA(model, rxnSet, group)

Significance analysis - Flux enrichment analysis using hypergeometric 1-sided test and FDR correction for multiple testing

USAGE:

resultCell = FEA (model, rxnSet, ‘subSystems’)

INPUTS:

• model – COBRA model structure with fields:

  • .rxns - n x 1 reaction identifiers

  • .subSystems - n x 1 subsystem assignment of each reaction, used for enrichment

• rxnSet – reaction set to be enriched (vector of reaction indices e.g. 1:10)

• group – name of the model field tested for enrichment, e.g. ‘subSystems’ : FEA looks for significantly enriched subsystems in rxnSet

OUTPUT:

resultCell – cell structure of enriched groups

Example

load ecoli_core_model; resultCell = FEA(modelEcore, 1:10, ‘subSystems’);