MLE of spherical rotational symmetric distributions

Estimates the parameters of a von Mises-Fisher or Kent distribution.

Usage

kent_mle(x)

vmf_mle(x)

Source

Adapted from Directional::kent.mle() and Directional::vmf.mle()

Arguments

x

object of class "Vec3", "Line", "Ray", or "Plane", where the rows are the observations and the columns are the coordinates.

Examples

x <- rkent(100, mu = Line(120, 50), k = 5, b = 1)
kent_mle(x)
#> $G
#> Line object (n = 3):
#>         azimuth   plunge
#> mean  116.06250 54.39736
#> major  10.63589 10.78332
#> minor 273.40633 33.45465
#> 
#> $param
#>      kappa       beta        psi 
#> 4.00875580 0.41012182 0.05130077 
#> 
#> $logcon
#> [1] 4.470396
#> 
#> $loglik
#> [1] -145.04
#> 
vmf_mle(x)
#> $loglik
#> [1] -145.6211
#> 
#> $mu
#> Line object (n = 1):
#>   azimuth    plunge 
#> 116.06250  54.39736 
#> 
#> $kappa
#> [1] 3.96884
#>