r - 如何在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？我需要获得n = floor(abs(y1 - y2))相邻y1和之间的点y2。

标签： rinterpolationpurrr

Answer 1

您可以使用purrr::pmap组合 with在和dplyr::lag之间进行插值。 对于两点之间的线性插值，将间隔削减或并不重要，因此仍然有效：x[i-1]x[i]
xyapproxfun(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