Rcpp：从矩阵中删除列和行[重复]

Question

这个问题已经在这里有了答案：

• RCPP NumericMatrix-如何删除行/列？ 2个答案

我正在尝试创建一个函数，该函数采用矩阵nxp和索引e ，并返回通过从X消除第e列和第e行获得的子矩阵。我认为最简单的方法是创建一个n -1 xp-1矩阵，并在其中插入第e行和第e列形成的交叉周围的角。 使用Corner语法，同样的方法对于EigenRcpp也可以正常工作。 似乎Rcpp不喜欢分配给Ranges。 我收到以下错误消息：

error: non-static reference member 'Rcpp::SubMatrix<14>::MATRIX& Rcpp::SubMatrix<14>::m', can't use default assignment operator

我复制下面的代码。

我的问题：如何为矩阵块分配值？ 有一个更好的方法吗？

#include <Rcpp.h>
using namespace Rcpp;

NumericMatrix elimMat(NumericMatrix X, int e){

int p = X.cols() - 1;
int n = X.rows() - 1;
int f = e - 1;
NumericMatrix M = NumericMatrix(n, p);

M(Range(0, f - 1), Range(0, f - 1)) = X(Range(0, f - 1), Range(0, f - 1)); //TopLeft same
M(Range(0, f - 1), Range(f, p - 1)) = X(Range(0, f - 1), Range(f + 1, p)); //TopRight
M(Range(f, n - 1), Range(0, f - 1)) = X(Range(f + 1, n), Range(0, f - 1)); //BottomLeft
M(Range(f, n - 1), Range(f, p - 1)) = X(Range(f + 1, n), Range(f + 1, p)); //BottomRight

return M;
}

谢谢您的帮助，Marco

Answer 1

关于这个错误，

错误：非静态引用成员'Rcpp :: SubMatrix <14> :: MATRIX＆Rcpp :: SubMatrix <14> :: m'，无法使用默认赋值运算符

我不认为您（或客户端上的任何人）可以做任何事情，因为问题来自SubMatrix的类定义 。 根据编译器，

private:
    MATRIX& m ;
    vec_iterator iter ;
    int m_nr, nc, nr ;

m应该是static成员，不是，因此是错误。

---

尽管我无法为您提供以这种方式使用Range分配的一般解决方法，但是对于解决有关“ NumericMatrix的一部分的特定问题，这是一种可行的方法：

#include <Rcpp.h>

// [[Rcpp::export]]
Rcpp::NumericMatrix crossout(Rcpp::NumericMatrix X, int e) {
  Rcpp::NumericMatrix result = Rcpp::no_init_matrix(X.nrow() - 1, X.ncol() - 1);
  Rcpp::NumericMatrix::iterator src = X.begin(), dst = result.begin();

  while (src != X.end()) {
    if (((src - X.begin()) % X.nrow()) != e && ((src - X.begin()) / X.nrow()) != e) {
      *dst = *src;
      ++dst;
    }
    ++src;
  }

  return result;
}

---

/*** R

M <- matrix(1:9 + .5, nrow = 3)
R> all.equal(M[c(1, 3), c(1, 3)], crossout(M, 1))
#[1] TRUE

R> crossout(M, 1)
#     [,1] [,2]
#[1,]  1.5  7.5
#[2,]  3.5  9.5

M2 <- matrix(1:20 + .5, nrow = 4)
R> all.equal(M2[c(1:2, 4), c(1:2, 4:5)], crossout(M2, 2))
#[1] TRUE

R> crossout(M2, 2)
#     [,1] [,2] [,3] [,4]
#[1,]  1.5  5.5 13.5 17.5
#[2,]  2.5  6.5 14.5 18.5
#[3,]  4.5  8.5 16.5 20.5

R> any(unique(c(M2[3,], M2[,3])) %in% crossout(M2, 2))
#[1] FALSE

*/

---

上面，如果源迭代器的位置不在e指定的列或行之内，则将“源”迭代器的当前值复制到“目标”迭代器。 以(src - X.begin()) % X.nrow()获得行索引，以(src - X.begin()) / X.nrow()获得列索引。

---

编辑：关于下面的评论，为行和列尺寸适应非平方矩阵/不同索引将非常简单：

#include <Rcpp.h>

// [[Rcpp::export]]
Rcpp::NumericMatrix crossout2(Rcpp::NumericMatrix X, int r, int c) {
  if (r >= X.nrow() || c >= X.ncol()) {
    Rf_error("Invalid row or column index.\n");
  }

  Rcpp::NumericMatrix result = Rcpp::no_init_matrix(X.nrow() - 1, X.ncol() - 1);
  Rcpp::NumericMatrix::iterator src = X.begin(), dst = result.begin();

  while (src != X.end()) {
    if (((src - X.begin()) % X.nrow()) != r && ((src - X.begin()) / X.nrow()) != c) {
      *dst = *src;
      ++dst;
    }
    ++src;
  }
  return result;
}

---