Discussion :

[nicolas_p356]

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

Code :

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

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

J'obtiens

Code :

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

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

[pitipoisson]

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 :

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

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
>

[nicolas_p356]

Cher pitipoisson

Merci pour votre réponse, je suis convaincu.


cordialement

nicolas