Discussion :

[boubouchon]

Bonjour,

J'ai une matrice de ce type :
A =
1 2 3
4 5 6
7 8 9
9 8 7
6 5 4
3 2 1

Et j'aimerais obtenir cette matrice (produit des termes pour des couples de lignes) :
A =
4 10 18
63 64 63
18 10 4

Est il possible d'effectuer cette opération sans boucle ?

Merci d'avance !

[ced]

Bonjour,

Il y a certainement plusieurs façons de faire ça sans passer par des boucles.
En voici une :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

A[seq(1:(nrow(A)/2))*2 - 1, , drop=F] * A[seq(1:(nrow(A)/2))*2, , drop=F]

[boubouchon]

Merci beaucoup pour ta réponse !
J'ai bien compris le principe pour deux lignes.
Est-il possible de faire la même opération avec non pas deux lignes mais n lignes ?

[Sengar]

Hello,

Avec n lignes en reutilisant le code de ced:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
n<-3 #met ce que tu veux
s<-lapply(0:(n-1),function(x){A[seq(1:(nrow(A)/n))*n - x, , drop=F]})
res<-s[[1]]
sapply(2:length(s),function(x){res<<-res*s[[x]]})
res