simulatorDT: Discrete time simulation of a boolean model

Description Usage Arguments Value Author(s) See Also Examples

View source: R/simulatorDT.R

Description

Simulates multiple time points within C (for speed).

Usage

1
2
simulatorDT(CNOlist, model, simList, indices, boolUpdates,
prevSim=NULL)

Arguments

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.

indices

An indexList as created by indexFinder.

boolUpdates

The number of update rounds the simulator should run for.

prevSim

The results from simulatorT0 can be used here as initial conditions.

Value

A 3-dimensional array that gives the value of all species under each condition at each update (conditions, species, update).

Author(s)

A. MacNamara

See Also

gaBinaryDT, getFitDT

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
# 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.

library(CellNOptR)
library(CNORdt)
data(CNOlistPB, package="CNORdt")
data(modelPB, package="CNORdt")

indexOrig <- indexFinder(CNOlistPB, modelPB, verbose=TRUE)
fields4Sim <- prep4sim(modelPB)
boolUpdates=10

simResults <- simulatorDT(
  CNOlist=CNOlistPB,
  model=modelPB,
  simList=fields4Sim,
  indices=indexOrig,
  boolUpdates=boolUpdates
)
simResults = convert2array(simResults, dim(CNOlistPB$valueSignals[[1]])[1],
length(modelPB$namesSpecies), boolUpdates)

saezlab/CNORdt documentation built on April 22, 2018, 12:07 p.m.