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