Joint Mean And Covariance Modeling Of Matrix-Variate Data