Updated NLL calculation using dnorm()