R - 函数 F(x) 的导数

问题描述

我在 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)

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

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）