exp2flux (version 0.1)

fluxDifferences: Report the fold change of fluxes between two models

Description

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

Usage

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.

Examples

## 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)