Formally it is nothing but distance multivariance applied to the Monte Carlo emprical transform of the data. Hence its values vary for repeated runs.
copula.multivariance(x, vec = 1:ncol(x), type = "total", ...)either a data matrix or a list of doubly centered distance matrices
if x is a matrix, then this indicates which columns are treated together as one sample; if x is a list, these are the indexes for which the multivariance is calculated. The default is all columns and all indexes, respectively.
default: "total.lower.upper", for details and other options see below
these are passed to cdms (which is only invoked if x is a matrix)
For the theoretic background see the reference [5] given on the main help page of this package: multivariance-package.