Anyone have an idea how to go about calculating the covariance between elements of a dirichlet distributed vector?