varmlm function

Residual variance of a model