Discussion :

[couanp]

Bonjour

j'ai la matrice xPrice (fichier matlab) et je cherche à ecrire une fonction qui me retourne dans une matrice de dimension (1 ligne* nombre de colonnes) , le nombre de NaN ( données manquantes) à la fin de chaque colonne. Si par exemple j'ai 3 valeurs manquantes successives à la fin de la colonne 2, je devrais avoir dans ma matrice retournée 3 au niveau de cette colonne. Si il n ya pas de NaN à la fin de la colonne, la valeur 0 est retournée.

merci de votre aide.

[Jerome Briot]

Voici la solution vectorisée :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
M=[  nan  2  nan 
     4  5  nan
     7  nan  9
    10 11 12
    13 nan 15
    nan nan nan
    19 nan nan]
 
idx=isnan(M);
 
NaNs_au_debut=sum(cumprod(idx))
 
NaNs_a_la_fin=sum(cumprod(idx(end:-1:1,:)))

A étudier en détail... comme d'habitude

[couanp]

Bonjour

message d'erreur suivant:

??? Undefined function or method 'cumprod' for input arguments of type 'logical'.

à quoi cela est il du?

Merci

[Jerome Briot]

Ajoute

Code :

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

idx=double(isnan(M));

[couanp]

Ca marche, par contre j'aimerais bien savoir en quoi l'utilisation du double a permis de régler le probleme.

Aurais tu une solution concernant la fonction qui renvoie la date de début pour les premieres valaurs non NaN dans chaque colonne et les dates de fin dans chaque colonne pour les dernieres valeurs non NaN?

Merci de ton aide.

[Jerome Briot]

La majeur partie des fonctions Matlab ne fonctionne bien qu'avec les variables de type Double, c'est tout.

Pour étudier les non-NaN, tu fais simplement cette modification :

Code :

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

idx=double(~isnan(M));

[couanp]

Pour les non-NaN, je ne cherche pas les sommes, mais plus tot, les dates correspondantes, donc j'utilise aussi la matrice xDates.

Ce que je cherche c'est comment récupérer la date correspondante à cet index, alors que MATLAB stocke les dates dans une matrice (xDates) et mes prix dans une autre(xPrice).

Merci

[coolzy]

si tes deux matrices sont organisees de la meme maniere, c'est a dire que tu peux passer d'un element Nan a la date correspondante en connaissant l'indice de ligne et de colonne, tu peux utiliser find pour trouver les coordonnees/indices de tes elements Nan dans l'une, puis aller lire la date correspondante dans l'autre matrice en utilisant ces coordonnees.

[coolzy]

Si tu as d'un cote une matrice A avec dedans tes dates. Si tu as de l'autre cote une matrice B avec dedans les valeurs, dont les NaN.
Si tu sais que B(i,j) = NaN, tu n'as plus qu'a recuperer A(i,j) pour avoir la date correspondante.

pour avoir les indices i,j de l'element Nan :
regarde Find dans Help de matlab.
[r,c] = find (X) renvoie dans r et c les indices des elements non nuls de X. Mais tu peux chercher d'autres elements par exemple find (X ==0) pour trouver les elements nuls. Je ne sais pas si find (X == NaN) marche, il faut tester.

[couanp]

je cherche à afficher les dates de debut pour chaque colonne de la matrice xPrice lorsque je suis à la premiere cellule qui n'est pas NaN.

je cherche egalement à afficher la date fin de ces meme colonnes lorsque je tombe sur la derniere valeur de chaque colonne qui n'est pas NaN.

Je n'y arrive toujours pas car, je n'arrive pas à travailler sur les deux matrices malgré les explications.

Je joins les deux fichiers avec mes deux matrices: xPrice qui contient mes données avec des NaN et xDates et textdata qui contiennent les dates.

Merci de votre aide.

[Caro-Line]

ça devrait marcher...si j'ai bien compris ton problème...

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
 
%initialisations
dates_debut=zeros(1,size(xPrice,2));
dates_fin=dates_debut;
%boucle sur les colonnes de xPrice
for iCol=1:size(xPrice,2)
    %la colonne étudiée
    xCol =xPrice(:,iCol);
    %Les emplacements des non-nan
    idx = find(~isnan(xCol));
    %le premier indice donne la première date
    if length(idx)>0 %on verifie quand meme qu'il y en a au moins un
        dates_debut(iCol)=xDates(idx(1));
        dates_fin(iCol)=xDates(idx(end));
    end    
end

Mais l'information qui manquait depuis le début de ton post c'est que xDates est un vecteur pas une matrice !

[couanp]

Cela marche. Par contre j'ai remplacé xDates par textdata pour avoir directement les vraies dates: par exemple "21/12/2004"

Merci encore