bdARsimCondEnd: Conditional Simulation of BD via Accept-Reject

Description Usage Arguments Details Value Author(s) Examples

View source: R/accRejSim.R

Description

Simulates linear birth-death processes conditional on observing the end time (or a series of discrete observations), via simple accept reject (ie marginal simulation and accepting if it has the right end state).

Usage

1
2
3
4
5
bdARsimCondEnd(Naccepted = NULL, Ntotal = NULL,Nmax=NULL,
 bd.PO = new("CTMC_PO_1", states = c(5, 7, 3),
 times = c(0, 0.4, 1)), L = 0.5, m = 0.7, nu = 0.4)
bdARsimCondEnd.1(Naccepted = NULL, Ntotal = NULL, Nmax=NULL,
T = 1.02, L = 0.3, m = 0.4, nu = 0.1, a = 8, b = 9)

Arguments

Naccepted

Number of accepted sims to have at the end. Naccepted overrides Ntotal. If you want to use Ntotal, Naccepted should be NULL. Note that the number of sims will be >= Naccepted, probably not exactly equal to Naccepted.

Ntotal

Number of marginal sims to do; no guarantee of how many sims you will get out, but a better guarantee of how long it will take. If it gets no sims, it returns list().

Nmax

Different than Ntotal; it works with Naccepted. The function quits when either it has Naccepted sims or when it has done Nmax attempts. If it hits the max,returns whatever has been simulated so far, possibly list() if nothing.

T

Length of time of the chain.

L

Linear Birth rate.

m

Linear death rate.

nu

Immigration rate.

bd.PO

For bdARsimCondEnd, this is a list of observations essentially; Either class "CTMC_PO_1" or the analogous list.

a

Starting state.

b

Ending state (when you have just one observation).

Details

Outputs a list of BDMC objects. If Naccepted is not NULL then the list will be at least Naccepted long.

Value

List of BDMC objects.

Author(s)

Charles Doss

Examples

1
2
bdARsimCondEnd.1(Naccepted=10);  #default parameters; simulates at least10.
bdARsimCondEnd.1(Ntotal=10);  #default parameters; maybe end with none.

DOBAD documentation built on May 2, 2019, 3:04 a.m.