Add function to calculate likelihood of saturated model

[e3bo]

This could be done using dtlmnet with a suitable design matrix. Initial work indicates that this will take a long time unless the speed of dtlmnet goes way up, and that the dev.ratio for our models with the flow rate predictors is not very close to one.