Я хочу найти частную производную f (x) по c, сигма, ню, с & сигма и с & ню. (т. е. df / d c, df / dsigma, df / dnu, df ^ 2 / dcdsigma, df ^ 2 / dcdnu)
f (x) имеет логарифмическое абсолютное значение и лог-модифицированную функцию Бесселя Третий вид. Я использовал функцию Deriv для вычисления частной производной f (x) по c, sigma, c & sigma. Я прочитал в документации пакета Deriv, что Drive может вычислять производную Бесселя по первому аргументу только по nu. Итак, как я могу рассчитать df / dnu и df ^ 2 / dcdnu?
fx=function(x,c,sigma1,nu1)
{
log (gamma(1/nu1))+log(abs(x-c)
+log ( besselK(x = ((1/sigma1^2) * abs(x- c) * sqrt(2 * sigma1^2/nu1)), nu = (1/nu1 ))))
}
dc=Deriv(fx,"c")
dcdc=Deriv(dc,"c")
dcdsig=Deriv(dc,"sigma1")
dsig=Deriv(fx,"sigma1")
dtheta=Deriv(fx,"theta1")