Discussion :

[Isabella83]

Bonjour à tous,

J'ai un problème avec ma boucle for + boucle of qui parcours une matrice.
J'ai une matrice du genre :

Code :

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

1
2
3
4
5
6
7
> mat
              m          n
[1,] -3.6989596  1.1571740
[2,] -0.4700084  1.1571740
[3,]  0.2969587  0.2198091
[4,]  0.3898892 -4.7359769
[5,]  5.2133679 -2.6597122

et je voudrais que si pour chaque ligne, la valeur des deux colonnes en valeur absolue est superieur à deux, alors je garde la ligne, sinon je ne la garde pas. J'ai donc creer une matrice, que je rempli si les conditions sont remplies, sinon, non de cette facon

Code :

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

1
2
3
4
5
6
7
8
c <- matrix(data = NA, nrow = 5, ncol = 2)
 
	for(i in 1 :nrow(mat))
		{
		if(abs(mat[i,1])>2 && abs(mat[i,2])>2) 
			c[i,1]<-mat[i,1] 
			c[i,2]<-mat[i,2]
		}

Mais ca ne fonctionne pas, j'obtiens :

Code :

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

1
2
3
4
5
6
         [,1]       [,2]
[1,]       NA  1.1571740
[2,]       NA  1.1571740
[3,]       NA  0.2198091
[4,]       NA -4.7359769
[5,] 5.213368 -2.6597122

quelqu'un pourrait me dire d'où mon erreur vient ?

Merci d'avance

[Hadrien35]

Bonjour.

Je n'ai pas testé sur mon PC mais selon moi, le problème pourrait venir du fait que tu ne défini pas (à l'aide d'accolades) les actions à faire dans le cas ou les conditions de ton if sont respectées.
Du coup, J'ai l'impression que si les conditions sont respectées, la première colonne ne sera pas collée dans ta matrice c, mais la deuxième le sera quand même.
Dans ton exemple, tu n'est sensé conserver que la dernière ligne, et tu as donc bien quatre cellules vides pour la premère colonne, mais bien 4 cellules remplies pour la secondes.
De plus, la matrice c que tu crée sera de même taille que ta matrice de départ, c'est à dire que tu aura une matrice vide sauf pour les lignes ou les conditions sont remplies, est-ce ce que tu cherches?
Voici un petit code qui pourrait, je pense, régler ton problème:

Code :

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

1
2
3
4
5
6
7
8
9
c <- matrix(data = NA, nrow = 5, ncol = 2)
 
for(i in 1 :nrow(mat))
{
	if(abs(mat[i,1])>2 && abs(mat[i,2])>2)
        {
			c[i,]<-mat[i,]
	}
}

Dans le cas ou la matrice C contient des NA pour les lignes ou les conditions ne sont pas remplies. Si tu veux supprimer ces lignes, je te conseille d'ajouter un compteur de ligne du style:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
c <- matrix(data = NA, nrow = 5, ncol = 2)
ligne=1
for(i in 1 :nrow(mat))
{
	if(abs(mat[i,1])>2 && abs(mat[i,2])>2)
        {
			c[ligne,]<-mat[i,]
                        ligne=ligne+1
	}
}
ligne <- NULL

En rajoutant à la fin une commande pour supprimer les lignes de c non remplies par ta fonction (la différence entre ligne et nrow(mat)).
Voila, je ne suis pas garanti du résultat, mais ça devrait marcher je pense.

Cordialement.

[dickoa]

Salut
Je pense que tu l'as lu quelque part ou entendu, mais quand on travaille avec R il faut éviter tant que possible les boucles, et préferer la vectorisation des calculs.
En faisant ça, tu gagne en rapidité (pas toujours mais assez souvent), et le code est plus élégant et simple.
Dans ton cas précis, avoir chaque élément de la colonne supérieure à deux tu peux faire :

Code :

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

1
2
3
4
5
6
7
8
 
 
mat[abs(mat[,1])>2 & abs(mat[,2])>2, ]
 
# mais je te conseille d'utiliser l'argument drop = FALSE pour conserver la structure de matrix
# au cas ou tu n'as qu'une ligne qui respecte ta condition
 
mat[abs(mat[,1])>2 & abs(mat[,2])>2, , drop = FALSE]

Tu comprendras que je n'ai pas envie de voir ce qui cloche dans ta boucle, parce que on est dans R et tu peux t'en passer.