solve.shift: Expected shifts ODE solver

In jasonxu90/bdsem: Analysis of discretely observed Birth-Death-Shift processes

Description

Numerically solves the ODE defined in de.shift using zvode from package deSolve given evaluation time, initial state values s1 and s2, and birth/shift/death rates

Usage

  solve.shift(time, dt, s1, s2, r, lam, v, mu)

Arguments

time	A number corresponding to the desired evaluation time of ODE
dt	A number giving the increment length used in solving the ODE
s1	A complex number giving the initial value of the ODE G
s2	A complex number
lam	Birth rate
v	Shift rate
mu	Death rate
r	A real number

Value

The function value of the shifts generating function

Examples

time = 5;  dt = 1; s1 = exp(2*1i); s2 = exp(3*1i); lam = .5; v = .2; mu = .4; r = 2
solve.shift(time,dt,s1,s2, r, lam,v,mu)

jasonxu90/bdsem documentation built on May 18, 2019, 5:54 p.m.