Add hrfs to glm()