Compute Full Deviance
devfun_gls(varpar, gls_obj)
varpar
gls_obj
gls
the full deviance, a numerical scalar.
This code is adapted from code in devfun_vp internal function of pbkrtest package.
devfun_vp
pbkrtest
Useful links