optimal.prob: Returns an estimate of the fundamental biodiversity number

optimal.probR Documentation

Returns an estimate of the fundamental biodiversity number

Description

Returns a maximum likelihood estimate for the fundamental biodiversity number \theta (function optimal.theta()) or the probability of mutation (function optimal.prob()) and optionally return information about the likely error

Usage

optimal.prob(x,  interval=NULL, N=NULL, like=NULL, ...)
optimal.theta(x, interval=NULL, N=NULL, like=NULL, ...)

Arguments

x

Ecosystem vector or species count table

interval

Bracketing interval for probability of mutation to be passed to the optimization routine (here optimize()). Default of NULL means to use a wide interval. Note that this argument is interpreted as an interval of \theta for both optimal.prob() and optimal.theta().

N

Integer; the number of parametric resampled estimates to give. Default of NULL means to return just the maximum likelihood estimate

like

Units of likelihood to calculate credible interval. Edwards recommends using 2

...

Further arguments passed to optimize()

Note

The fundamental biodiversity parameter \theta is 2\nu J, where \nu is the probability of mutation (ie, as estimated by optimal.prob()), and J is the size of the ecosystem.

For the general case of dispersal limitation, see functions etienne() and optimal.params().

Author(s)

Robin K. S. Hankin

See Also

etienne,optimal.params.sloss,optimal.params.gst

Examples

data(butterflies)
optimal.prob(butterflies)
optimal.theta(butterflies)

untb documentation built on Aug. 20, 2023, 9:08 a.m.