Regression for matrix-valued data via Kronecker products factorization