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,] 112.1300 54.548731
#> [2,] 204.0171  1.343155
#> [3,] 294.9725 35.417948
#> attr(,"class")
#> [1] "line"
#> 
#> $param
#>     kappa      beta       psi 
#> 4.7999605 1.2683185 0.3603214 
#> 
#> $logcon
#> [1] 5.144895
#> 
#> $loglik
#> [1] -129.1967
#> 
vmf.mle(x)
#> $loglik
#> [1] -136.1387
#> 
#> $mu
#>      azimuth   plunge
#> [1,]  112.13 54.54873
#> attr(,"class")
#> [1] "line"
#> 
#> $kappa
#> [1] 4.370741
#>