library(ggplot2) library(plotly) library(dplyr) data(airquality) data(cars) ggplotter <- ggplot(airquality, aes(x = Temp, y = Ozone)) + geom_point() + geom_smooth(method = "lm") + coord_cartesian(ylim = c(0,100)) ggplotter2 <- ggplot(cars, aes(x = speed, y = dist)) + geom_point() + geom_smooth(method = "lm") + coord_cartesian(ylim = c(0,100)) airquality_na <- na.omit(airquality) line <- lm(Ozone ~ Temp, data = airquality_na) line2 <- lm(dist ~ speed, data = cars) x = airquality_na$Temp; y = airquality_na$Ozone plotter = plot_ly(data = airquality_na,x= ~Temp,y = ~Ozone, type = "scatter" ,mode = "markers", width = 400, height = 400,name = "Data") %>% add_lines(x = x, y = fitted(line), name = "Regression")