插值以估算R

Question

每个ID都有两年数据，其中y2011和y2016分别代表2011年和2016年的数据。

id <- c(1:5)
y2011 <- c(100,200,150,121,20)
y2016 <- c(111,195,180,121,25)

dat <- data.frame(id, y2011, y2016)

现在，我想通过插值估计2012、2013、2014、2015年中期的近似数据。如何使用R做到这一点？我已经在R中尝试过approx函数。

approx(dat$y2011, dat$y2016)

但是没有得到适当的解决方案。

Answer 1

这是一种选择。实际上，您需要为每一行调用approx，所以我为此使用mapply。

middle_years <- t(mapply(function(a, b) approx(c(2011,2016), c(a,b), 2012:2015)$y, dat$y2011, dat$y2016))
colnames(middle_years) <- paste0("y", 2012:2015)

cbind(dat, middle_years)
#   id y2011 y2016 y2012 y2013 y2014 y2015
# 1  1   100   111 102.2 104.4 106.6 108.8
# 2  2   200   195 199.0 198.0 197.0 196.0
# 3  3   150   180 156.0 162.0 168.0 174.0
# 4  4   121   121 121.0 121.0 121.0 121.0
# 5  5    20    25  21.0  22.0  23.0  24.0