在构造时用单个值填充 NumericMatrix

Question

我试图在构造时用单个值填充 NumericMatrix。例如，请考虑以下内容：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
void test() {
  NumericMatrix res(1, 1, NA_REAL);
}

这是抛出以下错误：

error: call to constructor of 'Vector<14, PreserveStorage>' is ambiguous
        VECTOR( start, start + (static_cast<R_xlen_t>(nrows_)*ncols) ),
        ^       ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
file46e92f4e027d.cpp:6:17: note: in instantiation of function template specialization 'Rcpp::Matrix<14, PreserveStorage>::Matrix<double>' requested here
  NumericMatrix res(1, 1, NA_REAL);

/Library/Frameworks/R.framework/Versions/4.0/Resources/library/Rcpp/include/Rcpp/vector/Vector.h:88:5: note: candidate constructor [with T = double]
    Vector( const T& size, const stored_type& u,
    ^
/Library/Frameworks/R.framework/Versions/4.0/Resources/library/Rcpp/include/Rcpp/vector/Vector.h:211:5: note: candidate constructor [with InputIterator = double]
    Vector( InputIterator first, InputIterator last){
    ^

为什么 NumericMatrix 不能在固定尺寸的情况下用单个值实例化？

Answer 1

简而言之，这是有效的（一条较长的线被分成三段以供显示）：

> Rcpp::cppFunction("NumericVector fp() { 
+     NumericVector res(3,NA_REAL); 
+     return res;}")
> fp()
[1] NA NA NA  
>

但是没有匹配的构造函数使用 rows, cols 矩阵。所以你必须使用上面给你的向量，并手动设置尺寸。

例如通过（我把所有内容都放在一行中，我在这里分手进行说明）

> Rcpp::cppFunction("NumericMatrix fp(int n, int k) { 
+         NumericVector res(n*k,NA_REAL); 
+         res.attr(\"dim\") = IntegerVector::create(n,k); 
+         return NumericMatrix(res);}")
> fp(2,3)
     [,1] [,2] [,3]
[1,]   NA   NA   NA
[2,]   NA   NA   NA
>

Answer 2

不是为了篡夺德克，但没必要用.attr()设置矩阵的维度。

与向量不同，填充矩阵需要为构造函数提供具有 n * p 个元素和维度的迭代器。

Matrix(const int& nrows_, const int& ncols, Iterator start)

其他构造函数见：inst/include/Rcpp/vector/Matrix.h

考虑到这一点，原来的例子可以改成：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
Rcpp::NumericMatrix matrix_fill_by_vec(int n, int p) {

  // fill matrix using a vector
  Rcpp::NumericVector A = Rcpp::NumericVector(n * p, NA_REAL); 
  Rcpp::NumericMatrix B = Rcpp::NumericMatrix(n, p, A.begin());

  return B;
}

试驾一下，我们得到：

matrix_fill_by_vec(3, 2)
#      [,1] [,2]
# [1,]   NA   NA
# [2,]   NA   NA
# [3,]   NA   NA

在构造时用单个值填充 NumericMatrix

Fill a NumericMatrix with a single value on construction

rcpp