getBirthMeans.timeList: Compute expected birth over a grid at a list of evaluation...
In jasonxu90/bdsem: Analysis of discretely observed Birth-Death-Shift processes

Description Usage Arguments Value Examples

Transforms the output matrices from makeGrid.birth.r1 and makeGrid.birth.partial to a list of expected sufficient statistics matrices using fft. Does this at several input times given in tList. Returns a list of matrices, where each list entry corresponds to a time in tList. The i,j entry of each matrix corresponds to the number of expected births when process has i type 1 particles and j type 2 particles, beginning with initNum type 1 particles, at the end of corresponding time interval.

1 2	getBirthMeans.timeList(tList, lam, v, mu, initNum, s1.seq, s2.seq, dt)

`tList`	A vector of numbers corresponding to the desired evaluation times
`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
`initNum`	An integer giving the number of initial particles
`lam`	Per-particle birth rate
`v`	Per-particle shift rate
`mu`	Per-particle death rate

A list of matrices of dimension length(s1.seq) by length(s2.seq); each list entry corresponds to an evaluation time from tList

tList = c(1,2);  dt = 1; lam = .5; v = .2; mu = .4; initNum = 10
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)
getBirthMeans.timeList(tList, lam, v, mu, initNum, s1.seq, s2.seq, dt)