dynamo.sim.Gillespie
- dynamo.sim.Gillespie(a=None, b=None, la=None, aa=None, ai=None, si=None, be=None, ga=None, C0=array([[0.0], [0.0], [0.0], [0.0], [0.0]]), t_span=[0, 50], n_traj=1, t_eval=None, dt=1, method='basic', verbose=False)[source]
A simulator of RNA dynamics that includes RNA bursting, transcription, metabolic labeling, splicing, transcription, RNA/protein degradation
- Parameters:
a (
Optional
[float
]) – rate of active promoter switches to inactive oneb (
Optional
[float
]) – rate of inactive promoter switches to active oneaa (
Optional
[float
]) – transcription rate with active promoterai (
Optional
[float
]) – transcription rate with inactive promoterga (
Optional
[float
]) – gamma, the fraction of labeled u turns to unlabeled sC0 (
ndarray
) – A numpy array with dimension of 5 x n_gene. Here 5 corresponds to the five species (s - promoter state, ul, uu, sl, su) for each gene.t_span (
List
) – list of between and end time of simulationn_traj (
int
) – number of simulation trajectory to uset_eval (
Optional
[float
]) – the time points at which data is simulateddt (
float
) – delta t used in simulationmethod (
str
) – method to simulate the expression dynamicsverbose (
bool
) – whether to report running information
- Returns:
an Annodata object containing the simulated data.
- Return type:
adata