make a power curve with fitted data in ggplot2

power_curve_ggplot2.R

library(dplyr)

g  <- tibble(x = 1:10,
             y = rnorm(10, x^(0.2), 0.05))


library(ggplot2)

mod <- lm(log(y) ~ log(x), data = g)

pwr <- function(x) 
  exp(predict(mod, newdata = data.frame(x=x)))

ggplot(g,
       aes(x = x, y = y)) +
  geom_point() +
  stat_function(fun = pwr)