Take the numerical derivative for the first derivative of an ARMA using the 2 point rule.
jacobian_arma(theta, p, q, tau)
A vec
that contains all the parameter estimates.
A int
that indicates the number of AR coefficients
A int
that indicates the number of MA coefficients.
A vec
containing the scales e.g. \(2^{\tau}\)
A mat
that returns the first numerical derivative of the ARMA process.