Multiplicación de matrices usando NumericMatrix y NumericVector en Rcpp
Me pregunto si hay una manera de calcular la multiplicación de matrices usando las clases NumericMatrix y NumericVector. Me pregunto si hay alguna manera simple de ayudarme a evitar el siguiente ciclo para realizar este cálculo. Solo quiero 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];
}
¡Muchas gracias por adelantado!