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