Discussion :

[orland]

Bonsoir
j'ai un vecteur de données comme suit :

Code :

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

31  41  70  80  94 133 400 420 442

je veux détecter les points : 133 , 442
J'explique comment : le principe est le domaine seuil c'est-à-dire le passage d'un point à un autre par exemple de 31 à 41 doit dépasser un seuil de 70.
donc de 31 à 41 la différence est < 70 je passe au suivant, entre 70 et 41 est < 70 .. ainsi de suite jusqu'à entre 400 et 133 où la différence est > 70 donc je retiens le point 133. De la meme manière je détecte le point 442.
j'ai essayé avec ce code mais ça me donne pas ce que je cherche:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
det=function(data)
{
	n = length(data)
      seuil= 70
	pt = vector()
 
	for(i in 1:(n-1))
     {
		ang = data[i+1]-data[i]
		if(ang>seuil)
              pt=c(pt,data[i+1])
	}
return(pt)
}

Avez-vous des idées pour m'aider ?
Merci d'avance.

[hoccha]

Hello,

Dans votre exemple, je ne comprends pas pourquoi le 442 est "détecté" (442-420=22 < 70).

Si non, une piste avec diff:

Code :

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

1
2
> diff(c(31,41,70,80 ,94,133,400,420,442))
[1]  10  29  10  14  39 267  20  22

Bien à vous

[orland]

concernant 442: je dois détecter toujours le dernier point (cas particulier).

[hoccha]

re-bonsoir,
Il est facile alors de trouver une solution avec le diff. Si vous ne trouvez pas, je peux vous aider.
Bien à vous

[orland]

Re-bonsoir
j'arrive pas à trouver la solution ! veuillez m'aider svp

[hoccha]

Alors,

Code :

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

v <- c(31,41,70,80 ,94,133,400,420,442)

diff() donne les différences:

Code :

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

1
2
> diff(v)
[1]  10  29  10  14  39 267  20  22

Pour avoir les différences > que 70:

Code :

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

1
2
> diff(v)>70
[1] FALSE FALSE FALSE FALSE FALSE  TRUE FALSE FALSE

Pour avoir les valeurs où la différence est > que 70:

Code :

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

1
2
> v[diff(v)>70]
[1] 133

Pour avoir la dernière valeur du vecteur:

Code :

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

1
2
> v[length(v)]
[1] 442

Voici donc la solution final:

Code :

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

1
2
> c(v[diff(v)>70],v[length(v)])
[1] 133 442

J'espère avoir pu vous aider!

Bien à vous,

[orland]

Merci
Mais j'ai réussi à le faire de cette manière aussi :

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
det=function(data)
{
	n = length(data)
      seuil= 70
	pt = vector()
      pt1 = vector() 
 
	for(i in 1:(n-1))
     {
		ang = data[i+1]-data[i]
		if(ang>seuil)
              pt=c(pt,data[i])
 
	}
pt=c(pt,det=function(data)
{
	n = length(data)
      seuil= data[n]/3
	pt = vector()
      pt1 = vector() 
 
	for(i in 1:(n-1))
     {
		ang = data[i+1]-data[i]
		if(ang>seuil)
              pt=c(pt,data[i])
 
	}
pt1=c(pt,data[n])
return(pt1)
}