BHC projections on PCA plot

makeFoIs.R

f <- "clusteringlopit/dunkley/hc_labels_observations_a.txt"
library("pRolocdata")
library("pRoloc")

data(dunkley2006)


makeFoIs <- function(fn, obj) {
    x <- readLines(fn)[-1] ## first is 'x'
    i <- grep("---", x)
    n <- length(i)
    res <- vector("list", length = n)
    for (j in 1:(n-1)) 
        res[[j]] <- FeaturesOfInterest(x[(i[j]+1):(i[j+1]-1)],
                                       paste("Cluster", j))
    res[[n]] <- FeaturesOfInterest(x[(i[n]+1):length(x)],
                                   paste("Cluster", n))
    FoICollection(res)
}

foic <- makeFoIs(f, dunkley2006)


for (i in seq_len(length(foic))) {
    .foi <- foi(foic)[[i]]
    plot2D(dunkley2006, main = description(.foi))
    highlightOnPlot(dunkley2006, .foi)
    scan(n = 1)
}