R 在双对数轴绘图机中绘制幂函数线

我在将幂函数线绘制到对数刻度上具有两个轴的绘图时遇到问题。我需要将轴保持在原木刻度上。我不想对数据进行日志转换。这就是我在plot函数上使用log=xy参数的原因，如下所示

但是，我还不能将幂函数添加到这个图中。到目前为止，我尝试了不同的方法，包括下面的nls，但似乎没有任何效果

谁能帮我一下吗？ 非常感谢。 玛拉

---

在这里获得帮助的最佳方法是提供一个完全可复制的示例，您给了我们您的代码，但您没有提供slpca1的数据，并尽可能简单地保留该示例，不需要轴标签、字体等，除非这是问题的基本部分

所以，我创建了一些假数据，拟合你的非线性回归，并绘制了拟合结果。希望这将帮助您开始，并且您可以根据需要进行修改

# create some fake data so I can run your code
CA <- runif(20, 1, 20)
Slope <- 0.1*CA^-0.04 + rnorm(20, 0, 0.001)
slpca1 <- data.frame(CA, Slope)
rm(CA, Slope)

# fit a nonlinear regression
nls4 = nls(Slope ~ i*CA^-z, start=list(i=0.1,z=-0.04), data=slpca1)
# plot the observed values
plot(Slope~CA, data=slpca1, log="xy")
# plot the predicted line
ord <- order(slpca1$CA)
lines(slpca1$CA[ord], predict(nls4)[ord])

---

您似乎对R中的一切工作方式相当困惑。为什么你画了你的横线，然后你读了？阴谋和PAR学会如何在一个图表上放置多个情节？

# create some fake data so I can run your code
CA <- runif(20, 1, 20)
Slope <- 0.1*CA^-0.04 + rnorm(20, 0, 0.001)
slpca1 <- data.frame(CA, Slope)
rm(CA, Slope)

# fit a nonlinear regression
nls4 = nls(Slope ~ i*CA^-z, start=list(i=0.1,z=-0.04), data=slpca1)
# plot the observed values
plot(Slope~CA, data=slpca1, log="xy")
# plot the predicted line
ord <- order(slpca1$CA)
lines(slpca1$CA[ord], predict(nls4)[ord])