从 deriv 命令中提取梯度

Question

前几天，一位同事问了我以下问题。在下面的代码中，如何提取渐变：

> x=5
> a = eval(deriv(~ x^3, "x"))
> a
[1] 125
attr(,"gradient")
      x
[1,] 75

我的答案是

>  attr(a, "gradient")[1]
[1] 75

这种语法对我来说似乎很笨拙。有没有更好的方法来提取梯度？

Answer 1

不确定这些是否更好，但是：

with(attributes(a), gradient)

或

attributes(a)$gradient

是返回属性作为可供选择的列表的替代方案。

Answer 2

虽然它并不比您的方法更好，但您可以创建一个函数 grad，它接受带有渐变属性的数字并返回渐变值。

grad = function(x)attr(x,"gradient")[1]

grad(a)

现在可以重复使用。