getSWR: Calculate modes
In EMMIXskew: The EM Algorithm and Skew Mixture Distribution

Calculate the distances.

getSWR(dat,g,sigma,clust, tau)
intradist(dat,g,sigma, clust, tau) 
interdist(dat,g,sigma, clust, tau) 
mahalonobis(p, g, mu, sigma)

`dat`	The dataset, an n by p numeric matrix, where n is number of observations and p the dimension of data.
`p`	The dimension of the data
`g`	The number of components to be fit
`mu`	A numeric matrix with each column corresponding to the mean
`sigma`	An array of dimension (p,p,g) with first two dimensions corresponding covariance matrix of each component
`clust`	A vector of integers specifying the initial partitions of the data; the default is NULL.
`tau`	An n by g matrix of posterior probability for each data point