mydata = data.frame(
  fitness.minutes = c(5, 10, 15, 20, 25),
  cholesterol = c(210, 182, 166, 163, 145)
)

fit <- lm(cholesterol ~ fitness.minutes, data=mydata)
plot(cholesterol~fitness.minutes, data=mydata, pch = 16, cex = 1.3, col = "blue", main = "Cholesterol against daily activity", xlab = "Minutes of activities", ylab = "Cholesterol (mg/dl)")
fit
abline(fit)