Сравните четыре числовых вектора с интервалом допуска и сообщите общие значения

Question

У меня есть четыре больших вектора неравной длины . Ниже я предоставляю набор игрушечных данных, аналогичный моему исходному набору данных:

a <- c(1021.923, 3491.31, 102.3, 12019.11, 879.2, 583.1)
b <- c(21,32,523,123.1,123.4545,12345,95.434, 879.25, 1021.9,11,12,662)
c <- c(52,21,1021.9288,12019.12, 879.1)
d <- c(432.432,23466.3,45435,3456,123,6688,1021.95)

Есть ли способ сравнить все эти векторы один за другим с допустимым порогом ± 0,5 для совпадения? Другими словами, я хочу сообщить о числах, которые являются общими для всех четырех векторов, в то время как допускает отклонение 0,5 .

В случае набора данных игрушек, приведенного выше, окончательный ответ:

    Match1
a 1021.923
b 1021.900
c 1021.929
d 1021.950

Я понимаю, что это возможно для двух векторов, но как я могу сделать это для 4 векторов?

СВЯЗАННЫХ

• Все-во-всем setdiff для двух числовых векторов с числовым порогом для принятия совпадений
• Сравните два вектора чисел на основе порога допуска (±) 0,5

Answer 1

Вот решение для data.table.

Он масштабируется до n векторов, поэтому попробуйте подавать его сколько угодно. Он также хорошо работает, когда несколько значений имеют «совпадения» во всех векторах.

образец данных

a <- c(1021.923, 3491.31, 102.3, 12019.11, 879.2, 583.1)
b <- c(21,32,523,123.1,123.4545,12345,95.434, 879.25, 1021.9,11,12,662)
c <- c(52,21,1021.9288,12019.12, 879.1)
d <- c(432.432,23466.3,45435,3456,123,6688,1021.95)

код

library(data.table)

#create list with vectors
l <- list( a,b,c,d )
names(l) <- letters[1:4]
#create data.table to work with
DT <- rbindlist( lapply(l, function(x) {data.table( value = x)} ), idcol = "group")
#add margins to each value
DT[, `:=`( id = 1:.N, start = value - 0.5, end = value + 0.5 ) ]
#set keys for joining
setkey(DT, start, end)
#perform overlap-join
result <- foverlaps(DT,DT)

#cast, to check how the 'hits' each id has in each group (a,b,c,d)
answer <- dcast( result, 
             group + value ~ i.group, 
             fun.aggregate = function(x){ x * 1 }, 
             value.var = "i.value", 
             fill = NA )

#get your final answer
#set columns to look at (i.e. the names from the earlier created list)
cols = names(l)
#keep the rows without NA (use rowSums, because TRUE = 1, FALSE = 0 )
#so if rowSums == 0, then columns in the vactor 'cols' do not contain a 'NA'
answer[ rowSums( is.na( answer[ , ..cols ] ) ) == 0, ]

выход

#    group    value        a      b        c       d
# 1:     a 1021.923 1021.923 1021.9 1021.929 1021.95
# 2:     b 1021.900 1021.923 1021.9 1021.929 1021.95
# 3:     c 1021.929 1021.923 1021.9 1021.929 1021.95
# 4:     d 1021.950 1021.923 1021.9 1021.929 1021.95