Multiple Heckman Selection Model

Abstract: We introduce a novel matrix-variate extension of the Heckman selection model to accommodate multiple outcomes, providing a flexible and natural generalization of classical selection models for matrix-valued data. By relying on the matrix normal distribution, the proposed model captures dependencies across both rows and columns while accounting for selection bias. An Expectation/Conditional Maximization (ECM) algorithm is developed, yielding closed-form updates for all model parameters. We investigate key theoretical properties, including the connection between sample selection models and the recently developed multivariate unified skew-normal (SUN) distribution. The performance of the proposed approach is assessed through simulation studies, and its practical utility is illustrated using two real datasets. The proposed method is implemented in the R package mvHeckman.