Use `stat_function()` to plot a single function

Published

November 3, 2023

set.seed(45345)
fake_dat <- tibble(
  x = seq(-5, 5, length.out = 1000),
  y = x^3 + rnorm(1000, mean = 0, sd = 15)
)

point_plot <- fake_dat |> 
  ggplot(aes(x, y)) +
  geom_point(col = "#0072B2", alpha = 0.15, size = 3)

point_plot +
  stat_function(
    fun = \(x) x^3, # Put formula from model into a function
    col = "#009E73",
    linewidth = 2.5,
    lineend = 'round'
  ) +
  labs(title = 'Any function can be plotted with stat_function()')

Enjoyed this code snippet?