I need to estimate the parameters of a set of bivariate gaussians from data. The data take the form of an n x n identification/confusion matrix where n is the the number of response alternatives. In its most basic form, the stimuli are formed from two dimensions (e.g., height, width), with two levels of each dimension, giving four response alternatives. It is also assumed that there is a decision bound for each dimension, giving a total of a mean vector and covariance matrix for each response alternative, and two decision bounds. Could this code be adapted to this problem, and if so, what would your recommendation be for proceeding. My thanks in advance!