rej_samp

Rejection sampling from the (pseudo) pdf <code>f</code> in the interval between <code>x_min</code> and <code>x_max</code>. Returns <code>n</code> samples. Note that values of <code>f</code> below 0 are capped to zero

The fossil record is a joint expression of ecological, taphonomic,
evolutionary, and stratigraphic processes (Holland and Patzkowsky, 2012, ISBN:978-0226649382).
This package allowing to simulate biological processes in the time domain
(e.g., trait evolution, fossil abundance), and examine how their expression
in the rock record (stratigraphic domain) is influenced based on
age-depth models, ecological niche models, and taphonomic effects.
Functions simulating common processes used in modeling trait evolution or
event type data such as first/last occurrences are provided and can be used
standalone or as part of a pipeline. The package comes with example
data sets and tutorials in several vignettes, which can be used as a
template to set up one's own simulation.

Niklas Hohmann

StratPal

Stratigraphic Paleobiology Modeling Pipelines

rej_samp function

<dl><dt>f</dt>
<dd>function. (pseudo) pdf from which the sample is drawn</dd>
<dt>x_min</dt>
<dd>number, lower limit of the examined interval</dd>
<dt>x_max</dt>
<dd>number, upper limit of the examined interval</dd>
<dt>n</dt>
<dd>integer. number of samples drawn</dd>
<dt>f_max</dt>
<dd>number, maximum value of <code>f</code> in the interval from <code>x_min</code> to <code>x_max</code>. If f attains values larger than <code>f_max</code> a warning is throw, <code>f_max</code> is adjusted, and sampling is started again</dd>
<dt>max_try</dt>
<dd>maximum number of tries in the rejection sampling algorithm. If more tries are needed, an error is thrown. If this is the case, inspect of your function <code>f</code> is well-defined and positive, and if <code>f_max</code> provides a reasonable upper bound on it. Adjust <code>max_try</code> if you are certain that both is the case, e.g. if <code>f</code> is highly irregular.</dd></dl>

Arguments

random numbers from rejection sampling — rej_samp

<dl>

<dt>f</dt>
<dd>function. (pseudo) pdf from which the sample is drawn</dd>


<dt>x_min</dt>
<dd>number, lower limit of the examined interval</dd>


<dt>x_max</dt>
<dd>number, upper limit of the examined interval</dd>


<dt>n</dt>
<dd>integer. number of samples drawn</dd>


<dt>f_max</dt>
<dd>number, maximum value of <code>f</code> in the interval from <code>x_min</code> to <code>x_max</code>. If f attains values larger than <code>f_max</code> a warning is throw, <code>f_max</code> is adjusted, and sampling is started again</dd>


<dt>max_try</dt>
<dd>maximum number of tries in the rejection sampling algorithm. If more tries are needed, an error is thrown. If this is the case, inspect of your function <code>f</code> is well-defined and positive, and if <code>f_max</code> provides a reasonable upper bound on it. Adjust <code>max_try</code> if you are certain that both is the case, e.g. if <code>f</code> is highly irregular.</dd>

</dl>

rej_samp: random numbers from rejection sampling

Description

Usage

Value

Arguments

See Also

Examples