将回归线添加到对数 x 刻度

Question

我在非常宽的 x 范围内有一些数据，我正在使用 R 中的对数 x 轴绘制这些数据。当我使用 MS Excel 绘制数据并添加线性趋势线时，结果输出是我希望看到的给定我的数据集和预期结果（下图）。

我想使用 R 来实现这一点。我可以绘制数据并运行回归模型，但 R 似乎不想添加回归线/曲线。

这是我的数据：

data <-
  data.frame(
    x = c(
      6.37E-05,
      0.000584271,
      0.001044085,
      0.00011272,
      3.95E-05,
      6.05E-05,
      0.000277963,
      0.001800305,
      3.70E-05,
      0.002588335,
      0.000710123,
      0.000148309,
      0.001069184,
      0.001877917,
      0.01086776,
      0.043686462,
      0.000276426
    ),
    y = c(
      -0.3375,
      -0.35,
      -0.35,
      -0.35,
      -0.3,
      -0.35,
      -0.3625,
      -0.375,
      -0.3,
      -0.3625,
      -0.3,
      -0.3375,
      -0.29,
      -0.34,
      -0.4,
      -0.48,
      -0.3375
    )
  )

这是迄今为止我一直在尝试的代码。

attach(data)
    plot(data$x,data$y,log="x",xlab="independent variable",
     ylab="dependent variable")
model<-lm(data$y~data$x)
abline(model) # THIS STEP DOESN'T WORK

我有三个查询：

1)。如何在基本 R 中的 y~logx 图上显示回归线？

2)。我该如何重复 1）。使用ggplot？

3)。当回归线/曲线开始“下降”时，我如何找到 (log)x 值？

谢谢。

Answer 1

你能让这个例子重现吗？ 我认为这可能会奏效。

dataframe %>%
ggplot(aes(x = log(x), y = y) +
 geom_point() +
 geom_smooth(method = "loess") +
 xlab("independent variable")

这使用 dplyr 和 ggplot2。

Answer 2

ggplot(data, aes(x, y)) + geom_point() + stat_smooth(method = "lm")

ggplot(data, aes(x, y)) + geom_point() + stat_smooth(method = "lm") +
  coord_trans(x = "log10")

如您所见，回归线没有“开始‘下降’”的点。

---

更新

如果您在绘图之前对 x 值进行对数转换，您最终会得到如下图：

如果你真正感兴趣的是变化点分析，你可以适应它。

data.fit <- nls(100*y ~ pmin(plat, intx + slx*log(x)), 
                  data = data,
                  start = list(plat = -32, intx = -60, slx = -4), 
                  control = nls.control(warnOnly = T))
summary(data.fit)
coef(data.fit)