The mle algorithm for the matrix normal distribution

The mle algorithm for the matrix normal distribution The maximum likelihood estimation (MLE) of the parameters of the matrix normal distribution is considered. In the absence of analytical solutions of the system of likelihood equations for the among-row and among-column covariance matrices, a two-stage algorithm must be solved to obtain their maximum likelihood estimators. A necessary and sufficient condition for the existence of maximum likelihood estimators is given and the question of their stability as solutions of the system of likelihood equations is addressed. In particular, the covariance matrix parameters and their maximum likelihood estimators are defined up to a positive multiplicative constant; only their direct product is uniquely defined. Using simulated data undertwovariance-covariancestructures that, otherwise, are indistinguishable by semivariance analysis, further specific aspects of the procedure are studied: (1) the convergence of the MLE algorithm is assessed; (2) the empirical bias of the direct product ofcovariance matrix estimators is calculated for various sample sizes; and (3) the consistency of the estimator is evaluated by its mean Euclidean distance from the parameter, as a function of the sample size. The adequacy of the matrix normal model, including the separability of the variance-covariance structure, is tested on multiple time series of dental medicine data; other applications to real doubly multivariate data are outlined. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of Statistical Computation and Simulation Taylor & Francis

The mle algorithm for the matrix normal distribution

, Volume 64 (2): 19 – Sep 1, 1999
19 pages

/lp/taylor-francis/the-mle-algorithm-for-the-matrix-normal-distribution-8OfC5c1v4S

References (25)

Publisher
Taylor & Francis
Copyright Taylor & Francis Group, LLC
ISSN
1563-5163
eISSN
0094-9655
DOI
10.1080/00949659908811970
Publisher site
See Article on Publisher Site

Abstract

The maximum likelihood estimation (MLE) of the parameters of the matrix normal distribution is considered. In the absence of analytical solutions of the system of likelihood equations for the among-row and among-column covariance matrices, a two-stage algorithm must be solved to obtain their maximum likelihood estimators. A necessary and sufficient condition for the existence of maximum likelihood estimators is given and the question of their stability as solutions of the system of likelihood equations is addressed. In particular, the covariance matrix parameters and their maximum likelihood estimators are defined up to a positive multiplicative constant; only their direct product is uniquely defined. Using simulated data undertwovariance-covariancestructures that, otherwise, are indistinguishable by semivariance analysis, further specific aspects of the procedure are studied: (1) the convergence of the MLE algorithm is assessed; (2) the empirical bias of the direct product ofcovariance matrix estimators is calculated for various sample sizes; and (3) the consistency of the estimator is evaluated by its mean Euclidean distance from the parameter, as a function of the sample size. The adequacy of the matrix normal model, including the separability of the variance-covariance structure, is tested on multiple time series of dental medicine data; other applications to real doubly multivariate data are outlined.

Journal

Journal of Statistical Computation and SimulationTaylor & Francis

Published: Sep 1, 1999

Keywords: Matrix normal distribution; separability of variance-covariance structure; maximum likelihood estimation; two-stage algorithm; existence and stability of estimators; test of model adequacy

There are no references for this article.

Access the full text.

Sign up today, get DeepDyve free for 14 days.