X が単数の場合、次のコードは警告をスローします。ミュートする方法はありますか?
「警告: solve(): システムが特異なようです; 近似解を試みています」
関数:
// [[Rcpp::depends(RcppArmadillo)]]
#include <RcppArmadillo.h>
using namespace Rcpp;
// [[Rcpp::export]]
NumericVector fastLM(const NumericVector & y_, const NumericMatrix& X_) {
const arma::vec & y = as<arma::vec>(wrap(y_));
const arma::mat & X = as<arma::mat>(wrap(X_));
int n = X.n_rows, k = X.n_cols;
arma::colvec coef = arma::solve(X, y);
return(
as<NumericVector>(wrap(coef))
);
}
ありがとうございました