如何将（一阶）导数作为函数写入R？

Question

我需要pt *（p ^ g）*（1-p）^ d）的导数作为p的函数。 t，d和g都是定义的。 我在努力：

firstder<-D(expression(p-t*(p^g)*(1-p)^d), "p")
firstderivative<-function(p){
firstder
}

然而，打电话

firstderivative(p=0.1)

我只给出了一阶导数的表达。 提前致谢！

Answer 1

firstder是类"call"的对象。

class(firstder)
#[1] "call"

你忘记评估电话了。

firstder <- D(expression(p-t*(p^g)*(1-p)^d), "p")

firstderivative <- function(p){
    eval(firstder)
}

g <- 1
d <- 1
t <- 1
firstderivative(p=0.1)
#[1] 0.2

Answer 2

使用mosaicCore包很容易。 鉴于您有F(x)=x^2并且您希望基于x计算导数，所以您有：

library(mosaicCore)


dx2x <- deriv(~ x^2, "x") 
x <- -1:2
eval(dx2x)

在你的情况下，它将是：

library(mosaicCore)

dx <- deriv(~expression(p-t*(p^g)*(1-p)^d), "p")

x <- -1:2
eval(dx)