Discussion :

[Ptinéwik]

salut,

Alors voilà.
Je possède des matrices avec de très nombreuses lignes, et 8 colonnes.
Dans cette matrice, j'ai des valeurs manquantes, nan.
Pour mes calculs (régression et corrélation), je souhaite oter ces nan.
Je veux donc créer un programme qui rechercherait les lignes où j'ai un nan, et qu'il me supprime la ligne, et ainsi de suite.
Le problème, c'est que parfois il n'y a qu'un nan sur la ligne (genre il est à la 8ème colonne) et du coup, il faut que je supprime toute la ligne...

Voilà ce que j'ai commencé à écrire, mais ça ne marche pas....

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
% retire les lignes avec coord nan
clear all
close all
load X2
 
[n,p]=size(X2)
 
for i=1:p
    clear f
    f=find(isnan(X2(:,p))==0);
    if f==0;
        X2(:,:)=[];
    elseif X2(:,:);
    end
end

Merci, j'ai passé un grand bout de ma soirée dessus, et j'ai pas trouvé....

[Caro-Line]

Déjà il y a une erreur dans ton code :

Code :

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

X2(:,:)=[]

vide toute ta matrice !!
Et en plus tu parcours les colonnes.
Pourtant dans le texte tu dis que tu veux supprimer les lignes

Une solution qui supprime les lignes dans X2 dès qu'il y a au moins 1 NaN :

Code :

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

1
2
3
4
5
6
7
%on parcours toutes les lignes de X2 depuis le "bas"
for i=n:-1:1
   if any(isnan(X2(i,:)) %il y a au moins un NaN dans la ligne i de X2
      %=> on la supprime
      X2(i,:)=[];
   end
end

PS : je n'ai pas testé, pas le temps de lancer MATLAB...

[Ptinéwik]

Ok merci, je vais tester.

en fait je pensais que

Code :

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

X2(:,:)=[]

me vidait la ligne avec toutes les colonnes de l'endroit où j'avais mon nan...

Donc je vais voir si ca marche et je te tiens au courant.
Merci de ton aide

[Ptinéwik]

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
21
22
23
24
25
26
27
 
%30 Janvier 2008
%
%On retire les lignes avec les valeurs nan
 
X=fichier;
 
function Y=isnan(X)
 
[n,p]=size(X);
%on parcours toutes les lignes de X depuis le "bas"
 
tt = 'Progress [  0%% ]';
h = waitbar(0,tt);
fprintf(tt);
 
fin=???
 
for i=n:-1:1
   fprintf('\b\b\b\b\b\b%2.0f%% ]\n', 100*i/fin);
   if any(isnan(X(i,:))) %il y a au moins un NaN dans la ligne i de X
      %=> on la supprime
      X(i,:)=[];
   end
    waitbar(i/fin,h,sprintf('Progress [  %2.0f%% ] \n', 100*i/fin));
    pause(0.1);
end

Alors voilà ma fonction qui fonctionne bien merci.

Par contre j'ai voulu intégrer une barre de progression, mais j'ai un petit soucis, ptètre tout bête, c'est que je ne sais pas comment définir fin, ce qui fait que du coup cela ne fonctionne pas.....

[vinc-mai]

On peux eviter la boucle, je pense:

Code :

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

1
2
3
 
l=any(isnan(X2),2)
X2(l,:)=[]

[Ptinéwik]

Je vais essayer sans la boucle,
Mais est-ce que quelqu'un a une idée pour la progression?

[Caro-Line]

Comme i va de n à 1 il faut changer le pas pour prendre ça en compte :

Code :

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

fin=n;%tout simplement ;-) donc inutile....

et

Code :

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

1
2
3
4
5
for i=n:-1:1
   %blabla
   pas = (n-i+1)/n;
   waitbar(pas,h,sprintf('Progress [  %2.0f%% ] \n', pas*100));
end

Et sinon effectivement la solution de vinc-mai est la plus jolie

[Ptinéwik]

Merci, vous êtes géniaux