Nothing
DuplicateClusterRemove <- function(dpobj){
cp <- dpobj$clusterParameters
dup_array <- matrix(cp[[1]], ncol=ncol(dpobj$data), byrow = TRUE)
dup <- duplicated(dup_array)
if (all(!dup)){
return(dpobj)
}
dupLabels <- seq_len(dpobj$numberClusters)
inds <- match(data.frame(t(dup_array)),
data.frame(t(unique(dup_array))))
oldLabs <- dpobj$clusterLabels
newLabs <- oldLabs
for(i in seq_along(dupLabels)){
if(dupLabels[i] == inds[i]){
next
}
else{
newLabs[which(oldLabs == dupLabels[i])] <- inds[i]
}
}
dpobj$clusterLabels <- newLabs
newCP <- cp
newCP[[1]] <- cp[[1]][,,which(!dup), drop=FALSE]
newCP[[2]] <- cp[[2]][,,which(!dup), drop=FALSE]
newPointsPerCluster <- vapply(1:max(newLabs),
function(x) sum(newLabs==x),
numeric(1))
while (any(newPointsPerCluster==0)){
newLabs[newLabs > which(newPointsPerCluster == 0)] = newLabs[newLabs > which(newPointsPerCluster == 0)] - 1
newPointsPerCluster <- sapply(1:max(newLabs),
function(x) sum(newLabs == x),
numeric(1))
#print(oldLabs)
#print(newLabs)
print(c(newCP[[1]]))
}
dpobj$clusterLabels <- newLabs
dpobj$clusterParameters <- newCP
dpobj$numberClusters <- length(unique(newLabs))
dpobj$pointsPerCluster <- vapply(1:max(newLabs),
function(x) sum(newLabs == x),
numeric(1))
return(dpobj)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.