如何绘制线性回归到双对数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")