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