Here is a linear model:

x <- 1:10
y<-rnorm(10)
plot(x,y, pch=16,col="blue")

Fitting a linear maodel:

fit <- lm(y~x)