Return negative loglikelihood function, i.e., - sum (log pdf(x, theta),axis=0) where theta are the parameters (including loc and scale)
