Utilisation de numericDeriv

Bonjour,
Je cherche à comprendre l'utilisation de la fonction numericDeriv. si je tape

Code:

---

1 2	x<-as.numeric(1:5) numericDeriv(quote(x^3), "x")

---

J'obtiens

Code:

---

1 2 3 4 5 6 7 8

[1] 1 8 27 64 125 attr(,"gradient") [,1] [,2] [,3] [,4] [,5] [1,] 3 0 0 0 0 [2,] 0 12 0 0 0 [3,] 0 0 27 0 0 [4,] 0 0 0 48 0 [5,] 0 0 0 0 75

---

Or ce n'est pas ce que j'attendais (!) à la place de "12" (qui est la dérivée quand x=2) j'attendais 13 parce que [Y(t+1)-Y(t-1)]/[X(t+1)-X(t-1)]=(3^3-1^3)/[3-1]=13. et ainsi de suite (a la place de 27 j'attendais 28...)

Pourquoi donc ne calcule-t-il pas la dérivée numérique ?

merci pour votre aide.

nicolas

Bonjour et bienvenue,

Prendre un delta x de 2 pour l'estimation numérique d'une dérivée, c'est vraiment très grossier. La fonction se base sur la précision machine pour déterminer le delta à utiliser (tout ceci est expliqué dans writing R extensions).

Un petit exemple simplifié :

Code:

---

1 2 3 4 5 6 7 8 9 10

> eps <- sqrt(.Machine$double.eps) > eps [1] 1.490116e-08   > ((2 + eps)^3 - (2 - eps)^3)/((2 + eps) - (2 - eps)) [1] 12   > ((2 + eps)^3 - (2 - eps)^3)/(2 * eps) [1] 12 >

---

Cher pitipoisson

Merci pour votre réponse, je suis convaincu.


cordialement

nicolas