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

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