examples, unit tests

do we get the same results as with HetGLS if family is gaussian?

essentially do WLS, or feasible WLS, with GLM.

(GLM has a lot of mix-and-estimate options, and I don't think we have tried most of the "weird" combinations yet.)