Discussion :

[Anna444]

Bonjour,

J'ai à traiter des résultats d'analyses réalisés à différentes températures (et qui doivent donc être corrigés par un facteur multiplicatif donné par le constructeur de l'appareil de prise de mesure).

Etant donné le nombre de mesures, j'ai voulu créer une fonction sur R, qui prend en argument une température et à partir d'un tableau fourni contenant les facteurs de corrections aux différentes températures retourne le bon facteur multiplicatif.
Seulement mon programme retourne des résultats plus qu'étranges. Pour certaines températures il retourne le bon facteur, pour d'autres, il retourne le facteur de la colonne d'à côté et je ne comprend vraiment pas d'où cela vient.
Peut être que je n'importe pas le tableau de facteurs de correction de la bonne façon ?
D'autant plus que si je demande de retourner l'élément contenu dans le tableau pour la même ligne et même colonne que celle retournée par la fonction mais depuis la console ou bien dans le script mais en dehors de la fonction là j'obtiens le bon résultat.

Mon tableau de facteurs de corrections se présente de la façon suivante :

La première colonne contient la partie entière de la température (donc varie de 3 à 28 pour 3°C à 28°C) les colonnes suivantes contiennent les facteurs corrections en fonction de la partie décimale de la température (précision au dixième), donc la ligne 1 colonne 2 contient la facteur de correction pour la température de 3,0°C, la ligne 1 colonne 6 le facteur de température pour 3,5°C ...



(en jaune sont surlignées 2 valeurs avec lesquelles le problème se pose, en effet pour la température 17,4°C la fonction retourne le facteur correspondant à 17,3°C, par contre pour 17,3°C il retourne le bon facteur)
Même problème avec 17,2°C, il retourne le facteur correspondant à 17,1°C. Étrangement le problème ne se pose pas pour 17,8°C ni pour 8,2°C ou 8,4°C
Le script est le suivant :




Et ci-dessous les retours que j'obtiens :





Du coup si quelqu'un d'entre vous a déjà rencontré ce type de problème ou bien a une idée d'où il pourrait venir cela m'aiderait beaucoup.
Merci d'avance et bonne soirée

[mgdondon]

Bonjour,

Il est préférable de copier-coller le code et de fournir des données pour le tester.

• Le problème vient du codage des nombres :

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
> T <- 17.3
> int_T=floor(T)
> dec_T=(T-int_T)*10
> dec_T
[1] 3
> as.character(dec_T)
[1] "3.00000000000001"
> 
> T <- 17.4
> int_T=floor(T)
> dec_T=(T-int_T)*10
> dec_T
[1] 4
> as.character(dec_T)
[1] "3.99999999999999"

• Votre programme fonctionne correctement en prenant l'arrondi : dec_T=round((T-int_T)*10).

• Une alternative pour accéder au tableau :

Code :

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

1
2
3
4
5
6
7
8
9
10
> Tab_corr_T <- data.frame(Int_T=17,C0=1.196,C1=1.193,C2=1.191,C3=1.188,C4=1.185,
+                                   C5=1.182,C6=1.179,C7=1.177,C8=1.174,C9=1.171)
> row.names(Tab_corr_T) <- Tab_corr_T$Int_T
> names(Tab_corr_T) <- c("Int_T","0","1","2","3","4","5","6","7","8","9")
> 
> fact_corr_T <- function(T,tableau_corrections=Tab_corr_T){
+   int_T <- floor(T)
+   dec_T  <- round((T-int_T)*10)
+   return(tableau_corrections[as.character(int_T),as.character(dec_T)])
+ }

Code :

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

1
2
3
4
5
6
7
8
> fact_corr_T(17.1)
[1] 1.193
> fact_corr_T(17.2)
[1] 1.191
> fact_corr_T(17.3)
[1] 1.188
> fact_corr_T(17.4)
[1] 1.185

Cordialement,

[Anna444]

Oh merci beaucoup pour cette réponse !

(Et désolée pour le copier-coller je n'ai pas l'habitude des forums de questions mais je note pour la prochaine fois )