R/matrixfisher.mle.R

In Directional: A Collection of Functions for Directional Data Analysis

Defines functions matrixfisher.mle

Documented in matrixfisher.mle

## Maximum likelihood estimates for Matrix Fisher parameter F (3X3)
matrixfisher.mle <- function(X) {
   N <- dim(X)[3]
   #Xbar <- apply(X, 1:2, sum)
   #svd( Xbar / N )
   Xbar <- rowMeans(X, dims = 2)
   svd( Xbar )
}

Try the Directional package in your browser