1

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))
  );
}

ありがとうございました

4

1 に答える 1