makeGrid.birth.partial: Evaluates partial derivative of expected births ODE over 2D...
In jasonxu90/bdsem: Analysis of discretely observed Birth-Death-Shift processes

Description Usage Arguments Value Examples

Numerically partially differentiates the solution given in solve.birth over a grid of input values s1, s2 at one fixed time, and r=1

1 2	makeGrid.birth.partial(time, dt, s1.seq, s2.seq, lam, v, mu)

`time`	A number corresponding to the desired evaluation time of ODEs
`dt`	A number giving the increment length used in solving the ODE
`s1.seq`	A vector of complex numbers; initial values of the ODE G
`s2.seq`	A vector of complex numbers as inputs of s2.seq
`lam`	Birth rate
`v`	Shift rate
`mu`	Death rate

A matrix of dimension length(s1.seq) by length(s2.seq) of the function values

time = 5;  dt = 5; lam = .5; v = .2; mu = .4
gridLength = 32
s1.seq <- exp(2*pi*1i*seq(from = 0, to = (gridLength-1))/gridLength)
s2.seq <- exp(2*pi*1i*seq(from = 0, to = (gridLength-1))/gridLength)
makeGrid.birth.partial(time,dt,s1.seq,s2.seq,lam,v,mu)