library(lattice)

library(mlbench)
data(BreastCancer)
data <- na.omit(BreastCancer)
n <- dim(data)[1]
p <- dim(data)[2]
X <- matrix(unlist(data[,2:10]),n,p)+runif(n=dim(X)[1])
name <- data[,11]


c.obj <- cmdscale(dist(X),k=3)
c.obj <- sammon(dist(X),k=3)$points
c.obj <- isoMDS(dist(X),k=3)$points

y <- c(c.obj[,2],c.obj[,3],c.obj[,3])
x <- c(c.obj[,1],c.obj[,1],c.obj[,2])
panl <- rep(c("1 vs 2","1 vs 3","2 vs 3"),each=n)
xyplot(y~x|panl,groups=rep(name,times=3),pch=16,col=1:6)