Ecart entre deux courbes - Interpolation

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Aroniaina RANDRIAMANANJARA
Messages : 2
Enregistré le : 04 Juil 2019, 08:56

Ecart entre deux courbes - Interpolation

Messagepar Aroniaina RANDRIAMANANJARA » 04 Juil 2019, 12:47

Bonjour à tous

C'est un peu long désolée. J'ai un problème que je n'arrive pas à régler. J'ai une base de donnée contenant les colonnes :
Col
Paquet
Tub
NUM_BOITE
Zmmvert
Denshum
codevert
Zmmsec
Denssec
codesec

Mon but est de tracer la différence entre les courbes de Denshum et Denssec. Le problème est que la longueur de leurs axes des abcisses n'est pas la même (Zmmvert pour Denshum et Zmmsec pour Densec). Je les ai mis en graphe et il apparait qu'elles sont systématiquement décalées. Je voudrais alors étirer les valeurs de Zmmsec et faire une interpolation si cela est la solution pour avoir des nouvelles valeurs de Denssec et pour obtenir à la fin la différence entre les 2 courbes que je représenterait sur un graphe.
Ensuite à chaque fois que le contenu de la colonne "codevert" change (exemple de 0.0.6067 à 0.2.6067), je recommence la même chose, c'est à dire la différence entre les deux courbes (Denshum et Densec) + interpolation si c'est la solution pour que leurs axes des abcisses aient la même valeur de debut et fin.
Voici le tableau pour illustrer :

Code : Tout sélectionner

  ï..Col Paquet Tub NUM_BOITE Zmmvert Denshum codevert Zmmsec Denssec  codesec
1       0    354   0      6067       0 218,895 0.0.6067      0 336,571 0.0.6067
2       0    354   0      6067   0,625 572,143 0.0.6067  0,625 596,619 0.0.6067
3       0    354   0      6067    1,25 720,286 0.0.6067   1,25 620,238 0.0.6067
4       0    354   0      6067   1,875     720 0.0.6067  1,875 586,095 0.0.6067
5       0    354   0      6067     2,5 670,143 0.0.6067    2,5 526,286 0.0.6067
6       0    354   0      6067   3,125 609,714 0.0.6067  3,125 485,762 0.0.6067
7       0    354   0      6067    3,75 619,857 0.0.6067   3,75 481,571 0.0.6067
8       0    354   0      6067   4,375 551,952 0.0.6067  4,375 376,524 0.0.6067
9       0    354   0      6067       5 390,381 0.0.6067      5 276,524 0.0.6067
10      0    354   0      6067   5,625 334,571 0.0.6067                       
11      0    354   0      6067    6,25 343,524 0.0.6067                       
12      0    354   0      6067   6,875 501,476 0.0.6067                       
13      0    354   0      6067     7,5 634,857 0.0.6067                       
14      0    354   0      6067   8,125 558,667 0.0.6067                       
15      0    354   0      6067    8,75 501,571 0.0.6067                       
16      0    354   0      6067   9,375 478,238 0.0.6067                       
17      0    354   0      6067      10 484,952 0.0.6067                       
18      0    354   0      6067  10,625 463,095 0.0.6067                       
19      0    354   0      6067   11,25  420,19 0.0.6067                       
20      0    354   1      6067       0 209,952 0.1.6067      0   184,5 0.1.6067
21      0    354   1      6067   0,625 608,905 0.1.6067  0,625 507,952 0.1.6067
22      0    354   1      6067    1,25 911,619 0.1.6067   1,25 739,381 0.1.6067
23      0    354   1      6067   1,875  944,19 0.1.6067  1,875 739,762 0.1.6067
24      0    354   1      6067     2,5 908,619 0.1.6067    2,5 704,238 0.1.6067
25      0    354   1      6067   3,125 898,952 0.1.6067  3,125     716 0.1.6067
26      0    354   1      6067    3,75  901,81 0.1.6067                       
27      0    354   1      6067   4,375 897,667 0.1.6067                       
28      0    354   1      6067       5 895,095 0.1.6067                       
29      0    354   1      6067   5,625 898,381 0.1.6067                       
30      0    354   1      6067    6,25     896 0.1.6067                       
31      0    354   2      6067       0 250,143 0.2.6067      0 224,056 0.2.6067
32      0    354   2      6067   0,625     386 0.2.6067  0,625 315,476 0.2.6067
33      0    354   2      6067    1,25  516,19 0.2.6067   1,25     405 0.2.6067
34      0    354   2      6067   1,875 614,762 0.2.6067  1,875 464,619 0.2.6067
35      0    354   2      6067     2,5 673,286 0.2.6067    2,5 503,476 0.2.6067
36      0    354   2      6067   3,125 682,143 0.2.6067  3,125 553,286 0.2.6067
37      0    354   2      6067    3,75 791,667 0.2.6067   3,75 477,619 0.2.6067
38      0    354   2      6067   4,375 547,429 0.2.6067  4,375  306,81 0.2.6067
39      0    354   2      6067       5 360,048 0.2.6067      5 252,762 0.2.6067
40      0    354   2      6067   5,625 348,714 0.2.6067                       
41      0    354   2      6067    6,25  416,19 0.2.6067                       
42      0    354   2      6067   6,875  544,81 0.2.6067                       
43      0    354   2      6067     7,5 711,524 0.2.6067                       
44      0    354   2      6067   8,125 767,048 0.2.6067                       
45      0    354   2      6067    8,75 740,571 0.2.6067                       
46      0    354   2      6067   9,375 791,524 0.2.6067                       
47      0    354   2      6067      10 899,905 0.2.6067


Merci beaucoup si vous pouvez m'aider

Facundo Muñoz
Messages : 156
Enregistré le : 04 Juil 2019, 09:58
Contact :

Re: Ecart entre deux courbes - Interpolation

Messagepar Facundo Muñoz » 08 Juil 2019, 09:08

Bonjour Aroniaina,

Je ne comprends pas complètement ta question. Concrètement :

1. Quand tu dis "tracer la différence", tu veux dire faire un graphique de la courbe différence ?

2. Les valeurs des abcisses (Zmmvert et Zmmsec) sont toujours exactement les mêmes, sauf pour les observations manquantes du Denssec, ou les abscisses Zmmsec sont également manqantes. Serait-il possible d'utiliser une seule variable d'abcisses (disons, Zmm) égale à l'actuelle Zmmvert ?

3. En fait tu as la première moitié de la courbe Denssec, et tu voudrais la continuer, c'est ça ? C'est une extrapolation plutôt qu'une interpolation. Je crois pas qu'il y ait une méthode universelle de le faire. On est obligés d'utiliser des hypothèses et modéliser la courbe. Mais je peux pas te conseiller sans savoir de quoi il s'agit. C'est une série temporelle ? Y a t'il une forme paramétrique attendue (lineaire, quadratique, logaritmique, ... ) ? Peut-on faire des hypothèses sur l'ecart par rapport à Denshum (constant, paramétrique) ?

ƒacu.-
ƒacu.-

Aroniaina RANDRIAMANANJARA
Messages : 2
Enregistré le : 04 Juil 2019, 08:56

Re: Ecart entre deux courbes - Interpolation

Messagepar Aroniaina RANDRIAMANANJARA » 09 Juil 2019, 06:03

Bonjour,

Déjà je vous remercie beaucoup pour votre réponse.
En fait, Denshum c'est la densité d'une carotte de bois vert selon le plan radial. Zmmvert c'est la longueur de la carotte à l'état vert et tous les 0.625mm on a une valeur de densité, ce qui permet d'avoir un profil de densité le long de la carotte.
Et Densec c'est la même carotte mais après séchage, voilà pourquoi Zmmsec est plus court par rapport à Zmmvert car après le séchage, il y avait un retrait.
Mon but est d'avoir un profil d'humidité le long de la carotte en faisant la différence de la densité humide et la densité sec le long de la carotte, mais comme à l'état Sec, il y a eu un retrait, ils auront plus les même nombre de valeurs. Voilà pourquoi je me dis qu'il faut une interpolation ou extrapolation.

Et la colonne codevert et codesec sont les codes de la carotte, et quand le code change, on change d'espece. Et ce que je cherche à faire c'est avoir ce profil issu de la différence entre humide et sec de manière automatique à chaque fois que je change d'essence, c'est à dire à chaque fois que le code change

J'espère que j'ai pu éclaircir et encore merci.

Facundo Muñoz
Messages : 156
Enregistré le : 04 Juil 2019, 09:58
Contact :

Re: Ecart entre deux courbes - Interpolation

Messagepar Facundo Muñoz » 09 Juil 2019, 12:19

Bonjour Aroniaina,

Oui, c'est plus claire. Par contre ce n'est pas une question sur R, sinon sur comment analyser tes données. Ce qui est au-delà le propo de ce forum.
Je te suggère de réflechir et chercher de l'aide si besoin pour bien définir la manipulation que tu souhaites faire, et ensuite on pourra t'aider pour la mettre en place avec R.

Merci t à bientôt
ƒacu.-
ƒacu.-


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité