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}\).
f.nonconst(t, f0, s, d)
the RNA abundance at time t
time in h (can be a vector)
the abundance at time t=0
the synthesis rate (see details)
the degradation rate (see details)
Both rates can be either (i) a single number (constant rate), (ii) a data frame with names "offset", "factor" and "exponent" (for linear functions, see ComputeNonConstantParam; only one row allowed) or (iii) a unary function time->rate. Functions
f.nonconst.linear