
MLE of the ESAG distribution.
ESAGmle(y, tol = 1e-07)
A matrix with the data expressed in Euclidean coordinates, i.e. unit vectors.
The tolerance to accept that the E-M algorithm used to estimate the concentration parameter has converged.
A list including:
The mean vector in
The two gamma parameters.
The log-likelihood value.
The log-likelihood value of the isotropic angular Gaussian distribution. That is, the projected normal distribution which is rotational symmetric.
MLE of the MLE of the ESAG distributiontribution, on the sphere, is implemented. ESAG stands for Elliptically Symmetric Angular Gaussian and it was suugested by Paine et al. (2017). Unlike the projected normal distribution this is rotationally symmetric and is a competitor of the spherical Kent distribution (which is also non rotational symmetric).
Mardia, K. V. and Jupp, P. E. (2000). Directional statistics. Chicester: John Wiley & Sons.
Paine P.J., Preston S.P., Tsagris M. and Wood A.T.A. (2017). An Elliptically Symmetric Angular Gaussian Distribution. Statistics and Computing (To appear).
# NOT RUN {
m <- colMeans( as.matrix( iris[,1:3] ) )
y <- ESAGsim(1000, c(m, 1,0.5) )
ESAGmle(y)
# }
Run the code above in your browser using DataLab