Existence and Uniqueness of the Kronecker Covariance MLE

03/12/2020
by   Mathias Drton, et al.
0

In matrix-valued datasets the sampled matrices often exhibit correlations among both their rows and their columns. A useful and parsimonious model of such dependence is the matrix normal model, in which the covariances among the elements of a random matrix are parameterized in terms of the Kronecker product of two covariance matrices, one representing row covariances and one representing column covariance. An appealing feature of such a matrix normal model is that the Kronecker covariance structure allows for standard likelihood inference even when only a very small number of data matrices is available. For instance, in some cases a likelihood ratio test of dependence may be performed with a sample size of one. However, more generally the sample size required to ensure boundedness of the matrix normal likelihood or the existence of a unique maximizer depends in a complicated way on the matrix dimensions. This motivates the study of how large a sample size is needed to ensure that maximum likelihood estimators exist, and exist uniquely with probability one. Our main result gives precise sample size thresholds in the paradigm where the number of rows and the number of columns of the data matrices differ by at most a factor of two. Our proof uses invariance properties that allow us to consider data matrices in canonical form, as obtained from the Kronecker canonical form for matrix pencils.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

11/07/2020

Maximum likelihood estimation for tensor normal models via castling transforms

In this paper, we study sample size thresholds for maximum likelihood es...
01/20/2018

Joint CLT for eigenvalue statistics from several dependent large dimensional sample covariance matrices with application

Let X_n=(x_ij) be a k × n data matrix with complex-valued, independent a...
06/27/2020

Generic canonical forms for perplectic and symplectic normal matrices

Let B be some invertible Hermitian or skew-Hermitian matrix. A matrix A ...
04/11/2020

Covariance Estimation for Matrix-valued Data

Covariance estimation for matrix-valued data has received an increasing ...
05/15/2022

Covariance Model with General Linear Structure and Divergent Parameters

For estimating the large covariance matrix with a limited sample size, w...
03/22/2015

Relaxed Leverage Sampling for Low-rank Matrix Completion

We consider the problem of exact recovery of any m× n matrix of rank ϱ f...
10/31/2017

Bayesian Learning of Random Graphs & Correlation Structure of Multivariate Data, with Distance between Graphs

We present a method for the simultaneous Bayesian learning of the correl...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.