在RCPP/RcppArmadillo中，我想要做以下事情:从一个n x n矩阵A中，我想提取一个子矩阵A[-j，-j]，其中j是一个索引向量:在R中，它可以是这样的
A = matrix(1:16, 4, 4)
j = c(2, 3)
A[-j, -j]
似乎这个功能在RCPP或RcppArmadillo中不可用-抱歉，如果我忽略了什么.R中的一种方法是
pos = setdiff(1:nrow(A), j)
A[pos, pos]
这将延续到RcppArmadillo，但似乎必须创建向量pos作为j的补充-我不确定如何有效地做到这一点.
有没有人有一个有效的实现或一段代码可以分享的 idea ？
Aradillo docs具有函数.shed，该函数接受一个或多个自变量，该自变量为"……包含要移除的行/列/片的索引".从我的读数中删除行和列都将接受.shed()的两个调用.
使用您的示例
// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
// [[Rcpp::export]]
arma::mat fun(arma::mat X, arma::uvec row, arma::uvec col) {
X.shed_cols(col); // remove columns
X.shed_rows(row); // remove rows
return(X);
}
/***R
A = matrix(1:16, 4, 4)
j = c(2, 3)
A[-j, -j]
# minus one for zero indexing
fun(A, j-1, j-1)
*/