d   <- read.table("out_MRA.csv",skip=10,sep=",")
par(pty="s")
plot(
  d[,1],d[,2],
  log="xy",
  xlim=c(2000,200000),
  ylim=c(2000,200000),
  xlab="Observed values",
  ylab="Estimated values"
  )
abline(0,1)