starting.values: returns the null starting values used in model.test.sisters

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/EvoRAG_Code.R

Description

returns the null starting values used in model.test.sisters

Usage

1

Arguments

MODEL

any one of the evolutionary models implemented by model.test.sisters except BM_2rate and OU_2rate.

Details

because nlm (and other optimization functions in R) often gets stuck on local likelihood optima, model.test.sisters uses a variety of starting parameters in combination with nlm optimization and reports the results for the best set of starting parameters as the maximum likelihood estimate. The null starting parameters have been optimized for rates of Beta and Alpha that are typically less than 1. The user can also provide their own matrix of starting parameters. This function is provided here so the user can determine if the starting parameters are suited to their particular dataset. Alternatively, values of Euclidean distances and of L can also be transformed (i.e. by dividing large values by a constant) so they lie within a range acceptable for the starting parameters.

Value

returns a matrix with the starting values. Each column is a different parameter, and the last column is NA

Author(s)

Jason T. Weir

References

Weir JT, D Wheatcroft, & T Price. 2012. The role of ecological constraint in driving the evolution of avian song frequency across a latitudinal gradient. Evolution 66, 2773-2783.

Weir JT, & D Wheatcroft. 2011. A latitudinal gradient in rates of evolution of avian syllable diversity and song length. Proceedings of the Royal Society of London, B 278, 1713-1720.

See Also

model.test.sisters

Examples

1
2
starting.values(MODEL = "OU_linear")
    

EvoRAG documentation built on May 2, 2019, 8:57 a.m.