如何在我现有的散点图上以y=ax^b
的形式绘制功率趋势线,并在散点图上显示该趋势线的等式和R平方值?
此post提到了除电源线以外的其他每条趋势线。
我目前拥有的plot是线性回归线(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
我最终的目标是拥有类似this的图表。