CellNOptR (version 1.18.0)

simulatorTN: Simulation of a boolean model at any time points dependent on a previous one.

Description

This is a the simulator at TN using a C implementation. The computation relies on the time TN-1. T1 is a special case that is solved by using simulatorT2.

Usage

simulatorTN(simResultsPrev, CNOlist, model, simList, indexList, timeIndex=3)

Arguments

simResultsPrev

a matrix that is the output of simulatorTN (i.e. one row per condition and one column per species IN THE MODEL)

CNOlist

a CNOlist

model

a model that only contains the reactions to be evaluated

simList

a simList as created by prep4sim, that has also already been cut to contain only the reactions to be evaluated

indexList

an indexList as created by indexFinder

timeIndex

3 by default to behave like deprecated function simulatorT2. This is the timeIndex at which the simulation is requested.

Value

This function outputs a single matrix of format similar to valueSignals in the CNOlist but that contains an output for each species in the model. This matrix is the simulated equivalent of valueSignals at time 1, if you consider only the columns given by indexSignals.

Details

Differences from the BoolSimEngMKM simulator include: the valueInhibitors has not been previously flipped; the function outputs the values across all conditions for all species in the model, instead of only for the signal species. This is because then the output of this function can be used as initial values for the version of the simulator that works on time point 2 (not implemented in this version).

If you would like to compute the output of a model that contains some of the gates in the model but not all, we suggest that you use the function SimulateT1 and specify in the bStringT1 argument which gates you want to be included. Indeed, SimulateT1 is a wrapper around simulatorT1 that takes care of cutting the model for you before simulating it.

References

1. J. Saez-Rodriguez, L. G. Alexopoulos, J. Epperlein, R. Samaga, D. A. Lauffenburger, S. Klamt and P. K. Sorger. Discrete logic modeling as a means to link protein signaling networks with functional analysis of mammalian signal transduction, Molecular Systems Biology, 5:331, 2009.
2. M. K. Morris, J. Saez-Rodriguez, D. Clarke, P. K. Sorger, D. A. Lauffenburger. Training Signaling Pathway Maps to Biochemical Data with Constrained Fuzzy Logic: Quantitative Analysis of Liver Cell Responses to Inflammatory Stimuli, PLoS Comp. Biol., 7(3): e1001099, 2011.

See Also

simulateTN, cutAndPlotResultsT1

Examples

#This computes the output of the full model, which is normally not done on a stand alone basis, but if you have a model and would like to visualise its output compared to your data, then this is what you should do

data(CNOlistToy,package="CellNOptR")
data(ToyModel,package="CellNOptR")

indicesToy<-indexFinder(CNOlistToy,ToyModel,verbose=TRUE)
ToyFields4Sim<-prep4sim(ToyModel)

#Sim<-simulatorTN(
#	CNOlist=CNOlistToy,
#	model=ToyModel,
#	simList=ToyFields4Sim,
#	indexList=indicesToy)