计算R中某点的函数导数

Question

通过R中的这个文档阅读，我不明白如何在特定点获取函数的导数。

他们在这里用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