j'ai une norme qui est décrite par la matrice sigma
multiplication norme/matrice vectoriced
sigma <- matrix(c(1,0.5,0,0.5,1,0,0,0,1),3,3))
Pour calculer la norme d'un vecteur I COMPUTE
t(x) %*% sigma %*% x
qui fonctionne très bien pour un vecteur, par exemple x = 1:3
.
Cependant, je veux calculer la norme de nombreux vecteurs en même temps, ce qui est je
x <- t(matrix(rep(1:3, 10),3,10))
(bien sûr rempli d'entrées différentes).
Existe-t-il un moyen de calculer la norme de chaque vecteur simultanément? I.e. quelque chose comme
lapply(1:10, function(i) t(x[i,]) %*% sigma %*% x[i,])
Cela me donne "Erreur dans t (x)% *% sigma: arguments non concordantes" –
Ok, alors x% *% sigma% *% t (X). Cela a fonctionné pour moi. – EugenR
@EugenR Vous calculez des éléments inutiles. – jogo