library(rggobi)
library(MASS)

fgl0 <- fgl[-10]
gg <- ggobi(fgl)
d <- display(gg$fgl, "2D Tour")
variables(d) <- list(X=0:8) # 0-based
glyph_colour(gg$fgl) <- fgl$type

dd <- dist(fgl0[-40, ])
fgl.iso <- isoMDS(dd)
eqscplot(fgl.iso$points, type = "n", xlab = "", ylab = "", axes = FALSE)
points(fgl.iso$points, col=unclass(fgl$type), pch=19)

fgl.iso3 <- isoMDS(dd, k = 3)
gg$iso3 <- fgl.iso3$points
glyph_colour(gg$iso3) <- fgl$type
d2 <- display(gg$iso3, "Rotation")

mkEdges <- function(ds)
{
    nm <- row.names(ds)
    A <- outer(nm, nm, function(x, y) x)
    a <- which(lower.tri(A), arr.ind=TRUE)
    cbind(nm[a[,1]], nm[a[,2]])
}
gg$dist <- as.vector(dist(fgl0))
edges(gg$dist) <- mkEdges(fgl0)

close(gg)