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

numeric. Can be three element vector, three column array, or an object of class "line" or "plane"

Examples

x <- rkent(100, mu = Line(120, 50), k = 5, b = 1)
kent.mle(x)
#> $G
#>         azimuth   plunge
#> [1,] 116.626856 54.87378
#> [2,]   6.078073 13.87059
#> [3,] 267.361734 31.53770
#> attr(,"class")
#> [1] "line"
#> 
#> $param
#>       kappa        beta         psi 
#>  5.15963154  1.11871818 -0.04061578 
#> 
#> $logcon
#> [1] 5.410598
#> 
#> $loglik
#> [1] -121.7777
#> 
vmf.mle(x)
#> $loglik
#> [1] -126.5135
#> 
#> $mu
#>       azimuth   plunge
#> [1,] 116.6269 54.87378
#> attr(,"class")
#> [1] "line"
#> 
#> $kappa
#> [1] 4.816487
#>