Метод R для сравнения NA между двумя векторами

Я пытаюсь сравнить два вектора для идентичных элементов, скажем, для векторов a и b ниже.

# vectors a and b
a <- c(1, 2, 3, NA)
b <- c(1, 3, NA, NA)

# compare for identity
a == b

## [1]  TRUE FALSE    NA    NA

Однако я бы хотел, чтобы сравнение включало NA, чтобы, когда оба элемента в данной позиции имеют значение NA, результат будет TRUE, но если только один из них равен NA, то это будет FALSE. Для этого я написал функцию ниже.

# function to include NA's in the comparison
compare_incl_NA <- function(vec1, vec2){
  res = ifelse(is.na(vec1) | is.na(vec2), 
               is.na(vec1) & is.na(vec2), 
               vec1 == vec2)
  res
}

# run function
compare_incl_NA(a,b)

## [1]  TRUE FALSE FALSE  TRUE

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