#r #derivative
Вопрос:
У меня есть некоторые проблемы с выводом функций в 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)
но я беспокоюсь, что это не очень точно.
Спасибо!
Комментарии:
1. Вы видели этот пост RPubs?
2. @RuiBarradas интересно, хотя я думаю, что проблема здесь в том, что F(x) не является выражением, потому что это результат интеграции!
3. Вы правы, R не может вывести эту функцию. R-не лучший инструмент для символьных вычислений, разве вы не можете вычислить нужные вам производные вручную?