这是我最初发布的问题的延续 http://r.789695.n4.nabble.com/subset-between-data-table-list-and-single-data-table-object-tp4673202.html。马修曾建议我在这里发表我的问题所以我现在正在这样做。
这是我在下面的输入：
library(data.table)
library(pracma) # for the interp1 function
tempbigdata1 <- data.table(c(14.80, 14.81, 14.82), c(7900, 7920, 7930), c("02437100", "02437100", "02437100"))
tempbigdata2 <- data.table(c(9.98, 9.99, 10.00), c(816, 819, 821), c("02446500", "02446500", "02446500"))
tempbigdata3 <- data.table(c(75.65, 75.66, 75.67), c(23600, 23700, 23800), c("02467000", "02467000", "02467000"))
tempsbigdata <- rbind(tempbigdata1, tempbigdata2, tempbigdata3)
setnames(tempsbigdata,c("y", "x", "site_no"))
setkey(tempsbigdata, site_no)
tempsbigdata
y x site_no
1: 14.80 7900 02437100
2: 14.81 7920 02437100
3: 14.82 7930 02437100
4: 9.98 816 02446500
5: 9.99 819 02446500
6: 10.00 821 02446500
7: 75.65 23600 02467000
8: 75.66 23700 02467000
9: 75.67 23800 02467000
aimsmall <- data.table(c("02437100", "02446500", "02467000"), c(3882.65, 819.82, 23742.37), c(1830.0, 382.0, 10400.0))
setnames(aimsmall,c("site_no", "mean", "p50"))
setkey(aimsmall, site_no)
aimsmall
site_no mean p50
1: 02437100 3882.65 1830
2: 02446500 819.82 382
3: 02467000 23742.37 10400
我使用此代码使用tempsbigdata$y的{{1}}值生成插值aimsmall$mean：
site_no这是函数 meanpre <- tempsbigdata[,if(aimsmall$mean > min(tempsbigdata$x){
interp1(tempsbigdata$x, tempsbigdata$y,
xi = aimsmall$mean, method ="linear")},by=site_no]
的输出，但不正确。
meanpre这是我想要的：
meanpre
site_no V1
1: 02437100 12.07599
2: 02437100 9.99410
3: 02437100 19.56813
4: 02446500 12.07599
5: 02446500 9.99410
6: 02446500 19.56813
7: 02467000 12.07599
8: 02467000 9.99410
9: 02467000 19.56813
有什么建议吗？谢谢。
更新1：
休，我过去使用过近似函数，对我的数据不准确;但是，meanpre
site_no V1
1: 02446500 9.99
2: 02467000 75.66
中的interp1函数是准确的。 pracma＆amp;中的mean和p50列aimsmall中的x值是放电值。 tempsbigdata中的y代表量具高度。我使用tempsbigdata函数来确定放电值的适当量具高度或interp1值，或y（以及mean）。
弗兰克，谢谢您的建议和建议的代码。这是您建议的代码的输出：
p50当我运行以下代码时，我得到以下结果：
tempsbigdata[aimsmall][,if(mean[1] > min(x)){interp1(tempsbigdata$x,tempsbigdata$y, xi = aimsmall$mean, method ="linear")},by=site_no]
site_no V1
1: 02446500 12.07599
2: 02446500 9.99410
3: 02446500 75.66424
4: 02467000 12.07599
5: 02467000 9.99410
6: 02467000 75.66424
有没有办法得到这个回报？谢谢。
interp1(tempsbigdata$x, tempsbigdata$y, xi = aimsmall$mean, method ="linear")
[1] 12.07599 9.99410 75.66424
更新2
弗兰克，谢谢你，我添加了代码，以便更容易在R中获取数据.Pracma是一个R包的数值方法例程，从GNU Octave [类似于MATLAB（R）]移植到R. site_no V1
1: 02446500 9.99
2: 02467000 75.66
函数是一维插值函数。
弗兰克，这很完美（关于“做东西”的R代码的最后评论）：
interp1