Discussion :

[Isabella83]

Bonjour à tous.
Je vous expose mon problème :
J'ai un tableau de cette forme que j'ai rempli après une étape de simulation (500x):
En gros, chaque i correspond a la valeur tirée, et proba[i] correspond au nombre de fois où i a été tiré lors des 500 fois.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
proba[0]=0
proba[1]=0
proba[2]=49
proba[3]=133
proba[4]=110
proba[5]=90
proba[6]=47
proba[7]=34
proba[8]=20
proba[9]=10
proba[10]=5
proba[11]=2

Je veux ensuite calculer une proba, je fais donc :

Code :

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

1
2
3
4
5
6
for(i=0;i<repeat;i++){
		if(proba[i] != 0){
			double pval = (proba[i])/(float)repeat;
			printf("proba[%d]=%d, p-val = %.3f \n",i,proba[i],pval);
		}
	}

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
proba[0]=0 pval non calculée
proba[1]=0 pval non calculée
proba[2]=49, p-val = 0.098
proba[3]=133, p-val = 0.266
proba[4]=110, p-val = 0.220
proba[5]=90, p-val = 0.180
proba[6]=47, p-val = 0.094
proba[7]=34, p-val = 0.068
proba[8]=20, p-val = 0.040
proba[9]=10, p-val = 0.020
proba[10]=5, p-val = 0.010
proba[11]=2, p-val = 0.004

Je cherche en fait à avoir l'indice de proba pour lequel p-val < 0.05, dans l'exemple, 8.

J'avais donc fait :

Code :

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

1
2
3
4
5
6
7
8
for(i=0;i<repeat;i++){
		if(proba[i] != 0){
			double pval = (proba[i])/(float)repeat;
			if(pval < 0.05){
				break;
				return(i);
			}
		}

Mais cela ne fonctionne pas ....
Comment puis je procéder ?

[Invité]

Je cherche en fait à avoir l'indice de proba pour lequel p-val < 0.05, dans l'exemple, 8.

J'avais donc fait :

Code :

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

1
2
3
4
5
6
7
8
for(i=0;i<repeat;i++){
		if(proba[i] != 0){
			double pval = (proba[i])/(float)repeat;
			if(pval < 0.05){
				break;
				return(i);
			}
		}

Mais cela ne fonctionne pas ....
Comment puis je procéder ?

Bien le bonjour,

Le break a pour effet de briser une séquence. Là en l'occurrence, à cause de ton break, jamais tu ne rentreras dans le return.

[Isabella83]

J'ai modifié en faisant :

Code :

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

1
2
3
4
5
6
7
8
9
10
for(i=0;i<500;i++){
		if(proba[i] != 0){
			double pval = (proba[i])/500;
			if(pval < 0.05){
				return(i);
				break;
			}
		}
 
	}

Mais j'obtiens :

Code :

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

1
2
 
warning: control reaches end of non-void function

[Invité]

Dans ces cas-là, on lit le message d'erreur.... (et on les comprends)

Il est assez explicite pour que je n'ai besoin de te le traduire.

[Isabella83]

Et bien non, je ne comprends pas mon erreur. Enfin je pense que la fonction rencontrera un problème lorsque je n'aurai aucune pval <0.05. et dans ce cas comment puis-je procéder ?
J'ai déclarer la fonction (simulation )contenant cette portion de code comme

Code :

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

int simulation(mes arguments)

Ensuite, dans mon main que j'ai déclaré comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
int main(void)
{	
	srand(time(NULL));
	int mrna_length = 2349 ;
	int exp_tag_nb = 42.08;
	int repeat = 500;
	double mRNA_coverage = 12.2;
	int j;
	j=simulation(mrna_length, mRNA_coverage, exp_tag_nb, repeat);
	printf("valeur de j = %d\n",j);	
	return 0;
 
}

[Isabella83]

Je suis en train de penser, je vais peut etre creer un tableau avec mes proba en valeur.
Ensuite, je vais créer une fonction qui me retourne le premier élément du tableau où valeur < seuil.
Je lui passerai en paramètre un tableau et sa taille, et le seuil .
Je vais essayer, ce serait peut être plus simple.

[supersnail]

Bonjour,

Ton warning vient du fait que ta fonction ne renvoie pas de valeur si elle ne trouve rien (ce qui fait couiner le compilateur évidemment).
Bref en comprenant un minimum l'anglais ça peut servir pour faire de la programmation

Puis ta seconde méthode est effectivement "meilleure" étant donné qu'elle généralise ta recherche à un tableau de n'importe quelle taille et à n'importe quelle valeur de seuil.

[Isabella83]

J'essai de créer une fonction comme je l'ai dit dans mon précédent poste.
Mon problème c'est lorsque aucune de mes valeurs n'est <0.05, et je ne vois vraiment pas comment faire :

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
28
29
30
31
32
33
34
35
36
37
38
#include <stdio.h> 
 
int find(double * ptr, size_t n_elements) 
{ 
	size_t i; 
	int result;
	for(i = 0; i < n_elements; i++){ 
		if(ptr[i]<0.05){
			result=i;
			break;
		}
		else{
			result=0;
		}
	}
	return(result);
 
} 
 
int main() 
{ 
	int f;
	double t[10]; 
	t[0]=0.5;
	t[1]=0.5;
	t[2]=0.5;
	t[3]=0.5;
	t[4]=0.5;
	t[5]=0.5;
	t[6]=0.05;
	t[7]=0.03;
	t[8]=0.02;
	t[9]=0.01;
 
	size_t n_elements = sizeof(t) / sizeof(t[0]); 
	f=find(t,n_elements);
	printf("%d",f); 
}

[supersnail]

Tu peux renvoyer une valeur "impossible" (par exemple un indice négatif), qu'il te suffit de vérifier dans la fonction appelante pour savoir si la fonction a trouvé un élément ou non.

[Isabella83]

J'ai fait comme ceci :

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
#include <stdio.h> 
 
int find(double * ptr, size_t n_elements, double target) 
{ 
	int i=0;
	while((i <= (int) n_elements) && (ptr[i] >= target)){
		i++;	
	}
	if(i < (int)n_elements){
		return i;
	}
	else{
		return -1;
	}
}

Qu'en pensez vous ?

[Invité]

Et toi, qu'en penses-tu ?
As-tu testé/essayé ?

[Isabella83]

je n'avais pas testé tous les cas possibles .....
Finalement j'ai fait :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
int find_index(double * ptr, size_t n_elements, double target) 
{ 
	int i;
	for(i=0;i<(int) n_elements;i++){
		if(ptr[i] < target){
			return(i);
			break;
		}	
 
	}
return(-1);
 
 
}

Et je pense que cette fois ci, c'est la bonne !!!

[Sve@r]

Code c :

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

1
2
3
while((i <= (int) n_elements) && (ptr[i] >= target)){
	i++;	
}

je n'avais pas testé tous les cas possibles

Bonjour
C'est surtout que tu tapais dans ptr[i] quand i = n_elements donc en dehors des valeurs licites du tableau (qui vont de 0 à n-1 !!!)...

Finalement j'ai fait :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
int find_index(double * ptr, size_t n_elements, double target) 
{ 
	int i;
	for(i=0;i<(int) n_elements;i++){
		if(ptr[i] < target){
			return(i);
			break;
		}	
 
	}
return(-1);
 
 
}

Et je pense que cette fois ci, c'est la bonne !!!

As-tu réfléchi au sens du mot clef return et de l'utilité du break placé juste après ???

PS: pas besoin de parenthèses pour return. Ce n'est pas une fonction...

[Isabella83]

Bonjour,

As-tu réfléchi au sens du mot clef return et de l'utilité du break placé juste après ???

Si je ne me trompe pas :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
int find_index(double * ptr, size_t n_elements, double target) 
{ 
	int i;
	for(i=0;i<(int) n_elements;i++){
		if(ptr[i] < target){
			return(i); /* me retourne la premiere valeur de mon tableau pour laquelle j'ai la valeur < 0.05 (mon seuil) */
			break;  /* on sort de la boucle dès qu'on trouve une valeur <0.05 car je veux la première */
		}	
 
	}
return(-1);  /* si jamais aucune valeur n'est <0.05, alors la fonction retourne -1 */

Je peux ainsi vérifier dans mon main si la valeur retournée par find_index() est >0.

[ternel]

En fait, si, tu te trompes.
return provoque la sortie de la fonction.
Toute instruction placée après un return (dans le flux de controle) est inatteignable.
ton break n'est même pas atteint.

Par ailleurs toute fonction (non void) doit se terminer par l'exécution un return.

[Isabella83]

En fait, si, tu te trompes.
return provoque la sortie de la fonction.
Toute instruction placée après un return (dans le flux de controle) est inatteignable.
ton break n'est même pas atteint.

D'accord.
J'ai supprimé le break et effectivement cela me retourne pareil.
Merci pour cette explication

[Bktero]

Merci de cliquer sur le bouton en bas de page si vous n'avez plus de question.