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,] 111.2018 52.587689
#> [2,] 202.6413  1.100846
#> [3,] 293.4828 37.390396
#> attr(,"class")
#> [1] "line"
#> 
#> $param
#>     kappa      beta       psi 
#> 4.4381643 1.1139463 0.3305089 
#> 
#> $logcon
#> [1] 4.850385
#> 
#> $loglik
#> [1] -136.2913
#> 
vmf.mle(x)
#> $loglik
#> [1] -142.1328
#> 
#> $mu
#>       azimuth   plunge
#> [1,] 111.2018 52.58769
#> attr(,"class")
#> [1] "line"
#> 
#> $kappa
#> [1] 4.112692
#>