fluxDifferences: Report the fold change of fluxes between two models

Description Usage Arguments Author(s) Examples

View source: R/fluxDifferences.R

Description

This functions calculates the fold change "(fluxModel2/fluxModel1)-1" for the fluxes of two given metabolic models.

Usage

1
fluxDifferences(model1, model2, foldReport = 2)

Arguments

model1

A valid model for the 'sybil' package.

model2

A valid model for the 'sybil' package. Must have the same reactions (reaction number and reaction identifiers) as "model1" with different restrictions.

foldReport

A threshold value to be reported. All reactions with a greater or equal fold change than the given threshold are reported.

Author(s)

Daniel Camilo Osorio <[email protected]>

Examples

 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)

exp2flux documentation built on May 30, 2017, 3:50 a.m.