d  <- read.table("inp_iris.txt",header=TRUE)
library(stats)
d.dist <- dist(d[,-5])        #Distance matrix
d.cmd <- cmdscale(d.dist,k=4) #MDS
x<-d.cmd
km <- kmeans(x,3) #k-means: number of cluster is 3
panel.text<-function(x, y, labels, cex, font, ...){
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(0,1,0,1) )
  txt <- c("cmdscale")
  cex<-2
  font<-"Helvetica"
  text(0.5, 0.5, txt, cex = cex)
} 
pairs(
  x[,1:4],
  diag.panel=NULL,
  text.panel=panel.text,
  pch=21,
  bg=c("red","green","blue")[unclass(km$cluster)]
  )