内插 NA 值

我有两组与时间无关的样本。我想合并它们并计算缺失值 在我不具备两者价值观的时代。简化示例：

A <- cbind(time=c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
           Avalue=c(1, 2, 3, 2, 1, 2, 3, 2, 1, 2))
B <- cbind(time=c(15, 30, 45, 60), Bvalue=c(100, 200, 300, 400))
C <- merge(A,B, all=TRUE)

   time Avalue Bvalue
1    10      1     NA
2    15     NA    100
3    20      2     NA
4    30      3    200
5    40      2     NA
6    45     NA    300
7    50      1     NA
8    60      2    400
9    70      3     NA
10   80      2     NA
11   90      1     NA
12  100      2     NA

通过假设每个样本之间的线性变化，可以计算缺失的 NA 值。 直观上很容易看出时间15和45处的A值应该是1.5。但对 B 的正确计算 例如在时间 20 将是

100 + (20 - 15) * (200 - 100) / (30 - 15)

等于 133.33333。 第一个括号是估计时间和最后一个可用样本之间的时间。 第二个括号是最近样本之间的差异。 第三个括号是最近样本之间的时间。

如何使用 R 计算 NA 值？

使用zoo包裹：

library(zoo)
Cz <- zoo(C)
index(Cz) <- Cz[,1]
Cz_approx <- na.approx(Cz)