Matrix-Multiplikation mit NumericMatrix und NumericVector in Rcpp

Ich frage mich, ob es eine Möglichkeit gibt, die Matrixmultiplikation mithilfe der NumericMatrix- und der NumericVector-Klasse zu berechnen. Ich frage mich, ob es eine einfache Möglichkeit gibt, die folgende Schleife zur Durchführung dieser Berechnung zu umgehen. Ich möchte nur X% *% Beta berechnen.

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

}

Vielen Dank im Voraus

Antworten auf die Frage(2)

Ihre Antwort auf die Frage