Как использовать Rcpp для реализации аналогичной функции для регрессионного анализа функций lm в R? ⇐ C++

[Anonymous]

Чтобы ускорить расчет очень больших матриц, я использую Rcpp для выполнения вычислений, которые включают в себя регрессионный анализ. В документации r я заметил, что lm использует метод декомпозиции qr. Я пробовал программировать, и код регрессионного анализа выглядит следующим образом.
// [[Rcpp::depends(RcppArmadillo)]]
// [[Rcpp::depends(bigmemory)]]

#include
#include
#include

using namespace std;
using namespace Rcpp;
using namespace arma;

List lm_armadillo_with_p_values(arma::mat mat) {
// seprate matrix to y and x
arma::mat X_mat = mat.cols(1, mat.n_cols - 1);
arma::colvec y_vec = mat.col(0);
// Rcpp::Rcout

Подробнее здесь: https://stackoverflow.com/questions/792 ... ysis-of-lm