loglog_model_to_raw_plot

p <- 10^(seq(-4,0,0.2))
x <- 1:21
d <- cbind(x,p)
d <- d %>% as.data.frame()

loglogmodel <- lm(data = d, formula = log(p) ~ log(x))

ggplot(d, aes(x=x, y=p)) +
  geom_point() +
  geom_line(aes(x = exp(loglogmodel$model$`log(x)`), y = exp(loglogmodel$fitted.values)), color = "blue")