[valse.git] / pkg / tests / testthat / helper-clustering.R

# Compute the sum of (normalized) sum of squares of closest distances to a medoid.
computeDistortion <- function(series, medoids)
{
	n <- ncol(series)
	L <- nrow(series)
	distortion <- 0.
	for (i in seq_len(n))
		distortion <- distortion + min( colSums( sweep(medoids,1,series[,i],'-')^2 ) / L )

	sqrt( distortion / n )
}
Commit	Line	Data
5ce95f26 BA	1	# Compute the sum of (normalized) sum of squares of closest distances to a medoid.
	2	computeDistortion <- function(series, medoids)
	3	{
	4	n <- ncol(series)
	5	L <- nrow(series)
	6	distortion <- 0.
	7	for (i in seq_len(n))
	8	distortion <- distortion + min( colSums( sweep(medoids,1,series[,i],'-')^2 ) / L )
	9
	10	sqrt( distortion / n )
	11	}