rcpp - 在 Rcpp 中将矩阵初始化为 NA

Question

有一种方法可以使用 NA 值初始化数值向量，例如。

NumericVector x(10,NumericVector::get_na())

有没有类似的方法可以将矩阵初始化为 NA 值？

Answer 1

这是一个不浪费内存的版本。

#include <Rcpp.h>
using namespace Rcpp ;

// [[Rcpp::export]]
NumericMatrix na_matrix(int n){
  NumericMatrix m(n,n) ;
  std::fill( m.begin(), m.end(), NumericVector::get_na() ) ;
  return m ;
}

FWIW，在 Rcpp11 ，您可以使用一些更具表现力的语法:

NumericMatrix m(n,n, NA) ;

感谢这个 constructor