Умножение матриц с использованием NumericMatrix и NumericVector в Rcpp

Мне интересно, есть ли способ расчета умножения матриц с использованием классов NumericMatrix и NumericVector. Я задаюсь вопросом, есть ли какой-нибудь простой способ помочь мне избежать следующего цикла для проведения этого вычисления. Я просто хочу рассчитать X% *% бета.

// assume X and beta are initialized and X is of dimension (nsites, p), 
// beta is a NumericVector with p elements. 
for(int j = 0; j < nsites; j++)
 {
    temp = 0;

    for(int l = 0; l < p; l++) temp = temp + X(j,l) * beta[l];

}

Заранее большое спасибо!

Ответы на вопрос(1)

Ваш ответ на вопрос