optimal.prob {untb} | R Documentation |
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
optimal.prob(x, interval=NULL, N=NULL, like=NULL, ...) optimal.theta(x, interval=NULL, N=NULL, like=NULL, ...)
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() |
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()
.
Robin K. S. Hankin
data(butterflies) optimal.prob(butterflies) optimal.theta(butterflies)