d   <- read.table("inp_MRA.csv",header=TRUE,sep=",")
mm  <- ncol(d)
yd  <- d[,1]
xd  <- cbind(1,d[,2:mm])
xd  <- as.matrix(xd)
xy  <- t(xd) %*% yd
xx  <- t(xd) %*% xd
bb  <- solve(xx,xy)
ed  <- xd %*% bb
rro <- cor(d)
rri <- solve(rro)
mcc <- sqrt(1-1/rri[1,1])
bb
rro
rri
mcc
par(pty="s")
plot(
  yd,ed,
  log="xy",
  xlim=c(2000,200000),
  ylim=c(2000,200000),
  xlab="Observed values",
  ylab="Estimated values"
  )
abline(0,1)