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)})