reslm <- lm(d$y~d$x)
summary(reslm)