计算R中某点的函数导数

Question

在R中阅读this documentation，我不明白如何在特定点中获取该函数的衍生物。

They do it here in C（使用gsl_deriv_central / forward / backward），但我想知道R中是否有等价物？

Answer 1

只需安装包numDeriv并使用grad功能即可。以下是一些易于检查的简单示例。

library(numDeriv)

grad(sin, 1:3)
[1]  0.5403023 -0.4161468 -0.9899925
cos(1:3)
[1]  0.5403023 -0.4161468 -0.9899925


f = function(x) x^2 + 2*x +3
grad(f, 1:3)
[1] 4 6 8
2*(1:3) + 2
[1] 4 6 8