Соедините два data.tables одной числовой переменной с динамическим допуском - PullRequest
Новая
       37

Соедините два data.tables одной числовой переменной с динамическим допуском

0 голосов
/ 09 июля 2019

У меня есть два data.tables, к которым я хотел бы присоединиться через одну числовую переменную (двойная точность).Тем не менее, числовая переменная имеет недостатки с неопределенностью.Поэтому я должен учесть определенный допуск, который отличается в зависимости от переменной.

В приведенном ниже примере «mz» - это переменная, с помощью которой я хотел бы объединить DT1 и DT2.Допуск, рассчитанный из переменной iso_mz: iso_mz * 5e-6.

DT1 <- data.table(mz = c(433.231512451172, 451.091953822545, 454.347605202415, 490.167234693255, 518.225894504123), 
Var1 = c(433.231018066406, 451.091430664062, 454.347015380859, 490.166381835938, 518.22509765625), 
Var2 = c(433.232147216797, 451.092559814453, 454.34814453125, 490.168273925781, 518.2265625))


DT2 <- data.table(iso_mz = c(451.0900, 490.1651, 518.2281, 433.2335), 
comp = c("m1", "m2", "m3", "m4"))

Если бы мне не пришлось использовать допуск, я бы использовал функциональность on =. () Пакета data.table.Я пытался адаптировать код из Соединение фреймов данных по измеренным значениям с диапазоном ошибок , но по какой-то причине я не смог запустить, ..

Желаемый вывод для моего исследования будет: 

Output <- data.table(
iso_mz = c(433.2335, 451.0900, 490.1651, 518.2281), 
comp = c("m4", "m1", "m2", "m3"),
mz = c(433.231512451172, 451.091953822545, 490.167234693255, 518.225894504123), 
Var1 = c(433.231018066406, 451.091430664062, 490.166381835938, 518.22509765625), 
Var2 = c(433.232147216797, 451.092559814453, 490.168273925781, 518.2265625))

Заранее спасибо!

1 Ответ

1 голос
/ 09 июля 2019

Вот подход, использующий foverlaps() из data.table.

tolerance = 5e-6
#create ranges to join on
DT1[, `:=`(min = mz - mz * tolerance, 
           max = mz + mz * tolerance) ]
DT2[, `:=`(min = iso_mz - iso_mz * tolerance, 
           max = iso_mz + iso_mz * tolerance) ]
#set keys
setkey(DT1, min, max )
setkey(DT2, min, max )
#perform overlap join, order, remove min-max columns
ans <- setorder( foverlaps( DT2, DT1 ), mz)[, `:=`(min=NULL,max=NULL,i.min=NULL,i.max=NULL)][]


# mz     Var1     Var2   iso_mz comp
# 1: 433.2315 433.2310 433.2321 433.2335   m4
# 2: 451.0920 451.0914 451.0926 451.0900   m1
# 3: 490.1672 490.1664 490.1683 490.1651   m2
# 4: 518.2259 518.2251 518.2266 518.2281   m3


#check
all.equal( setcolorder(ans, names(Output)), Output )
[1] TRUE
...