marginal_variance function

Marginal effect variance using the delta method