maxlikbat: Find maximum likelihood estimates for the inverse Batschelet...

Description Usage Arguments Value

View source: R/BatscheletOptimization.R

Description

Find maximum likelihood estimates for the inverse Batschelet distribution

Usage

1
2
3
maxlikbat(x, likfunbat_fun = likfuninvbat, weights, fixed_mu = NA,
  fixed_kp = NA, fixed_lam = NA, init_kp = 1, init_lam = 0,
  max_its = 20, kp_max = 1000)

Arguments

x

An set of angles in radians.

likfunbat_fun

A function that returns a function, which is the desired likelihood function to maximize over. Either likfuninvbat (the default) or likfunpowbat.

weights

A vector of length length(x), which gives importance weights to be used for x.

fixed_mu

If NA, the mu will be estimated as the mean direction. Else, mu will be fixed at fixed_mu.

fixed_kp

If NA, kp will be estimated. Else, a numeric giving the fixed value of kp.

fixed_lam

If NA, lam will be estimated. Else, a numeric giving the fixed value of lam.

init_kp

An initial value for κ, only used if neither κ or λ is fixed.

init_lam

An initial value for λ, only used if neither κ or λ is fixed.

max_its

Used only if neither fixed_kp nor fixed_lam is provided. The maximum number of iterations for the Nelder-Mead optimization.

kp_max

Used only if fixed_lam is provided. The maximum value for kappa to search for.

Value

The maximum likelihood estimates for the mu, kp, and lam.


keesmulder/flexcircmix documentation built on May 29, 2019, 3:02 a.m.