d <- read.table("inp_MRA.csv",skip=2,sep=",")
names(d) <- c("y","x1","x2","x3","x4")
pairs(
  d,
  diag.panel=NULL,
  pch=21
  )