Rで3階微分を求めたい

Question

Rで"-sum(a+bx)+sum(ylog(a+bx))"の3階微分を求めたいと思っています。numDerivパッケージで2階微分まで求めました。 次の3階微分まで求めたいのですが、方法はありますでしょうか。どなたかご存じでしたら、ご助言をお願いいたします。
なお、deriv()関数は「sum」関数があるとエラーがでてしまいます。

    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

Comments

[aokikenichi]

3階微分をD関数でして、その結果でsumとかでは駄目でしょうか