data.lm<-lm(data$y~data$x)
results<-summary(data.lm)
coefs<-coef(data.lm[1])
fit.lm<-cbind(data$x, coefs[2]*data$x+coefs[1])
plot(data$x, data$y, xlab=xtext, ylab=ytext)
lines(fit.lm)
print(summary(data.lm))