Add option to fit GLM output rather than mean, conditioned on estimate output.