Discussion :

[ju_cayenne]

Bonjour à tous,

J'ai une matrice 2 colonnes (une DATE sous la forme AAAAMMJJ et la seconde des valeurs)
J'aimerai que lorsque le nombre de jour dans le mois est inférieur à 25 on supprime toutes les données du mois

EXemple on a :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
20090301	6,15
20090302	14
20090303	14,4
20090304	7,3
20090305	14
20090801	1,65
20090802	0
20090803	0
20090804	8,7
20090805	5,5
20090806	0
20090807	0
20090808	3,55
20090809	2,3
20090810	0
20090811	0
20090812	11,25
20090813	0
20090814	0,2
20090815	0

du coup le mois de mars ne compte que 5 jours j'aimerai supprimer de ma matrice ces cinq lignes. Il faudrait trouver un moyen d'automatiser car j'ai plusieurs séries, certaines bonnes, d'autres où les mois à supprimer sont différents ...

Merci

[Invité]

Bonjour,

On va commencer par obtenir le mois seul : (data étant le tableau que tu présentes dans ton message)

Code :

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

mois = mod(round(data(:,1)/100),100);

Puis on calcule le nombre de jours dans chaque mois avec les fonctions unique + accumarray (voir FAQ Comment regrouper des éléments selon des valeurs communes ? pour plus de détails)

Code :

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

1
2
[unik, ~, ic] = unique(mois);
nbJours = accumarray(ic(:), mois, [numel(unik) 1], @numel);

On repère les mauvais mois :

Code :

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

1
2
minJours = 25;
mauvais_mois = unik(nbJours < minJours);

Et on les supprime des données initiales:

Code :

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

data(ismember(mois, mauvais_mois), :) = [];

[ju_cayenne]

Bonjour,

J'ai bien compris la logique du code, mais j'ai 10 ans donc 10 fois le mois de janvier, 10 fois le mois de février ....
"Unik" est un vecteur de 1 à 12, et donc "nbJours" somme tous les jours de janvier des 10 ans soit 310.
Donc lorsque j'applique

Code :

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

mauvais_mois = unik(nbJours < minJours);

j'ai rien inférieur à 25.

[Invité]

Effectivement, dans ce cas, on va prendre le mois, mais aussi l'année :

Code :

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

mois = mod(round(data(:,1)/100),100);

[ju_cayenne]

C'est super merci beaucoup