for later:
need to check this with other reweighted least squares, optional convergence criteria
  RLM, GLM saves basic history, 
  in FGLS, I haven't done it yet