Простая комбинаторика в R

Я хочу поделиться функцией R для нахождения всех возможных уникальных ненаправленных комбинаций между элементами одного вектора:

combi <- function(vec1)
{
  si <- length(vec1)
  first <- rep(vec1, (si-1):0)
  secR <- rev(vec1)
  second <- secR[sequence(1:(si-1))]
  second <- rev(second)
  combi <- matrix(cbind(first, second), ncol = 2)
  return(combi)
}

и спросить, есть ли более простой способ сделать это? (Мне нужно, чтобы результат был в 2-столбцовой матрице).