Computes the penalized value of a derivative by adding the derivative of the penalty to it.
McCullagh_penalized(derivative, i1, j1, n, psi, delta, alpha, c = 1)derivative
the base derivative
first index of psi
second index of psi
matrix of observed counts
matrix of symmetry parameters
scalar asymmetry parameter
vector of asymmetry parameters
normalizing factor to make pi sum to 1.0. Default is 1.0.