如何在R中用y插值向量

Question

我有一个带坐标的向量

x <- c(1, 3, 5, 7)
y <- c(2.6, 7.3, 6.8, 1.9)

我需要同时通过x 和y 对向量 (x, y) 进行插值。 x我可以做

approxfun(x, y)(1:7)

但是如何通过y 进行插值？我需要在相邻的y1 和y2 之间获取n = floor(abs(y1 - y2)) 点。

Answer 1

您可以使用purrr::pmap 结合dplyr::lag 在x[i-1] 和x[i] 之间进行插值。
对于两点之间的线性插值，将间隔削减x 或y 并不重要，这样approxfun(x, y) 仍然有效：

library(dplyr)
library(purrr)

x <- c(1, 3, 5, 7)
y <- c(2.6, 7.3, 6.8, 1.9)


df <- data.frame(x=x,y=y)

df %>% mutate(xprev = lag(x),yprev=lag(y)) %>% 
  purrr::pmap(~with(list(...),{
    if (is.na(xprev)) {
      data.frame(x=x,y=y)
    } else {
      n <-  floor(abs(yprev-y)) + 2
      if (n>0) {
        xstep <- (x-xprev) / n  
        xsteps <- seq(xprev+xstep,x ,by = xstep )       
        data.frame(x=xsteps, y=approxfun(c(xprev,x),c(yprev,y))(xsteps))
      }
    }
  })) %>% bind_rows 
#>     x     y
#> 1 1.0 2.600
#> 2 1.5 3.775
#> 3 2.0 4.950
#> 4 2.5 6.125
#> 5 3.0 7.300
#> 6 5.5 5.575
#> 7 6.0 4.350
#> 8 6.5 3.125
#> 9 7.0 1.900