model_variance function

Compute model variance