Returns a function that computes the score (gradient of log-likelihood) with respect to parameters. Uses user-provided score function if available, otherwise falls back to numerical differentiation via numDeriv::grad.
Usage
# S3 method for class 'dfr_dist'
score(model, ...)