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" or "Plane"

Examples

x <- rkent(100, mu = Line(120, 50), k = 5, b = 1)
kent_mle(x)
#> $G
#> Line object (n = 3):
#>        azimuth   plunge
#> mean  117.1742 55.15746
#> major 331.6208 29.85831
#> minor 231.9722 16.27603
#> 
#> $param
#>      kappa       beta        psi 
#>  3.6427651 -0.4280853 -0.6593137 
#> 
#> $logcon
#> [1] 4.201601
#> 
#> $loglik
#> [1] -154.1196
#> 
vmf_mle(x)
#> $loglik
#> [1] -154.9082
#> 
#> $mu
#> Line object (n = 1):
#>   azimuth    plunge 
#> 117.17419  55.15746 
#> 
#> $kappa
#> [1] 3.6065
#>