Likelihood-Based Inference with Separable Correlation Matrices

Abstract

This paper proposes methods for likelihood-based inference in multivariate linear regressions when the correlation matrix of the responses is separable; that is, it has a Kronecker product structure, but the variances are unrestricted. The methods are enabled by a block-coordinate ascent-like algorithm with closed-form updates that strictly increases the likelihood at every iteration until convergence. In the numerical experiments, the proposed algorithm is 300--2500 times faster than a general-purpose solver, making parametric bootstrap tests of correlation and covariance separability practical. Parameters are identifiable, and standard errors can therefore be obtained from the expected Fisher information, which can be computed efficiently using the Kronecker product structure. Simulations show that the proposed estimator has lower error than both separable covariance and unrestricted estimators when the model holds, and that bootstrap tests maintain nominal size where asymptotic tests fail. An application to dissolved oxygen data from the Mississippi River demonstrates that separable correlation captures location-specific variance patterns that separable covariance cannot.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…