r - R - 函数 F(x) 的导数

Question

我在 R 中导出函数时遇到了一些麻烦。给定函数 g(x,y)（我在这里给出 a 和 b 的简化示例）我定义

a <- Vectorize(function(x) -x^2/2)
b <- function(x,y) x*sin(y)
g(x,y) <- function(x,y) exp(a(x))*b(x,y)
F <-  Vectorize( function(y) {integrate(function(x) g(x, y), lower = 0, upper = 1, subdivisions=30000000)$value})

f <- Vectorize(function(x) deriv( function(x) CDF(x), namevec = "x") )

但它不起作用。我做了一个解决方法

h = 0.000001
f <- Vectorize(function(x) (F(x+h)-F(x))/h)

但我担心它不是很精确。谢谢！

标签： rderivative