与区间匹配并提取两个矩阵 R 之间的值

我的列表中有 n 个矩阵，还有一个附加矩阵，其中包含我想要在矩阵列表中查找的值。

为了获取矩阵列表，我使用以下代码：

setwd("C:\\~\\Documents\\R") 


import.multiple.txt.files<-function(pattern=".txt",header=T)
{
list.1<-list.files(pattern=".txt")
list.2<-list()
for (i in 1:length(list.1))
{
list.2[[i]]<-read.delim(list.1[i])
}
names(list.2)<-list.1
list.2

}


txt.import.matrix<-cbind(txt.import)

我的列表如下所示：（我仅显示 n=2 的示例）。每个数组中的行数不同（这里我只取 5 行和 6 行来简化，但我的真实数据有超过 500 行）。

txt.import.matrix[1]

    [[1]]
     X.     RT.     Area.  m.z.      
1     1     1.01   2820.1  358.9777  
2     2     1.03   9571.8  368.4238  
3     3     2.03   6674.0  284.3294  
4     4     2.03   5856.3  922.0094  
5     5     3.03   27814.6 261.1299  


txt.import.matrix[2]

    [[2]]
     X.     RT.     Area.  m.z.      
1     1     1.01    7820.1 358.9777  
2     2     1.06    8271.8 368.4238  
3     3     2.03   12674.0 284.3294  
4     4     2.03    5856.6 922.0096  
5     5     2.03   17814.6 261.1299
6     6     3.65    5546.5 528.6475  

我想在矩阵列表中找到另一个值数组。该数组是通过将列表中的所有数组合并到一个数组中并删除重复项而获得的。

reduced.list.pre.filtering

     RT.   m.z.
1    1.01  358.9777
2    1.07  368.4238
3    2.05  284.3295
4    2.03  922.0092
5    3.03  261.1299
6    3.56  869.4558

我想获得一个新的矩阵，其中写有Area.匹配结果RT. ± 0.02 and m.z. ± 0.0002对于列表中的所有矩阵。输出可能是这样的。

     RT.   m.z.        Area.[1]      Area.[2]
1    1.01  358.9777    2820.1        7820.1
2    1.07  368.4238                  8271.8      
3    2.05  284.3295    6674.0        12674.0
4    2.03  922.0092    5856.3             
5    3.03  261.1299    27814.6            
6    3.65  528.6475    

我只知道如何仅匹配一个数组中的一个精确值。这里的难点是在数组列表中查找值，并且需要找到该值±一个区间。如果您有任何建议，我将非常感激。

Using non-equi从 data.table v1.9.7 的当前开发版本加入（请参阅安装说明 https://github.com/Rdatatable/data.table/wiki/Installation)，它允许为连接提供非相等条件：

require(data.table) # v1.9.7
names(ll) = c("Area1", "Area2")
A = rbindlist(lapply(ll, as.data.table), idcol = "id")           ## (1)

B = as.data.table(mat)
B[, c("RT.minus", "RT.plus") := .(RT.-0.02, RT.+0.02)]
B[, c("m.z.minus", "m.z.plus") := .(m.z.-0.0002, m.z.+0.0002)]   ## (2)

ans = A[B, .(id, X., RT. = i.RT., m.z. = i.m.z., Area.), 
           on = .(RT. >= RT.minus, RT. <= RT.plus, 
                  m.z. >= m.z.minus, m.z. <= m.z.plus)]          ## (3)

dcast(ans, RT. + m.z. ~ id)                                      ## (4)
# or dcast(ans, RT. + m.z. ~ id, fill = 0)
#     RT.     m.z.   Area1   Area2
# 1: 1.01 358.9777  2820.1  7820.1
# 2: 1.07 368.4238      NA  8271.8
# 3: 2.03 922.0092  5856.3      NA
# 4: 2.05 284.3295  6674.0 12674.0
# 5: 3.03 261.1299 27814.6      NA

[1] 命名矩阵列表（称为ll此处）并将它们每个转换为数据表 using lapply()，并使用按行绑定它们rbindlist，并将名称添加为额外列 (idcol）。叫它A.

[2] 转换第二个矩阵（称为mat在这里）到数据表以及。添加与您要搜索的范围/间隔相对应的其他列（因为on=正如我们将在下一步中看到的，参数无法处理表达式yet）。叫它B.

[3] 执行条件连接/子集。对于中的每一行B，找到匹配的行A对应于提供给的条件on=参数，并提取列id, X., R.T. and m.z.对于那些匹配的行索引。

[4] 最好将其保留为 [3]。但如果您希望它如答案中所示，我们可以将其重新塑造为宽格式。fill = 0将取代NA结果中的 s 为0.