Check skewnormal residual deviance calculation

[nehill197]

We're wondering why there is a band that doesn't have as many values in the residual vs fitted values plot.

[nehill197]

Likely due to the cases where the deviance was slightly negative, which I set to 0.
See the dev_skewnorm() function.
Dig into this further to see why there are negative deviances.