In matrix-valued datasets the sampled matrices often exhibit correlations among both their rows and columns. A useful parsimonious model of such dependence is matrix normal model, in which covariances elements a random are parameterized terms Kronecker product two covariance matrices, one representing row column covariance. An appealing feature that structure allows for standard likelihood infe...