我有以下数据框
ref_dat
k Intensity Slope
1 0.021467214 33 16
2 0.012444759 33 8
3 0.006079156 33 4
4 0.003792025 33 2
5 0.022765086 47 16
6 0.017335375 47 8
7 0.007279113 47 4
8 0.003396341 47 2
9 0.034101951 75 16
10 0.023045525 75 8
11 0.009163758 75 4
12 0.003253822 75 2
13 0.039832629 110 16
14 0.029826826 110 8
15 0.010529787 110 4
16 0.002945170 110 2
17 0.042140630 155 16
18 0.032338534 155 8
19 0.013758171 155 4
20 0.002605399 155 2
我正在尝试使用interp中的akima函数获取给定斜率和强度集合的k的插值（和外推）值，所以我正在执行以下操作。
我的新坡度
s=10
我的新强度
i=c(seq(0,50,5),rev(seq(0,50,5)))
并计算新的k值
k_new=interp(x=ref_dat$Intensity,y=ref_dat$Slope,z=ref_dat$k,xo=i,yo=s,extrap = FALSE)$z
，结果是
k_new
[,1]
[1,] NA
[2,] NA
[3,] NA
[4,] NA
[5,] NA
[6,] NA
[7,] NA
[8,] 0.01539903
[9,] 0.01714568
[10,] 0.01850739
[11,] 0.01930460
[12,] 0.01930460
[13,] NA
[14,] NA
[15,] NA
[16,] NA
[17,] NA
[18,] NA
[19,] NA
[20,] NA
[21,] NA
[22,] NA
转向推断
k_new_ext=interp(x=ref_dat$Intensity,y=ref_dat$Slope,z=ref_dat$k,xo=i,yo=s,extrap = TRUE)$z
和reult
k_new_ext
[,1]
[1,] NA
[2,] NA
[3,] NA
[4,] NA
[5,] NA
[6,] NA
[7,] NA
[8,] 0.01539903
[9,] 0.01714568
[10,] 0.01850739
[11,] 0.01930460
[12,] 0.01930460
[13,] NA
[14,] NA
[15,] NA
[16,] NA
[17,] NA
[18,] NA
[19,] NA
[20,] NA
[21,] NA
[22,] NA
稍微探讨一下结果
cbind(i,k_new,k_new_ext)
i
[1,] 0 NA NA
[2,] 5 NA NA
[3,] 10 NA NA
[4,] 15 NA NA
[5,] 20 NA NA
[6,] 25 NA NA
[7,] 30 NA NA
[8,] 35 0.01539903 0.01539903
[9,] 40 0.01714568 0.01714568
[10,] 45 0.01850739 0.01850739
[11,] 50 0.01930460 0.01930460
[12,] 50 0.01930460 0.01930460
[13,] 45 NA NA
[14,] 40 NA NA
[15,] 35 NA NA
[16,] 30 NA NA
[17,] 25 NA NA
[18,] 20 NA NA
[19,] 15 NA NA
[20,] 10 NA NA
[21,] 5 NA NA
[22,] 0 NA NA
为什么我得到NA值，即使新的强度值在ref_dat的值范围内？让extrap=TRUE也无法解决我的问题。