d   <- read.table("out_MDS_iris.csv",skip=10,sep=",")
dat <- d[,2:ncol(d)]
names(dat) <- c("MDS","MDS","MDS","MDS")
pairs(
  dat,
  diag.panel=NULL,
  pch=21
  )