d <- read.table("inp_MRA.csv",header=TRUE,sep=",")
cor(d)
result<-lm(y~x1+x2+x3+x4,data=d)
summary(result)
coefficients(result)
predict(result)
residuals(result)
par(pty="s")
plot(
  d[,1],predict(result),
  log="xy",
  xlim=c(2000,200000),
  ylim=c(2000,200000)
)
abline(0,1)