Использование mapply () для сравнения строк столбцов

Рассмотрим вложенный семейный вызов apply:

• mapply() - внешний: попарная итерация между соответствующими столбцами Data_1 и Data_2
• sapply - inner: векторная итерация каждого значения в столбце Data_1 для сравнения элементов

Ниже проверяется, находится ли каждое значение Data_1 между двумя значениями каждого столбца Data_2:

Данные

txt = '  data1       data2       data3      data4
-0.710003   -0.714271   -0.709946   -0.713645
-0.710458   -0.715011   -0.710117   -0.714157
-0.71071    -0.714048   -0.710235   -0.713515
-0.710255   -0.713991   -0.709722   -0.713972'

Data_1 <- read.table(text=txt, header=TRUE)

txt = ' dataA       dataB       dataC      dataD
-0.71097    -0.714059   -0.70928    -0.714059
-0.710343   -0.714576   -0.709338   -0.713644'

Data_2 <- read.table(text=txt, header=TRUE)

Код

check_inbetween <- function(x,y){
  sapply(x, function(i) (i > y[1] & i < y[2]))
}

inbetween_matrix <- mapply(check_inbetween, Data_1, Data_2)

inbetween_matrix
#      data1 data2 data3 data4
# [1,] FALSE FALSE FALSE  TRUE
# [2,]  TRUE FALSE FALSE FALSE
# [3,]  TRUE FALSE FALSE FALSE
# [4,] FALSE FALSE FALSE  TRUE

Это одно решение, основанное на data.table, но если вы используете data.frame или матрицу, его будет легко адаптировать. Чтобы достичь того, чего вы хотите, вы должны использовать один lapply над другим lapply. Верхний из них выполняет итерацию по столбцам, вызывая другой цикл по строкам.

library(data.table)

# it gets all elements of column 'j' to do diff computation
get_column_diff <- function(dt_1, dt_2, j){

        get_point_diff <- function(i){
                # it should return a vector with all differences 
                # in comparison of the point (i,j)
                unlist(dt_1[i, ..j]) - unlist(dt_2[, ..j])
        }


        i_rows <- 1:nrow(dt_1)
        lapply(X=i_rows, FUN=get_point_diff)

}

j_cols <- 1:ncol(Data_1)
lapply(FUN=get_column_diff, dt_1=Data_1, dt_2=Data_2, X=j_cols)

Функция возвращает список из списка, каждый элемент которого является результатом столбца, а его элемент - списком с результатом строки.

Что касается прироста скорости, я не могу сказать, насколько он будет быстрым без сравнения тестов, но, вероятно, он будет быстрее на любом цикле.