FM.replicate: Replicate 'FM.run'

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

Description

Replicate FM.run

Usage

  FM.replicate(numReps, N, tList, lam, v, mu, initList,
    initGuess)

Arguments

numReps	The number of desired replications
N	An integer specifying the number of observation intervals/realization from the simple BSD process.
tList	A list of observation interval lengths. The number of datasets returned is equal to the length of tList
lam	Per-particle birth rate
v	Per-particle shift rate
mu	Per-particle death rate
initList	A vector containing possible initial population sizes
initGuess	Vector of numbers, initial guess for optim

Value

An array with entries of type returned by FM.run. Rows correspond to a dt value in tList, and columns correspond to replications

Examples

tList <- c(.2,.4,.6); initList <- c(1:15)
lam = .06; v = .02; mu = .11
trueParam <- c(lam,v,mu)
N = 50; numReps = 2
example <- FM.replicate(numReps,N,tList,lam,v,mu,initList, trueParam)

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