devfun_vp function

Compute deviance of a linear mixed model as a function of variance parameters