Description Usage Arguments Value Examples
Transforms the output matrices from
makeGrid.particleT.r0
and
makeGrid.particleT.partial
to a list of
expected sufficient statistics matrices using fft
.
Does this at several initial particle counts from
initList Returns a list of matrices, where each list
entry corresponds to a number of initial particles The
i,j entry of each matrix corresponds to the expected
particle time spent with i type 1 particles and j type 2
particles, beginning with initNum type 1 particles, by
the end of corresponding time interval.
1 2 | getParticleT.initList(u, initList, lam, v, mu, s1.seq,
s2.seq, dt)
|
initList |
A vector of integers corresponding to the desired initial particle counts |
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 |
u |
A number giving the observation interval length, equivalently the time to evaluate the ODEs |
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 initial number of particles from initList
1 2 3 4 5 6 7 | initList = c(10,11)
#gives matrices of expected sufficient statistics corresponding to 10 initial particles and 11 initial particles
u = 1; dt = 1; 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)
getParticleT.initList(u, initList, lam, v, mu, s1.seq, s2.seq, dt)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.