%!demo
%! ## Compute the negative loglikelihood for a fitted Burr distribution
%! pd = makedist ("Burr", "alpha", 1, "c", 2, "k", 1)
%! rand ("seed", 21);
%! data = random (pd, 100, 1);
%! pd_fitted = fitdist (data, "Burr")
%! nlogL = negloglik (pd_fitted)
%!
%! ## This is useful for assessing the fit of a Burr distribution to data;
%! ## lower values indicate a better fit, e.g., for income modeling.
