If scale is None, then the default scale will be calculated.
Default scale is defined by self.scaletype and set in fit.
If scale is not None, then it is used as a fixed scale.
Returns:
scorendarray_1d
The first derivative of the loglikelihood function calculated as
the sum of score_obs