This function permits estimation of total diversity based on a sample
frequency count table. Unlike breakaway
, it does not require
an input for the number of species observed once, making it an excellent
exploratory tool for microbial ecologists who believe that their sample may
contain spurious singletons. The underlying estimation procedure is similar
to that of breakaway
and is outlined in Willis & Bunge (2014).
The diversity estimate, standard error, estimated model coefficients and
plot of the fitted model are returned.
breakaway_nof1(
input_data,
output = NULL,
plot = NULL,
answers = NULL,
print = NULL
)
An object of class alpha_estimate
A category representing algorithm behaviour. code=1 indicates no nonlinear models converged and the transformed WLRM diversity estimate of Rocchetti et. al. (2011) is returned. code=2 indicates that the iteratively reweighted model converged and was returned. code=3 indicates that iterative reweighting did not converge but a model based on a simplified variance structure was returned (in this case, the variance of the frequency ratios is assumed to be proportional to the denominator frequency index). Please peruse your fitted model before using your diversity estimate.
The ``name'' of the selected model. The first integer represents the numerator polynomial degree and the second integer represents the denominator polynomial degree. See Willis & Bunge (2014) for details.
Estimated model parameters and standard errors.
The estimate of total (observed plus unobserved) diversity.
The standard error in the diversity estimate.
The chosen nonlinear model for frequency ratios.
An input type that can be processed by convert()
Deprecated; only for backwards compatibility
Deprecated; only for backwards compatibility
Deprecated; only for backwards compatibility
Deprecated; only for backwards compatibility
Amy Willis
Willis, A. (2015). Species richness estimation with high diversity but spurious singletons. arXiv.
Willis, A. and Bunge, J. (2015). Estimating diversity via frequency ratios. Biometrics.
breakaway
; apples
breakaway_nof1(apples[-1, ])
breakaway_nof1(apples[-1, ], plot = FALSE, output = FALSE, answers = TRUE)
Run the code above in your browser using DataLab