Как добавить силовую линию тренда на график рассеяния в R?
Как я могу построить силовую линию тренда в виде y=ax^b на моем существующем графике рассеяния и есть ли уравнение и значение R-квадрата этой линии тренда отображаются на графике рассеяния?
В этом посте упоминается, как сделать любую другую линию тренда, кроме линий электропередач.
График, который у меня сейчас есть, с линиями линейной регрессии (lm функция), созданный с использованием следующего (упрощенного) кода:
# plotting
plot(y1, x1, pch=22)
points(y2, x2, pch=19) # plot data from second year on the same axis
text(y1, x1, labels=locationname)
legend("topleft", legend=c("Jun13-May14","Jun14-May15"), pch=c(22,19))
# linear regression fit
lm_mod1 <- lm(y1~x1)
abline(lm_mod1,col="red") # add regression line to the plot
lm_coef1 <- round(coef(lm_mod1), 3) # extract coefficients
mtext(bquote(y == .(lm_coef1[2])*x + .(lm_coef1[1])) # display equation
mtext(bquote(R^2==.(summary(lm_mod1314)$adj.r.squared)) # display R^2
В конце концов, моя цель - создать подобный график.