R中点的线性插值

Question

这似乎是一个非常简单的问题，但这里有：

我有一个数据框：

test_df <- data.frame(x1 = c(277422033,24118536.4,2096819.0,
                               182293.4,15905,1330,105,16,1),
                      x2 = c(2.496e-3,2.495e-2,2.496e-1,
                               2.496e0,2.47e1,2.48e2,2.456e3,
                               3.7978e4,3.781e5))

我想对其进行线性插值以增加点数。 这些变量在对数刻度上线性相关，即

plot(log10(test_df[,1]),log10(test_df[,2]))

因此，我的问题是，如何线性内插这些值以增加值的数量？

这是我尝试使用线性模型（而不是近似函数）的尝试：

我将线性模型定义为：

test.lm <- lm(log10(x1) ~ log10(x2), data = test_df)

然后为新点定义一个新变量：

ss <- seq(min(test_df$x2),max(test_df$x2),length.out = 100) # new x1

然后预测新值并绘制点

newY <- predict(test.lm, newdata = data.frame(x2 = ss)) # interpolated values

test_df2 <- data.frame(x1 = 10^newY,
                       x2 = ss)

points(newY,log10(ss),col = "red")

这按我预期的那样工作，即最终的图形符合我的预期。

我想增加test_df2中的点数，这可以通过增加length.out来完成。

ss <- seq(min(test_df$x2),max(test_df$x2),length.out = 10000000)

但这使我的机器上的运行时间很长，以至于我必须重新启动R。

有没有一种方法可以线性插值均匀分布的点数，从而扩展ss中指定的全部点数？

Answer 1

只需使用

ss <- 10^seq(log10(min(test_df$x2)),log10(max(test_df$x2)),length.out = 1000)

使您的新数据在对数刻度上均匀分布。