R语言regularize.values警告问题：保留重复v2值的解决方法

问题原因 #

你遇到的警告是因为splinefun和approx要求输入的自变量（这里是v2）必须是唯一值。当存在重复的v2时，函数会自动合并重复项并抛出警告，但原数据的重复值其实并没有被修改——只是插值过程中用到了去重后的值。

解决方案 #

我们可以手动提取v2和对应y_s的唯一组合，用这个唯一组合去做插值，同时保留原始的v2变量供后续使用。这样既消除了警告，又不影响后续对原重复值的操作。

修改后的完整代码 #

a = 1.1
b = 1.3

C = c(0, 0.3, 1.5, 3.5, 19.5)
v2 = c(0.00, 0.56, 1.00, 1.00, 1.00) 

# 定义原函数
fun <- function(C, a, b) {
  1 / (1 + (C / a)^b) 
}

x_d = seq(0, 10, by = 0.1)
y_s <- fun(v2, a, b)  

# 提取v2和y_s的唯一组合，避免插值函数警告
unique_pairs <- unique(data.frame(x = v2, y = y_s))

# 使用唯一值做样条插值
y_d_spline <- splinefun(unique_pairs$x, unique_pairs$y)(x_d)
# 使用唯一值做线性近似（修正原代码取$x的错误，近似结果应为$y）
y_d_approx <- approx(unique_pairs$x, unique_pairs$y, xout = x_d)$y

# 后续查找操作基于修正后的插值结果
which(y_d_approx <= 0.5)
which(y_d_approx <= 0.5)[1]
x_d[which(y_d_approx <= 0.5)[1]]

关键说明 #

• unique(data.frame(x = v2, y = y_s))会自动保留v2的唯一值，以及对应的y_s（重复的v2计算出的y_s完全相同，取唯一不会丢失有效信息）。
• 原始的v2变量完全未被修改，你后续可以正常使用它的所有重复值。
• 修正了原代码中approx(...)取$x的问题：approx返回的$x是输入的xout（即x_d），实际的近似结果存在$y中，这个调整能让你得到正确的插值结果。

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