fit_function function

Maximum Likelihood estimator for a Cauchy model

Find the maximum likelihood.

fit_function(
  minus_like,
  tree,
  trait,
  X,
  model,
  start.values,
  lower.values,
  upper.values,
  optim,
  rootTip
)

Arguments

• minus_like: a function giving the minus likelihood of the model.
• trait: named vector of traits at the tips.
• model: a model for the trait evolution. One of "cauchy" or "lambda" (see Details).
• optim: if "local", only a local optimization around the initial parameter values is performed (the default). If "global", a global maximization is attempted using the "MLSL" approach (see nloptr).
• rootTip: root tip for the reml.

Returns

A list, with the maximum likelihood rate parameter, and the likelihood value.