如何在R中使用approx函数对含NA的数据进行线性插值？

代码失效的核心原因 #

• 冒号运算符使用错误：R中:仅能生成整数序列，比如1:5会生成1、2、3、4、5，但你写的95.4968:101.491只会保留95到101的整数部分，丢失了小数精度；而103.856 : 44.7562因起始值大于结束值，会直接生成空向量，导致approx没有有效数据可计算。
• 未关联带缺失值的原始数据：你完全没用到带NA的id1向量，而是手动构造了无关的x、y参数，和填充id1中NA的需求完全脱节，自然无法完成缺失值填充。

正确的线性插值填充方法 #

针对带NA的id1向量，用approx的正确方式是基于原始数据的非NA值进行插值：

# 获取非NA值的位置和对应数值
non_na_pos <- which(!is.na(id1))
non_na_values <- id1[non_na_pos]

# 对所有位置进行插值，自动填充NA
filled_id1 <- approx(x = non_na_pos, y = non_na_values, xout = 1:length(id1), method = "linear")$y

# 验证填充结果
anyNA(filled_id1)

也可以用zoo包的na.approx函数，更简洁地处理缺失值插值：

library(zoo)
filled_id1 <- na.approx(id1)

内容的提问来源于stack exchange，提问作者Barta Boucherrougui