aide pour création de functions simples mais utiles

Version imprimable

22/05/2007, 14h58
couanp

aide pour création de functions simples mais utiles

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.

Voici la solution vectorisée :

Code:

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 ;)

22/05/2007, 15h12
couanp

Bonjour

message d'erreur suivant:

Citation:

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

à quoi cela est il du?

Merci
22/05/2007, 15h23
Jerome Briot

Ajoute

Code:

idx=double(isnan(M));
22/05/2007, 15h36
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.
22/05/2007, 15h38
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:

idx=double(~isnan(M));
22/05/2007, 15h44
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
22/05/2007, 16h00
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.
22/05/2007, 16h32
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.
22/05/2007, 16h57
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.

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

Code:

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 !

23/05/2007, 10h53
couanp

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

Merci encore