Discussion :

[jeff6868]

Bonjour,

Voici mon petit soucis. J'utilise actuellement une boucle dans une fonction afin d'identifier des valeurs comme fausses. Voici mon code:

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
out2NA <- function(x,seuil){
    st1 = NULL
    # Variable temporaire mémorisant la dernière valeur "correcte" relevée #
    temp <- st1[1] <- x[1]
    # Conserver l'indice de la valeur tampon #
    ind_temp <- 1
    # Définir ecart_temps max entre 2 comparaisons #
    ecart_temps <- 15
    # Conserver les temps de chaque élément de la série temporelle #
    tps <- time(x)   
 
    # Suppression des valeurs absolues dont l'écart avec la précédente est >= au seuil défini
    # et suppression des valeurs suivantes si la valeur aberrante se répète #
    for (i in 2:length(x)){
    if((!is.na(x[i])){
    if((tps[i]-tps[ind_temp] < ecart_temps) & (abs(x[i]-temp) > seuil)){  
    st1[i] <- NA
    }
    else {    
    temp <- st1[i] <- x[i]
    ind_temp <- i
    }
    }    
    }
 
    return(st1)
    }

Pour l'instant, il fonctionne bien ainsi. Il compare la valeur x[i] avec la valeur temporaire mémorisée, et regarde si l'écart temporel entre les 2 données n'est pas trop important. Mais comme il ne détecte pas tout à fait ce que je veux, j'aimerais lui rajouter une condition supplémentaire dans cette ligne:

Code :

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

1
2
if((tps[i]-tps[ind_temp] < ecart_temps) & (abs(x[i]-temp) > seuil)  
    &(abs(x[i+1]-x[i])<1)){

Il ferait donc la même chose qu'avant, mais avec en plus une comparaison entre la valeur x[i] et la valeur qui la suit juste derrière x[i+1].

En testant ainsi, j'ai la réponse: "valeur manquante là où true/false est requis".
Je rajoute donc dans cette ligne:

Code :

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

if((!is.na(x[i])&(!is.na(x[i+1])){

et si je teste à nouveau, j'ai l'erreur suivante: "le nombre d'objets à remplacer n'est pas multiple de la taille du remplacement "

J'ai donc essayé en changeant length(x) ou en utilisant dim(x), pensant que c'est le x[i+1] qui pose problème pour le calcul (et donc un problème de longueur du vecteur pour les comparaisons), mais je n'arrive pas pour l'instant.
Quelqu'un pourrait-il m'aider?
Merci d'avance!

[A. D.]

Bonjour,

Je n'ai pas trop le temps de me plonger dans votre code mais à première vue comme ça, le problème n'apparaît-il pas lorsque i=length(x) ? En effet, à ce stade vous allez faire un test sur x[i+1] qui n'existe pas.

Bonne continuation


Cordialement,

A.D.

[jeff6868]

Re-bonjour,

Merci à votre tour pour votre réponse. En effet c'est fort possible. J'ai effectivement essayé de jouer sur le paramètre length(x) dans:

Code :

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

for (i in 2:length(x)){

En essayant notamment de remplacer length(x) par: length(x)+1, ou bien de remplacer: 2:length(x) (en sachant que le nombre total de lignes est 35100) par: 2:35101, mais ça ne fonctionne pas. Et puis faire encore un: "for i+1 in..." est encore moins possible.
Comment faire dans ce cas justement pour que le x[i+1] existe et soit accepté?

[A. D.]

Bonjour,

Le fait de faire aller i jusqu'à length(x)+1 ne va pas résoudre le problème, bien au contraire, puisque vous aurez en plus une étape (lorsque i = length(x)+1 ) où ni x[i] ni x[i+1] n'existent...
Je ne suis pas trop rentrée dans votre algorithme, mais à quoi sert le test sur la valeur x[i+1]?

Une solution ici est d'arrêter la boucle à i = length(x) - 1, et de faire le dernier calcul (sur lequel vous n'effectuerez pas de test sur x[i+1]) à part (après la boucle).

Bonne continuation


Cordialement,

A.D.