function(explan, depend, klepm)
{
	d.stats <- 0
	K <- ncol(explan)
	dmat <- cbind(explan, depend)
	dmat <- as.matrix(na.omit(dmat))
	tmat <- as.matrix(cov.wt(dmat)$cov)
	hmat <- tmat
	hmat[K + 1, K + 1] <- klepm * tmat[K + 1, K + 1]
	ihmat <- solve(hmat)
	d.stats <- NULL
	for(i in 1:K) {
		dstat <- 1/(ihmat[i, i] * tmat[i, i])
		d.stats <- c(d.stats, dstat)
	}
	names(d.stats) <- dimnames(explan)[[2]]
	list(d.stats = d.stats, mbound = klepm)
}
