51 lines
No EOL
1.2 KiB
R
51 lines
No EOL
1.2 KiB
R
cluster_own <- function(Z, nclust) {
|
|
true <- logical(1)
|
|
false <- logical(0)
|
|
maxclust <- nclust
|
|
# % Start of algorithm
|
|
m <- size(Z, 1) + 1
|
|
T <- zeros(m, 1)
|
|
# % maximum number of clusters based on inconsistency
|
|
if (m <= maxclust) {
|
|
T = t((1:m))
|
|
} else if (maxclust == 1) {
|
|
T <- ones(m, 1)
|
|
} else {
|
|
clsnum <- 1
|
|
for (k in (m - maxclust + 1):(m - 1)) {
|
|
i = Z(k, 1) # left tree
|
|
if (i <= m) { # original node, no leafs
|
|
T(i) = clsnum
|
|
clsnum = clsnum + 1
|
|
} else if (i < (2 * m - maxclust + 1)) { # created before cutoff, search down the tree
|
|
T <- clusternum(Z, T, i - m, clsnum)
|
|
clsnum <- clsnum + 1
|
|
}
|
|
i <- Z(k, 2) # right tree
|
|
if (i <= m) { # original node, no leafs
|
|
T[i] <- clsnum
|
|
clsnum <- clsnum + 1
|
|
} else if (i < (2 * m - maxclust + 1)) { # created before cutoff, search down the tree
|
|
T <- clusternum(Z, T, i - m, clsnum)
|
|
clsnum <- clsnum + 1
|
|
}
|
|
}
|
|
}
|
|
return(T)
|
|
}
|
|
|
|
clusternum <- function(X, T, k, c) {
|
|
m <- size(X, 1) + 1
|
|
while (!isempty(k)) {
|
|
# Get the children of nodes at this level
|
|
children <- X[k, 1:2]
|
|
children <- children
|
|
|
|
# Assign this node number to leaf children
|
|
t <- (children <= m)
|
|
T[children[t]] <- c
|
|
# Move to next level
|
|
k <- children(!t) - m
|
|
}
|
|
return(T)
|
|
} |