Learn R Programming

TRES (version 0.1.0)

kroncov: The covariance estimation of a tensor random variable

Description

This function estimates the covariance of a tensor random variable. We assume the covariance of the tensor random variable has a seperable Kronecker covariance structure, i.e. \(\boldsymbol{\Sigma}=\boldsymbol{\Sigma}_{m}\otimes\cdots\otimes\boldsymbol{\Sigma}_{1}\). This algorithm is described in Manceur, A. M., & Dutilleul, P. (2013).

Usage

kroncov(Tn)

Arguments

Tn

A \(p_1\times\cdots p_m\times n\) data array, where \(n\) is the sample size.

Value

lambda

The normalizing constant.

S

A matrix lists with each element being the individual estimation of the seperable Kronecker covariance element \(\boldsymbol{\Sigma}_m,\ldots,\boldsymbol{\Sigma}_1\).

References

Manceur, A. M., & Dutilleul, P. (2013). Maximum likelihood estimation for the tensor normal distribution: Algorithm, minimum sample size, and empirical bias and dispersion. Journal of Computational and Applied Mathematics, 239, 37-49.