我有以下数据:
someFactor = 500
x = c(1:250)
y = x^-.25 * someFactor
我以双对数图显示:
plot(x, y, log="xy")
现在我使用线性模型“找出”数据的斜率:
model = lm(log(y) ~ log(x))
model
这使:
Call:
lm(formula = log(y) ~ log(x))
Coefficients:
(Intercept) log(x)
6.215 -0.250
现在我想将线性回归绘制为红线,但 abline 不起作用:
abline(model, col="red")
将回归线添加到图中的最简单方法是什么?
lines(log(x), exp(predict(model, newdata=list(x=log(x)))) ,col="red")
在对数刻度上绘制的 x 值范围和用作自变量的 log(x) 值范围实际上有很大不同。这将为您提供全方位:
lines(x, exp(predict(model, newdata=list(x=x))) ,col="red")
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)