Discussion :

[Newenda]

Bonjour, imaginons :

Code :

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

1
2
3
 
temps = [1 1 2 2 3 3 3 4 4] %(en jour)
data= [741 nan 458 42 48 78 nan 85 95]

J'aimerais faire la somme des éléments de data du même jour, ie me retouver avec une matrice :

Code :

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

1
2
 
sum_data_jour = [741 458+42 48+78 85+95]

donc j'utilise UNIQUE :

Code :

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

1
2
3
4
5
temps_uni = unique(temps)

for i_uni = 1:length(temps_uni)
data(i_uni) = sum(data(temps==temps_uni(i_uni)));
end

le problème étant que la SUM des éléments d'un matrice contenant un NaN est toujours égale à NaN...pas terrible.

Je pourrais sûr rajouter cette commande avant pour transformer mes NaN en 0 :

Code :

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

1
2
 
data(isnan(data))=0;

Mais n'y a t il pas un moyen de faire un double indexage logique au lieu de tout ça? un truc du genre :

Code :

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

data(i_uni) = sum(data(data==temps_uni(i_uni) & isnan(data)=0));

Merci

[duf42]

C'est presque çà:

Code :

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

1
2
 
data(i_uni) = sum(data(data==temps_uni(i_uni) & ~isnan(data)));

Ou en une seule ligne:

Code :

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

1
2
 
arrayfun(@(x) (sum(data(data==x & ~isnan(data)))), unique(temps));

Enfin faut tester...

[Newenda]

je vais essayer merci !

[Jerome Briot]

Code :

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

1
2
idx = ~isnan(data);
data = full(sparse(temps(idx),1,data(idx)));

Il est souvent inutile d'utiliser FULL pour revenir à une matrice pleine...

[Newenda]

Code :

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

1
2
idx = ~isnan(data);
data = full(sparse(temps(idx),1,data(idx)));

Il est souvent inutile d'utiliser FULL pour revenir à une matrice pleine...

utile?

Et effectivement, cette petite ligne me simplifie bien la vie !

merci

[Jerome Briot]

utile?

j'ai bien dit "inutile"...