inst/unit_tests/runit.storedistance.R

In Rclusterpp: Linkable C++ clustering

compare.hclust <- function(h1, h2) {
	checkEquals(h1$merge, h2$merge, msg="Agglomerations don't match")
	checkEquals(h1$height, h2$height, msg="Agglomeration heights are not equal")
	checkEquals(h1$labels, h2$labels, msg="Cluster labels do not match")
	checkEquals(h1$order, h2$order, msg="Cluster orders do not match")
}

test.storedistance.average.euclidean <- function() {
	h <- hclust(dist(USArrests, method="euclidean"), method="average")
	r <- Rclusterpp.hclust(dist(USArrests, method="euclidean"), method="average")
	compare.hclust(h, r)
}

test.storedistance.average.manhattan <- function() {
	h <- hclust(dist(USArrests, method="manhattan"), method="average")
	r <- Rclusterpp.hclust(dist(USArrests, method="manhattan"), method="average")
	compare.hclust(h, r)
}

test.storedistance.single.euclidean <- function() {
	h <- hclust(dist(USArrests, method="euclidean"), method="single")
	r <- Rclusterpp.hclust(dist(USArrests, method="euclidean"), method="single")
	compare.hclust(h, r)
}

test.storedistance.complete.euclidean <- function() {
	h <- hclust(dist(USArrests, method="euclidean"), method="complete")
	r <- Rclusterpp.hclust(dist(USArrests, method="euclidean"), method="complete")
	compare.hclust(h, r)
}