Use the numpy package. numpy.mean and numpy.cov will give you a gaussian parameter estimate. Assuming you have 13 attributes and N is the number of observations, you need to set rowvar=0 when calling numpy.cov for your matrix N x 13 (or pass the transpose of your matrix as an argument to the function).
If your data is in a numpy data array:
mean = np.mean(data, axis=0) cov = np.cov(data, rowvar=0)
source share