Shape inference based on multivariate normal matrix distribution