Nothing
td_loglik_bw_rhs = function(t,x,pars)
{
lx = length(x)
la = pars[1]
mu = pars[2]
K = pars[3]
ddep = pars[length(pars)]
r = (ddep == 5) * pars[4]
n0 = (ddep == 2 | ddep == 4)
nn = 0:(lx + 1)
lnn = length(nn)
p = c(0,x[1:lx],0)
if(ddep == 1)
{
lavec = pmax(rep(0,lnn),la - (la-mu)/K * nn)
muvec = mu * rep(1,lnn)
} else {
if(ddep == 2 | ddep == 2.1 | ddep == 2.2)
{
y = -(log(la/mu)/log(K+n0))^(ddep != 2.2)
lavec = pmax(rep(0,lnn),la * (nn + n0)^y)
muvec = mu * rep(1,lnn)
} else {
if(ddep == 3)
{
lavec = la * rep(1,lnn)
muvec = mu + (la - mu) * nn/K
} else {
if(ddep == 4 | ddep == 4.1 | ddep == 4.2)
{
y = (log(la/mu)/log(K+n0))^(ddep != 4.2)
lavec = la * rep(1,lnn)
muvec = mu * (nn + n0)^y
} else {
if(ddep == 5)
{
lavec = pmax(rep(0,lnn),la - 1/(r+1)*(la-mu)/K * nn)
muvec = muvec = mu + r/(r+1)*(la-mu)/K * nn
}}}}}
dx = lavec[(2:(lx+1))] * nn[(2:(lx + 1))] * p[(2:(lx + 1))+1] + muvec[(2:(lx+1))] * nn[(2:(lx+1))] * p[(2:(lx+1))-1] - (lavec[(2:(lx+1))-1] * nn[(2:(lx+1))-1] + muvec[(2:(lx+1))+1] * nn[(2:(lx+1))+1]) * p[(2:(lx+1))]
dG = p[2]
return(list(c(dx,dG)))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.