通读R 中的这个文档,我不明白如何在特定点取函数的导数。
他们在这里用 C 来做(使用 gsl_deriv_central/forward/backward),但我想知道 R 中是否有等价物?
只需安装软件包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
Run Code Online (Sandbox Code Playgroud)