vMF_robust {scorematchingad}

Robust Fitting of von Mises Fisher

Description

Robust estimation for von Mises Fisher distribution using Windham().

Usage

vMF_robust(Y, cW, ...)

Arguments

Y

A matrix of observations in Cartesian coordinates.

cW

Tuning constants for each parameter in the vMF parameter vector. If a single number then the constant is the same for each element of the parameter vector.

...

Passed to Windham() and then passed onto vMF().

See Also

Other directional model estimators: Bingham(), FB(), vMF()

Other Windham robustness functions: Windham(), ppi_robust()

Examples

if (requireNamespace("movMF")){
  Y <- movMF::rmovMF(1000, 100 * c(1, 1) / sqrt(2))
} else {
  Y <- matrix(rnorm(1000 * 2, sd = 0.01), ncol = 2)
  Y <- Y / sqrt(rowSums(Y^2))
}
vMF_robust(Y, cW = c(0.01, 0.01), method = "smfull")
vMF_robust(Y, cW = c(0.01, 0.01), method = "Mardia")