detecter des points spéciaux à partir d'un vecteur des données

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

Code:

---

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:

---

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 ? :oops:
Merci d'avance.

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:

---

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

---

Bien à vous

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

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

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

Alors,

Code:

---

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

---

diff() donne les différences:

Code:

---

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

---

Pour avoir les différences > que 70:

Code:

---

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:

---

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

---

Pour avoir la dernière valeur du vecteur:

Code:

---

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

---

Voici donc la solution final:

Code:

---

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

---

J'espère avoir pu vous aider!

Bien à vous,

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

Code:

---

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

---