Multiplicação de matrizes usando NumericMatrix e NumericVector no Rcpp

Gostaria de saber se existe uma maneira de calcular a multiplicação de matrizes usando as classes NumericMatrix e NumericVector. Gostaria de saber se existe alguma maneira simples de me ajudar a evitar o seguinte loop para realizar esse cálculo. Eu só quero calcular X% *% beta.

// 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];

}

Muito obrigado antecipadamente!

questionAnswers(1)

yourAnswerToTheQuestion