RPANDA (version 2.3)

remove.model: Removing a model from shift.estimates output

Description

Allows to remove a model from the model comparisons of shift.estimates output.

Usage

remove.model(shift.res, model)

Value

the same output resulting from shift.estimates but without the chosen model in model comparisons.

Arguments

shift.res

the output resulting from shift.estimates.

model

character. Specifies the model to remove from the set of model of diversification applied in shift.res.

Author

Nathan Mazet

Details

This function allow to remove model one at a time. The idea is to remove a model without having to reanalyse the phylogeny and all the combinations of shifts if a model (e.g. BVAR_DVAR) behaves strangely on the studied phylogeny.

References

Mazet, N., Morlon, H., Fabre, P., Condamine, F.L., (2023). Estimating clade‐specific diversification rates and palaeodiversity dynamics from reconstructed phylogenies. Methods in Ecology and in Evolution 14, 2575–2591. https://doi.org/10.1111/2041-210X.14195

See Also

shift.estimates

Examples

# loading data
data("shifts_cetacea")

# Removing "BVAR_DCST" model for the example
shifts_cetacea_noBVAR_DCST <- remove.model(shift.res = shifts_cetacea,
                                           model = "BVAR_DCST")