r - 如何绘制线性回归到双对数 R 图？

Question

我有以下数据:

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")

在我的图中添加回归线的最简单方法是什么？

Answer 1

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")