Discussion :

[ju_cayenne]

Bonjour,

J'ai une matrice 100 (lignes) x 20 (colonnes): DATA

Je moyenne les colonnes entre elles pour avoir 2 vecteurs:
- le vecteur 1 correspond à la moyenne de la colonne 1 à 10
- le vecteur 2 correspond à la moyenne de la colonne 11 à 20

Code :

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

1
2
3
 
MEAN_1=nanmean(DATA(:,1:10),2)
MEAN_2=nanmean(DATA(:,11:20),2)

J'aimerai calculer la moyenne en supprimant à chaque fois 2 colonnes
ca donnerai :

- le vecteur 1a correspond à la moyenne de la colonne 2 à 10
- le vecteur 2a correspond à la moyenne de la colonne 12 à 20

- le vecteur 1b correspond à la moyenne de la colonne 1 à 10 sans la colonne 2
- le vecteur 2b correspond à la moyenne de la colonne 11 à 20 sans la colonne 12

- le vecteur 1c correspond à la moyenne de la colonne 1 à 10 sans la colonne 3
- le vecteur 2c correspond à la moyenne de la colonne 11 à 20 sans la colonne 13

etc

Je pense à une boucle dans une autre, mais je n'arrive pas bien à voir comment faire. J'aimerai automatiser la suppression de colonne (la 1 avec la 11, la 2 avec la 12, ... la 10 avec la 20) car ma vraie matrice fait 200 colonnes
Merci

[le fab]

salut

tout simplement comme ca :

Code :

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

1
2
3
4
5
MEAN_1b=nanmean(DATA(:,[1 3:10]),2)
MEAN_2b=nanmean(DATA(:,[11 13:20]),2)
 
MEAN_1c=nanmean(DATA(:,[1 2 4:10]),2)
MEAN_2c=nanmean(DATA(:,[11 12 14:20]),2)

cf les chapitres indexage et concaténation de ce tuto

Fabien

[ju_cayenne]

Oui mais j'ai en réalité 200 colonnes, donc je ne vais pas faire 100 lignes de commande. Je voulais un moyen pour automatiser

[le fab]

comme ca alors (si je comprend bien)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
% vecteur d'index
vIdx = 1:100;
% pré allocation mémoire vecteurs résultats
MEAN_1 = nan(size(vIdx));
MEAN_2 = nan(size(vIdx));
% boucle de calcul
for ii = 1:vIdx
   % calcul moyenne des 100 1eres colonnes sans la iieme
   MEAN_1(ii) = nanmean(DATA(:,vIdx(vIdx ~=ii),2)
   % calcul moyenne des 100  colonnes suivantes sans la 100+iieme
   MEAN_2(ii) = nanmean(DATA(:,vIdx(vIdx ~=ii) + 100,2)
end