In many practical situations we would like to estimate the covariance matrix of a set of variables from an insufficient amount of data. More specifically, if we have a set of N independent, identically distributed measurements of an M dimensional random vector the maximum likelihood estimate is the sample covariance matrix. Here we consider the case where N
