Description Usage Arguments Author(s) Examples
This functions calculates the fold change "(fluxModel2/fluxModel1)-1"
for the fluxes of two given metabolic models.
1 | fluxDifferences(model1, model2, foldReport = 2)
|
model1 |
A valid model for the |
model2 |
A valid model for the |
foldReport |
A threshold value to be reported. All reactions with a greater or equal fold change than the given threshold are reported. |
Daniel Camilo Osorio <dcosorioh@unal.edu.co>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | ## Not run:
# Loading a model
library("sybil")
library("Biobase")
data("Ec_core")
# Generating expressionSets
expressionData <- matrix(data = runif(3*length(Ec_core@allGenes),min = 1,max = 100),
nrow = length(Ec_core@allGenes),
dimnames = list(c(Ec_core@allGenes),c()))
expressionData <- ExpressionSet(assayData = expressionData)
# Applying exp2flux
Ec_coreGE <- exp2flux(model = Ec_core,
expression = expressionData,
missing = "mean")
# Evaluating Differences
fluxDifferences(model1 = Ec_core,
model2 = Ec_coreGE,
foldReport = 0.5)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.