Я бы хотел вычислить третьи производные от "-sum (a + bx) + sum (ylog (a + bx))" в R. Я мог бы получить первую и вторую производные, используя пакет "numDeriv" как следующая программа. И еще одна производная, мне просто нужна третья производная. Как я должен сделать, чтобы получить это? Если вы знаете лучший способ, пожалуйста, дайте мне несколько советов.
library(numDeriv)
x<-c(1,2,3,4)
y<-c(2,3,5,4)
a<-1.278347
b<-0.888661
f<-function(a){return(-sum(a+b*x)+sum(y*log(a+b*x)))}
result<-genD(f, c(a))
result$D
#[1,]4.970284e-07 -1.239808