f.nonconst.linear: Function to compute the abundance of new or old RNA at time t...

In grandR: Comprehensive Analysis of Nucleotide Conversion Sequencing Data

Function to compute the abundance of new or old RNA at time t for non-constant rates.

Description

The standard mass action kinetics model of gene expression arises from the differential equation df/dt = s(t) - d(t) f(t), with s(t) being the synthesis rate at time t, d(t) the degradation rate at time t and f0=f(0) (the abundance at time 0). Here, both s and d have the following form s(t)=so+sf \cdot t^{se}.

Usage

f.nonconst.linear(t, f0, so, sf, se, do, df, de)

Arguments

t	time in h (can be a vector)
f0	the abundance at time t=0
so	synthesis date offset
sf	synthesis date factor
se	synthesis date exponent
do	degradation rate offset
df	degradation rate factor
de	degradation rate exponent

Value

the RNA abundance at time t

See Also

f.nonconst

grandR documentation built on April 4, 2025, 2:27 a.m.