library(gregmisc)

minor2 = function(mat, idx){
  little_mat = mat[idx,idx]
  return(det(little_mat))
}

index_vec = combinations(nrow(mat),2)
minors  = apply(index_vec, MARGIN=1, FUN=function(ix){ minor2(mat,ix)})