Discussion :

[powdaddy]

Bonjour, je souhaiterais effectuer un tri conditionnel sur un array de dimension 3

dans l'exemple ci dessous, j'effectue la somme des lignes de chacune des deux matrices (4x3) indicées par la 3ième dimension, puis j'ordonne la somme de ces lignes. Je souhaiterais ensuite permuter les colonnes de chaque matrice sur la base de l'ordre obtenu (et propre à chaque matrice)

Code :

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

1
2
3
4
5
 
t<-rnorm(24)
dim(t)<-c(4,3,2)
t1<-apply(t,c(2,3),sum)
cond<-t(apply(t1,2,order))

en dimension 2 je n'ai pas de difficultés je procèderais ainsi:

Code :

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

1
2
3
4
 
t2<-t[,,1]
cond2<-apply(t2,2,sum)
t2[,order(cond2)]

mon problème: je ne trouve pas d'expression equivalente à utiliser pour mon array. Le code suivant ne marche pas:

Code :

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

1
2
 
t[,order(cond)]

Merci de votre aide!

[powdaddy]

Je viens de trouver cette solution qui reduit la dimension de mon tableau
Ca marche bien mais si qqn a une solution en restant en dimension 3, je suis intéressé!

Code :

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

1
2
3
4
5
6
 
dim(t)<-c(4,6)
good<-rbind(t,apply(t,2,sum),rep(1:2,each=3))
good[,order(good[6,],good[5,])]
final<-good[-c(5,6),]
dim(final)<-c(4,3,2)

[vchouraki]

Bonsoir,

Une boucle for?

HTH

Vincent

[powdaddy]

Merci pour votre aide,
Je suis arrivé à la solution d'utiliser une boucle for.
En terme d'efficacité je ne suis pas certain que ce soit pire qu'un apply
En revanche c'est moins élégant...

[vchouraki]

Bonjour,

Oui c'est pareil en terme de performance sur les versions récentes de R.

Vincent