d   <- read.table("inp_iris.csv",skip=2,sep=",")
dat <- d[,1:4]
names(dat) <- c("x1","x2","x3","x4")
pairs(
  dat,
  diag.panel=NULL,
  pch=21
  )